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ABSTRACT 


FACULTY OF SCIENCE 
CHEMISTRY 

Doctor of Philosophy 

ON THE MOLECULAR THEORIES OF LIQUID CRYSTALS 
by Claudio Zannoni 


This thesis is concerned with the statistical theory of liquid crystals 
and the molecular reorientation in these anisotropic fluids. 

The first chapter reviews the definition and the experimental determination 
of the order parameters. We then develop, in the second chapter, a 
statistical theory of the nematic mesophase, starting from the total 
partition function for rod-like particles and arriving at an effective 
orientational partition function. The usual molecular field theory of 
nematics is obtained as the first term in a rather general expansion. 

An extension of the molecular field theory to uniaxial nematics formed by 
non-cylindrically symmetric molecules is discussed in Chapter 3 and 
comparison is made with available experimental data. 

The study of molecular reorientation in such anisotropic systems as 
liquid crystals or membranes by using electron resonance is the subject 
of the remainder of the thesis. In Chapter 4 some background material on 
spin relaxation and the theory of reorientation in the presence of an 
ordering potential is reviewed. Methods for the computation of the 
spectral densities are also developed. The molecular reorientation of a 
monoradical spin probe in a highly ordered smectic A phase is studied in 




Chapter 5 by analyzing the angular dependence of the electron resonance 
linewidths in terms of both a strong collision and a diffusion model. 

The determination of the order parameter is discussed. The spin 
relaxation theory employed is generalized in Chapter 6 to include the 
effect on the angular dependent linewidths of the pseudosecular terms 
in the static spin hamiltonian. The predictions of this theory are 
illustrated with model calculations. Finally, in Chapter 7 a theory is 
developed, and experimental data analyzed, for the angular dependence of 
the spectral linewidths of biradical spin probes. 





CHAPTER ONE 


ORIENTATIONAL ORDER IN LIQUID CRYSTALS 


1.1 Liquid Crystals 

A liquid crystal is an equilibrium phase intermediate between 
a liquid and a crystal. At a microscopic level crystals formed of 
non-spherical molecules are characterized by a long-range order both 
in the position and the orientation of their constituents and normal 
liquids by the absence of both. In contrast a liquid-crystalline phase, 
often called a mesophase, can still possess an orientational order, 
while the positional one is greatly reduced or absent altogether (1,2). 
More precisely, the amount of positional order helps to classify 
liquid crystals into various types: nematics and cholesterics having 
no positional order, and smectics having the molecular centres of mass 
arranged in layers (3). An idealized diagram of the molecular organi¬ 
zation of nematics, cholesterics, and smectics is pictured in figure 1.1. 
The common feature in all these three types, as well as in the various 
subtypes, is the existence of a preferred or average molecular orienta¬ 
tion, called the director. 

Liquid-crystalline phases can be produced in two basic ways: 
either by dissolving the mesogen in an appropriate solvent or by a 
thermal transition. Substances giving mesophases in just the first 
way are called lyotropics (4), an example being soaps; the others 
are termed thermotropics (5,6). We shall be only concerned with the 
latter type in this thesis, even if many of our arguments, mutatis 
mutandis, could be applied to lyotropics as well. 


Figure 1.1 Molecular organisation in (a) the nematic phase 
(b) smectic phase, and (c) the cholesteric phas 





















A typical example of a thermotropic nematic liquid crystal is 
4,4'-dimethoxyazoxybenzene. This is crystalline at room temperature 
and gives on melting a turbid phase which on further heating yields an 
ordinary, clear, liquid. The mesophase owes its characteristic turbid 
appearance to the fact that the director changes continuously orienta¬ 
tion from one point of the sample to the other, to preserve the 
spherical symmetry of the isotropic phase, and the resulting orienta¬ 
tional inhomogeneity causes light scattering (1,2). This full rota¬ 
tional symmetry of the nematic mesophase is easily broken. For example 
the director can be readily oriented by applying a magnetic field (7), 
because of the enhancement of the anisotropy in the susceptibility 
caused by the orientational correlation. The result is that the motion 
of every single molecule with respect to the director is relatively 
unaffected by the field, but the liquid crystal as a whole is aligned 
(2,8). The oriented mesophase is an unusual state of matter in that 
it has,- as an example, the optical and dielectric properties of a 
uniaxial crystal, the director baing the axis of symmetry, while still 
retaining the mobility of a fluid. In a sense it is then an anisotro¬ 
pic fluid. Much of both the theoretical interest and of the practical 
applications of liquid crystals (9,10,11) stem from this fact together 
with their ease of alignment by magnetic and electric fields(12) as 
well as properly treated surfaces (13). 

From a chemical point of view the molecules of mesogens, i.e. 
substances giving rise to liquid crystalline phases, can be quite simple 
and it is estimated that, on average, one out of two hundred of organic 
molecules has this property. To give some examples we have listed in 
Table 1.1 the structure and the transition temperatures of some 


TABLE 1.1 The structure and the transition temperature of 


some nematogens. 


NAME AND FORMULA 


T ... /£ T . . /£ Ref. 

solid-nem. nem.-isotropic 


1) 4,4'-dimethoxyazoxybenzene (PAA) 


118 


136 


14 


M»0- 




r 

% 



-O Ml 


2) 4-methoxybenzylidene-4'-n-butylaniline 22 48 15 



4) 2,4-nonadienic acid 



23 


49 14 



/ 






Figure 1.2 The molecular organisation in (a) the smectic A phase, 

(b) untilted smectic B phase and (c) the smectic C phase. 
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TABLE 1.2 The transition temperatures of some smectogens 


(a) 


Name and formula 

Temperatures (*C) for the transition to 

Ref. 

S B S C S A N 1 

4,-4'-di(n-heptyl)azoxybenzene 


• 

34 

54 

71 

23 

Diethyl 4,4'-azoxybenzendicarboxylate 

' ' V M 

Vi 



114 


122.6 

6 

4,4'-di-n-heptyloxyazoxybenzene 

* < 7 M IS O 

*S N .J 


74 


94 

124 

6 

Ethyl 4-ethoxybenzylidene-4'-aminocinnamate 

C 1 M J 

o 

81.7 


118.5 

156.5 

159 

6 

Terephthal-bis (butylaniline) 

113 

144 

172 

200 

236 

1 


(a) S A , Sg, S c . indicate the smectic A, B, C types; N and I indicate the nematic and isotropic phase. 
The first number from the left of every entry is the transition temperature from the solid. 




























typical nematics. Common features are a high length to breadth ratio, 
rigidity, and anisotropy in the polarizability tensor, usually ensu¬ 
red by conjugated it electron systems. They often contain permanent 
dipoles but these are not essential for the formation of the mesophase, 
as proved by the existence of virtually apolar mesogens (e,g, the 
last entry in Table 1.1) and from the fact that liquid crystals do 
not exhibit the ferroelectric’.ty expected if dipole-dipole interactions 
were important. 

Nematogens, in particular, are not optically active. If the 
molecule is chiral, then the hypothetical mesophase produced will be 
cholesteric, i.e. a naturally twisted version of nematics having the 
director wound in a helical shape, with a pitch of typically about 
2000 A (1,16,17). 

Molecules giving smectics, i.e. layered liquid crystals, have 
again the same general characteristics but moreover they usually have 
flexible alkyl chains at the ends and strong dipole moments more or 
less perpendicular to the central part. The different ordering of 
these dipoles in the smectic layer has been invoked (18,19) as a 
possible explanation of the existence of many different smectic types 
(20,21,22). In figure 1.2 we illustrate three of these types, the 
smectic A, B and C. Notice that the preferred orientation is perpendi¬ 
cular to the layer in the A type and tilted with respect to it in the 
C type. The smectic B, in contrast with the other two types, possesses 
some degree of spatial order inside each layer. The director can be 
parallel or tilted to the layer normal in analogy with smectic A and 
C types. The structure and transition temperatures of some smectogens 


are listed in Table 1.2. 
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To describe a liquid crystal two concepts are particularly 
relevant: that of a director and of a set of order parameters. 

The director is a unit vector defined at every point of the 
sample and giving the local preferred orientation. It is assumed 
to vary slowly and continuously from one point to another except in 
certain regions called disclinations, where the preferred orientation 
changes abruptly (1,24). The concept of the director is a macroscopic 
one, even though it can be given a molecular statistical definition, 
and a continuum theory of liquid crystals has been developed (1,24- 
29) which describes the behaviour of nematics under external perturba¬ 
tions in terms of the deformation energy of the director field. Ty¬ 
pical examples are the determination of the equilibrium conformation 
of a nematic under certain boundary conditions or in the presence of 
an electric or magnetic field, the study of its flow and surface pro¬ 
perties etc... While successful in explaining such effects, the conti¬ 
nuum theories give no direct information on the behaviour of liquid 
crystals at a molecular level, which is our major concern in this thesis 

From a molecular point of view we shall find it particularly 
useful to introduce a set of orientational order parameters describing 
to what extent the molecules align with respect to the director. 

The importance of the concept of an order parameter in studying phase 
transitionsis well known (30). In general it can be defined as a quanti 
ty (or a set of quantities) which is non zero below a certain critical 
temperature T^» and zero above it. As an example, for the common 
liquid-gas transition the order parameter can be taken as the diffe¬ 
rence between the density of the gas and the liquid P G -p L and for a 
ferromagnetic system it can be taken as the magnetization, which 


again vanishes above a certain temperature, the Curie point. The way 
the order parameter changes in the course of a phase transition is also 
important since it tells us something about the character of the transi¬ 
tion; a discontinuity indicates a first order phase transition while 
a continuous change indicates a second order one. The first type is 
also characterized of course by a discontinuity in the first derivati¬ 
ve of an appropriate thermodynamical potential (the Gibbs free energy, 
say) and the other by a discontinuity or an infinity in the second 
order derivative of this potential. 

In the rest of this first chapter we shall treat in some de¬ 
tail the set of order parameters relevant for the description of li¬ 
quid crystals. In Section 1,2 this set is defined in a general way 
and then specialized to the important case of uniaxial mesophases 
composed of cylindrically symmetric molecules. The way the order 
parameters can be determined experimentally is then discussed in 


Section 1.3. 


1.2 The Orientational Order Parameters 


It is usually assumed that the orientational order in uniaxial 

liquid crystals is the average of just the second Iegendre polynomial: 

— ~2 

P 2 = (3cos 6 - l)/2, where B is the angle between the axis of symmetry 
of the molecule, assumed to be rodlike, and that of the mesophase. 

The upper bar denotes an ensemble average (31,32). The argument for 
this choice is that P^ obviously reduces to the proper limits: ~ ® 

in the case of no alignment, i.e. an isotropic system, and P^ = 1 in 
the case of a perfect alignment. However, Humphries, James and 
Luckhurst (33) have shown that in the molecular field limit one should 
more properly have not just one, but an infinite set of order parameters, 
given by the ensemble averages of all the Legendre polynomials: 

P 2* P 4’ P 6’*’' 

In this section we shall define a set of orientational order 
parameters in a general way i.e. as the minimal set of quantities 
necessary to describe completely the singlet orientational distribution 
function f(aBy). This distribution is defined as the probability 
density that the orientation of a molecule is given by the three 
Euler angles (aBy) (34). The ensemble average of any single particle 
orientational function A(aBy) can be expressed in terms of f(aBy) as 


A 


l 


A(aBy) f(aBy) da sinfi dB dy. 


( 1.1 


As any other well behaved function of three Euler angles the singlet 

distribution can be expanded in a basis of Wigner rotation matrices 

D m as 

m,n 





0.2 


f(agy) = l C L] 


D (apy). 
,mn m,n 


The definition and the main properties of the Wigner rotation 

matrices are described in Appendix 1. Multiplying both sides of 

equation (1.2) by D L and integrating over the angles, we find, 

n^n 


Lrrtn 


2L+1 

8 ^~ 


L* 


D 


m,n 


(1.3 


because of the orthogonality of Wigner matrices. The averages 

^ are our set of order parameters (35,36). The expansion can be 
simplified using the symmetry property of the mesophase or of the 
constituent particles. For uniaxial phases, as the nematic or smectic 
A are found to be, the singlet distribution must be invariant under ro¬ 
tation about the z axis if this is chosen along the director i.a. the 

• . . . . l* 

axis of symmetry. This implies that m must be zero in D 

m,n 

Another general property, for a uniaxial mesophase with D 

°°h 

symmetry, is that only even L can appear in equation (1.2). To simplify 
further one has to investigate the symmetry of the molecules forming 
the mesophase. We shall discuss this in some detail in Chapter 3. 

By now we consider the simplest case of all, which arises of course 
when the molecules are cylindrically symmetric. Then a rotation about 
their symmetry axis should not modify the distribution function, which 
implies n=0 in equation (1.3). In other words in this limit the distri¬ 
bution function will depend only on the angle g between the director 
and the molecular symmetry axis: 


£(S) ’ 1 C L D 0,0 (6) ' 1 C L P L <COS6) > 

L L 


( L even ) (1.4 


8 - 


where ” (2L+1) ?^l 2, and the P^(x) are the Legendre polynomials. 
If we let x ■ cosB then 


and 


f ( x) - \ + | P 2 P 2 (x) + | P 4 P 4 (x) + ... 
f f(x) dx ■ 1. 


(1.5 


For the most common case of cylindrically symmetric molecules the 
order parameters are just the averages P 2 , P^ etc. and it is clear 
that a knowledge of all the P completely defines the singlet distri- 
bution; while knowing f(x) all the order parameters can be calculated. 

An alternative way of characterising the singlet distribution 
is through its moments y : 


v 


n 




f(x) dx. 


( 1.6 


These are sometimes easier to calculate as we shall see in Chapter 2, 

where we give a compact expression for the moments of the distribution 
2 

f(x) « exp(£x ), which is frequently used in studies of liquid crystals. 
The representation of the distribution function in terms of order pa¬ 
rameters is clearly equivalent to that in terms of moments, since 


p L (x) 


r-0 



L-2r 


(1.7 


where [n] means the largest integer contained in n. (Since L is even 

[L/2] - L/2). The coefficients a T are defined as (37) 

L»r 



in our case 


<-> r 0) (vp- 



The first few example are 


1 


X 

X 


X 

(3x 2 - 

l)/2 

x' 

(5x 3 - 

3x)/2 

X 

(35x^ ■ 

- 30x 2 + 3)/8 

1 

X 


(P Q + 2P 2 )/3 

(3P X + 2P 3 )/5 

(7P q + 20P 2 + 8P 4 )/35 


(1.8 


Averaging eq.(1.7) or eqs.(1.8) gives the desired relation between 
moments and order parameters. 

We have defined before the order parameters as averages of 
Wigner rotation matrices. It is worth mentioning that often in the 
literature an equivalent formulation in terms of cartesian, as opposed 
to spherical, coordinates and tensors has been used to introduce parame¬ 
ters that measure partial alignment. As an example, considering the case 
of second rank order parameters, an ordering matrix S is obtained (32). 
The matrix S is symmetric and traceless so that in its proper frame it 

can be described by just two parameters: S and (S - S ). The 
J J zz xx yy 

explicit relation between the elements of the ordering matrix and 

~2 

the set of D is 
m,n 


zz 

D 0,0 



- S = /3/2 

tiJ 2 + ! 

XX 

yy 

0,2 

xz 

= -/378«Vl 

- 

yz 

' ^ 8(a o,l 

* *o.-r 

xy 

=i/5/8{D 0,2 " 

B 0,-2> 




1.3 Determination of the Order Parameters 


Let us suppose we want to measure a molecular 2nd rank tenso- 
rial property F during a time sufficiently long to allow the observation 
of statistically averaged quantities. If F^ ,m ^ are the irreducible 
spherical components of F, as defined in Appendix l, and if F is symme¬ 
tric, then its only non zero components will be those of rank L=0 and 
L=2. If the molecule is rigid and if we neglect local field problems, 
the L-0 component, being a pure scalar, will be unchanged going from 
an isotropic to an oriented fluid, and we shall not consider it. We 
then note that F is measured in a laboratory (unprimed) frame and that 
the L*2 components of F can be written in terms of the molecule fixed 
components F as: 


F ( 2 ,n) 


i F<2,m) ’ 

m 


( 1.10 


Taking an ensemble average we have, for a uniaxial mesophase. 


F (2,n> , £ d 2 F (2,m)' 

t 0,m n ,0 


m 


= ? (2 ’ 0) 6 


n,0’ 


( 1.11 


where we have implicitly taken the z axis of our laboratory system 
parallel to the director. 

If the order parameters are different from zero, as in the 
nematic mesophase, it then turns out that the molecular anisotropy of 
F is not completely averaged out by the molecular motion as in an 
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isotropic liquid. To stress this point the tensor F as given by 
equation (1.11) is often called a partially averaged tensor (38). 

To illustrate the argument let us imagine the following 
gedankenexperiment (39): we have a rodlike molecule with only two 
protons interacting via a certain dipolar tensor, this being our 
tensor F. The proton magnetic resonance spectrum of this (imaginary) 
molecule in an isotropic, low viscosity, liquid would contain just one 
line, but on passing to the nematic state two lines, separated by an 
amount proportional to the partially averaged dipolar coupling tensor, 
will appear. The splitting increases as we lower the temperature, 
corresponding to the rather obvious fact that in doing so the orienta¬ 
tional order of the system increases. 

Note also that the partially averaged tensor F has to reflect 

the synmetry of the mesophase. If this is uniaxial, as for nematics 

and smectics A, F has to be cylindrically symmetric, as shown by 

equation (1.11), irrespective of the symmetry of F'. The component 

p( 2 , 0 ) _ ^ 2/3 (p - p) that we measure is a linear combination of the 
II J. 

molecular components of F and of three order parameters: 


p(2,0> . 


„2 GM.)- 

0,0 


+ 2Re (D 2 x F (2,1), }+ 2Re{D 2 2 F (2,2) ) 


( 1.12 


Because of the structure of equation (1.12) it is clear that some 
simplification can be achieved when we work either in the principal 


coordinate system of the ordering tensor, i.e. the one in which Dq + j= 0 , 

(2 ± 1 ) ' (2 2 )' 

or in the principal system of F', where F ’ =0 and F ’ equals 

(2 - 2 )' 

F ’ and is real. It is usual to make one of these two assumptions. 



guessing a principal system from the molecular structure if the molecu¬ 
le does not have sufficient symmetry elements to determine its location 
a priori. In this case we find: 


^ <2 ’ 0) ‘ B 0,0 F <2,Q) * F < 2 ' 2 >< D 0,2 * D 0 ,- 2 ) 


(1.13 


l 2 F (2,0) + 2F^ 2,2 ^ d 2 n cos 2y , (1.14 


0,0 


0,2 


where 


D 0,0 - d 0,0 - P 2 ' and the fact that 


D L (cxBy) = d L (B) exp(-ima - iny) 
m,n m,n 

(c.f. Appendix 1) has been used. In particular, if F' has cylindrical 
(2 2 )' 

symmetry, F =0 and we obtain 


p(2,0) 


P„ F 


( 2 , 0 )' 


(1,15 


which gives immediately 

P 2 = (F II " V /(F |. " F P = AF ' AF '- (1 - 16 

We find the same result if the molecule is cylindrically symmetric, 

~2 

in the sense that D „ = 0. In these two cases the order parameter 

O, z 

is then obtained as a ratio of the anisotropy in the partially averaged 
quantity F to the molecular one. It is worth noting here that the order 
as given by equation (1.16) is a measure of the alignment of the 
z axis of the molecular system we have chosen to work in. As we mentio¬ 
ned before one has usually two practical choices: the principal system 
of the molecular tensor F or the principal system of the 


ordering tensor, this last one having the molecular long axis parallel 
to z. Since we are usually interested in the order of the long axis, 
it is useful to go back to equation (1.10) and suppose that the 
molecular primed system is always the long-axis system. We then effect 
one more rotation, of(a'6' y '),from this primed system to the principal 
system of the tensor F (double primed system) and then take the ensemble 
average. 


p(2,0) 


1 

m,p 


2 * 


D 


0 ,m 


D (a'BV) F 
m,p 


(2,p) " 


(1.17 


If F" has cylindrical symmetry F^ 2 ’^ = p(2,0) 

tensor is cylindrically symmetric as well then 



If the ordering 



D 0,0 6 Qm 


and we find 


AF/ AF" = F (2 ’ o) / f (2 ’ 0) " = D 2 >0 D 2 0 (B’), 


(1.18 


Even with these assumptions of cylindrical symmetry we see that we 
need to know the angle 8' between the long axis and the principal axis 


of F to be able to determine meaningful values of D = P . 

w jU z 

An extreme example is the case that the two axis are perpendi¬ 
cular i.e. B' = tt/ 2. Then AF/ AF"*-?^^ and the measured apparent 
order parameter i.e. the quantity on the left hand side is negative and 
goes from -1/2 to 0. This case is often encountered in electron reso¬ 
nance studies for disclike molecules like vanadyl acetylacetonate (40,41), 


1<4 ~ 


Another limiting case is that the two axis are at the "magic" 

2 

angle (i.e. cosg’=l//3) so that q(B')= 0 and no measurable effect 

- , . 13 

AF is emerging. As an example if F is a C chemical shift tensor the 

corresponding lines will be unshifted changing the temperature and 

even going to the isotropic phase (42). Of course if the principal 

system of F is along the long axis we recover the familiar expression 

af/af"=p 2 . 

A more complicated case arises when the ordering tensor is not 

cylindrically symmetric (43). We then have two independent order para- 

2 2 

meters to determine: D and Re (D } and at least two independent 

U jU VJ y L. 

measures are obviously needed. 

Additional complications can arise if the molecules considered 
are not rigid (44,45). More assumptions are then required to keep 
the problem manageable, such as assuming that the internal motions are 
decoupled from the reorientation of the molecule as a whole and that 
they are so fast on the experimental time scale that only their average 
effect need be considered. We shall not treat this problem in more 
detail since in this thesis we deal mainly with rigid molecules and we 
shall instead return, for the sake of clarity, to the simplest case 
of a rigid cylindrically symmetric molecule with a cylindrically symme¬ 
tric tensor F. 

If we consider again equation (1.16) we have that the quantity 
AF' is a molecular constant, relatively independent of temperature and 
other state variables, so that it should only scale P^- if we plot AF 
against temperature we will have the same shape of the plot of vs.T. 
However, if the absolute magnitude of the order parameter has to be 
determined, the problem is how to find the value of the molecular 
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quantity F' in the nematic phase because, even if the free molecule 
one is known, local field effects can alter it in an unpredictable manner. 
Up to now we have in fact been talking of a single molecule quantity F, 
while it is clear that what we measure is a macroscopic quantity 
corresponding to it and containing contributions from all molecules in 
the sample. Our discussion is thus fitting to the case in which the 
intermolecular couplings have a negligible effect on the quantity F as, 
for example, can happen in practice for the magnetic susceptibility X,(l). 
In this case the macroscopic tensor is simply 

(2,m) = p* ( 2 *®) 

X * P F X 

where p is the number density. 

In other cases it is clearly possible to measure a macroscopic 
quantity (dielectric constant, refractive index tensor etc...) but the 
relation of this to the corresponding molecular property is much more 
complicated and usually unknown. What is often done is to invoke some 
model for the local field or some empirical expressions. See as an 
example the use of Vuks (46) and Neugeubauer (47) equations by 
Chandrasekhar and his group when F is the electric polarizability (48). 

The experimental method mainly discussed in this thesis will be electron 
spin resonance spectroscopy (38). We shall deal with very dilute 
0v» 10 Si) solutions of stable free radicals (spin probes). In our case 
then we can safely rule out intermolecular or local field problems but 
we have of course an additional complication in that we are measuring 
properties (ordering parameters and correlation times) of a solute 
molecule and not of the liquid crystal itself. We shall return to this 



problem in Chapter 4. 


It is maybe worth while to point out here that most of the 

techniques measure only the order parameter P^. Apart from NMR 

experiments involving solute molecules with the aim of determining 

their structure (43,44), the assumption of cylindrical symmetry is 

usually made and equation (1.18) is used. In this case the apparent 

(2 2 ) 1 2 (2 0) 1 

order parameter measured is in error of an amount *v.|F^ 2 ^ * I 

We shall study the effect of deviations from cylindrical symmetry in more 

detail in Chapter 3. 

Even if the molecule is cylindrically symmetric it would be 

of some importance to determine the higher rank order parameters P.,P,etc. 

4 o 

However, as directly measurable quantities of rank higher than two are 
rather difficult to find, one has to turn to indirect methods. As an 
example a technique able to measure the mean square value of a 2nd rank 
molecular quantity will yield values of P^. It is in fact easy to 
prove that 


( J ( 2 , 0 ) / f ( 2 , 0)' ) 2 . 


<f„ - v 2 /(f,; - f; > 2 


(1.19 


I + ip *18 f 
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There are at present only two techniques used to measure P^. 
One is Raman scattering, used by Priestley and Pershan (50,51). The 
molecular quantity F' is in this case the differential polarizability 
tensor for a certain localized Raman active mode. For example a probe 
like N-(p'-butoxybenzylidene)-p-cyanoaniline is dissolved in a nematic 
and the strong line coming from the Ci N stretching mode is studied. 



The other technique consists in analyzing the electron resonance 
linewidth of a radical dissolved in a liquid crystal and has been 
developed by our group (52-55). It is worth pointing out that both 
techniques have been used to study solutes and not the liquid crystal 
directly and also that both of them rely on certain assumptions. 

In particular the first one is complicated by local field and multiple 
scattering problems (1). 

We shall discuss in detail the second technique in this thesis; 
for the moment we note that measurements on exactly comparable 
systems, which would allow a cross check of the two different techniques, 
do not exist as yet. This comparison would be particularly interesting 
in view of some surprising results obtained by using the Raman scatte¬ 
ring method. For instance negative values of P. have been obtained 

4 

(50,51), in contrast with the predictions of the current theories of 
the nematic phase (c.f. Chapter 2). 


CHAPTER TWO 


A STATISTICAL THEORY OF THE NEMATIC PHASE 


2.1 Introduction 

The task of developing a molecular theory for phases of aniso¬ 
tropic particles seems at first sight a formidable one. To quote 
de Gennes (1): "The statistical mechanics of liquids is difficult; 
the statistical mechanics of nematics is still worse". While this is 
certainly true in general, fortunately we can take liquid crystals 
also as an illustration of the so-called inventor’s paradox (56) i.e. 
that sometimes a more complicated problem is easier to solve than a 
simple one. 

To be more specific, we have to define what we mean by a molecu 
lar theory of liquid crystals. Let us first state that we are not 
really aiming at a theory capable of predicting, for example, all the 
various phase transitions and in particular the solid - nematic or the 
liquid-gas transition for an assembly of anisotropic particles. 

As is well known these problems are practically untractable even for 
systems composed of much simpler molecules. What we are seeking 
instead is a theory able of describing the most characteristic 
properties of liquid crystals. That is, first of all, that they 
undergo an orientational transition from a disordered, isotropic phase 
i.e. an ordinary liquid, to a phase where an orientational order para¬ 
meter such as P2 is different from zero. Apart from this qualitative 
prediction of the existence of such a phase transition, a theory of 
the nematic mesophase is expected to predict, or at least reproduce 




the behaviour and the values of the relevant thermodynamical quanti¬ 
ties at this transition. In particular the following experimental 
observations should be explained by the theory: the isotropic-nematic 
(I-N) transitions are of first order type with latent heats typically 

of 1 KJ mol ^ (32), small volume changes i.e. ^ 0.3% (57), small 

(K) 

entropy changes i.e. S /R 'v-O.S (6). The value of the order parameter 

at the I-N transition is usually found to be about 0.4, (32). 

Finally values of the higher rank order parameters such as P^, though 
difficult to obtain experimentally (c.f. Chapter 1), are very valuable 
in verifying theoretical predictions. 

Having said what we expect of a theory of liquid crystals, we 
shall proceed to make a brief survey of the existing theories. These 
have a relatively long history dating back e.g. to Bom in 1916 (58) 
and Zwetkoff in 1942 (31), although the first really successful attempt 
at developing a theory was made by Maier and Saupe (M-S) in 1960 (57). 
To have a liquid crystal it is quite clear that a necessary condition 
is a strongly anisotropic intermolecular potential. Maier and Saupe 
based their approach on the assumption that the anisotropic potential 
between two molecules, considered rigid and cylindrically symmetric, 
comes solely from London dispersion forces. This attractive intera¬ 
ction is then treated in the molecular field approximation, generali¬ 
zing the well known Weiss model of ferromagnetism (59). The spirit 
of the molecular field treatment is to replace the detailed intera¬ 
ction between molecules (spins in the case of magnetism) by an 
effective potential or pseudopotential representing the average 
interaction of all other molecules with the one in question. 


The M-S anisotropic pseudopotential obtained in this way has the form 


U(x) = (A P 2 / V 2 ) P 2 (x), (2.1 

where x is the cosine of the angle between the symmetry axis of the 
molecule and that of the mesophase, A. is a molecular constant, propor¬ 
tional to the square of (a |( -a^), the anisotropy in the polarizability 
tensor of the molecules and V is the volume of the system. The re¬ 
maining quantity, P 2> is just the normal order parameter defined in 
Chapter 1. It is determined from the consistency condition as 

P 2 * |p 2 (x) exp{-$U(x) }dx / jexp{-3U(x) }dx, (2.2 

where 3 = -1/kT. 

It is easy to see that P 2 =0 is always a solution of equation (2.2), 
corresponding to a disordered system. As we shall discuss later 
(for a mathematically equivalent model) there is another solution with 
P 2 ^0, corresponding to the nematic mesophase, which exists below a 
certain critical temperature T . The value of P_ at the transition 

K. 4 

is predicted to be 0.43 (57) in reasonably good agreement with experi¬ 
ment. Another prediction of the M-S theory is that P 2 should be a 

universal function of the reduced temperature T* = T/T , where T is 

K. K. 

the absolute nematic isotropic transition temperature. In figure 2.1 
we compare a plot of P 2 vs.T* according to M-S theory with some typi¬ 
cal experimental results. We see that the agreement between theory 
and experiment is generally good but also that there is not really a 
universal curve. Another prediction is that the entropy at the 



Figure 2.1 The orientational order plotted as a function of the 
reduced temperature for a number of nematogens (2). The 
continuous line is the prediction of the Maier-Saupe theory 
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transition S /R should be 0.42.Again this isin qualitative agreement 
with experiment. 

Although very successful and widely used, the M-S theory has 
been subject to some criticism. For instance, M-S theory does not 
appear to consider steric or indeed any type of intermolecular forces 
except dispersion. Kaplan and Drauglis (60) have noticed that an 
order of magnitude calculation based just on dispersion forces gives 
a far too low value for the interaction parameter A as compared to 
the value determined from the transition temperature T . Their calcu- 

lx 

lation is itself open to criticism since the actual values of the 
molecular polarizability anisotropy cannot be extracted from the re¬ 
fractive index data with any certainty because of local field compli¬ 
cations None the less there is some doubt the polarizability alone 
is responsible for the existence of liquid crystals. Other discrepan¬ 
cies, between the Maier-Saupe theory and experiment this time, come 
from the fact that the volume dependence of the transition temperatu¬ 
re, as measured by the value of (dlnT/dlnV)— is nearer to four (61) 

2 

than to the value of two predicted from equation(2.1). 

Apart from dispersion forces the other obvious=candidate as 
an orienting factor and as responsible for liquid crystal formation 
is a steric (repulsive) one in combination with the anisotropy in the 
molecular shape i.e. with the fact that nematogen molecules are rather 
elongated and similar to rigid rods. Many authors, starting from 
Onsager in 1942 (62) have shown that systems of elongated objects 
interacting with just a hard potential i.e. 



°° if there is overlap 
0 if there is no overlap 


(2.3 


can give a transition to an oriented phase when the density of 
elongated objects becomes sufficiently high. The steric theories (62,63), 
however, have not been extremely successful. Firstly it is clear that 
purely steric models give a phase transition independent of temperatu¬ 
re and so they are rather inappropriate for the treatment of thermotro¬ 
pic liquid crystals. Apart from this difficulty, they predict values 
of the order parameter at the transition far too high (-'0.8) and 
also the entropy change is predicted to be much larger than observed. 

A Monte Carlo calculation (64) has been performed on a system of 
hard spherocylinders (65) but failed to exhibit an orientational 
phase transition even after over 50*10^ configurations had been gene¬ 
rated. 

While then it seems clear that steric forces as well as disper¬ 
sion are relevant to describe interactions between nematogen molecules 
it is likely that a blend of the two and possibly of other factors 
such as deviations from cylindrical symmetry, flexibility, effects 
of electrostatic multipoles will have to be taken into account before 
we have a proper understanding of the relation between the nematic 
behaviour of a molecule and its structure. In view of all these com¬ 
plications it may now become difficult to see why the Maier-Saupe 
theory is so successful. The reason is probably a very simple one, 
namely that the M-S potential has the correct functional form. A simple 
heuristic argument to illustrate this point is that the potential 
for a rod in a (molecular) field of uniaxial symmetry will be a function 
of the angle 6 between rod and field symmetry axis only. This means 



that we can expand U(x), where x=cos 6, as a sum of Legendre polynomial 
and if we approximate this series with the first symnetry allowed term 
then U(x) , vcP 2 (x). We do not know anything about the coefficient c 
but on the other hand, since U(x) has to vanish in the isotropic phase, 
it is obvious that the simplest correct assumption is that c will be 
proportional to the order parameter i.e. 

U(x) * c' P 2 P 2 (x), (2.4 

which is essentially the Maier - Saupe result but where now c' is 
essentially an adjustable parameter. 

A general mean field approch has been developed by Humphries, 
James and Luckhurst (H-J-L) (33). Their treatment consists of writ.’ng 
down the intermolecular potential between two molecules as a Pople 
expansion i.e. expanding the potential in a complete product basis 
set of spherical harmonics of the orientations of the two molecules 
(66,67). Three successive averages are then taken to obtain the 
pseudopotential for a single molecule. The first is over all the 
possible orientations of the intermolecular vector in space, the 
second over all the orientations of the second molecule in the molecu¬ 
lar field and the third average over the intermolecular distances. 

Their result is 


U(x) = l u L P L P L (x), (L even) (2.5 

la 

where P are the orientational order parameters determined self- 
L 

consistently by a set of equations similar to equation (2.2); 


Uj is an average interaction coefficient, which is a function of mole¬ 
cular properties as well as of the volume of the system. Humphries, 
James and Luckhurst studied in particular the effect of retaining terms 
up to L=4 on the properties of the mesophase. They found that in many 
cases the quartic term is necessary to reproduce quantitatively the 
temperature dependence of P ^ and in particular its value at the 
transition (68). In the H-J-L theory the volume dependence of u is 

JL 

assumed to be u^ = u^V Y . A fitting of the temperature dependence 
of gi ves a value of y=4 (69), in excellent agreement with the 

independent measurement of ( dlnT/dlnV )- by Me Coll and Shih (61). 

*2 

Both in M-S and H-J-L theory the calculation starts from the 
pairwise intermolecular potential which is then averaged in a certain 
way to obtain the molecular field. While this seems a physically 
plausible procedure, it is, in a sense, a somewhat unconventional 
statistical calculation and a formal justification should be provided. 
It is clear in fact that in general one has just a total potential 
energy for the system and that the various averages have to come from 
a proper partition of this total potential. It would be more 
satisfactory to start from the partition function for a system of 
anisotropic particles and proceed by successive, well-defined appro¬ 
ximations to obtain an effective orientational distribution and 
partition function. In Section 2.2 we provide this more rigorous 


treatment. 



2.2 A Molecular Theory of the Nematic Mesophase 


The configurational partition function for a system of N 
interacting rigid particles with positions defined by a vector R and 
orientation described by a set of three Euler angles (a,3,y)=fi can be 
written as (70) 


Q n = (1/N!)J{dR}{dO) 


exp{~6U(Rj...R^, 


( 2.6 


The notation (dR) and {dfi} is an abbreviation for the volume elements 

3 3 

in configurational space i.e. {dR} * d Rj...d and {dft}=dfij... df^. 

As usual, 6-1/(kT) where k is the Boltzmann constant and T the 
temperature. The potential energy U(Rj...R^,^ ... fi N ) of the system 
is a function of both the positions and the orientations of all the 
molecules. In general these will be coupled in a complicated way and 
the problem will be practically insoluble but, on the other side, we 
expect that liquid crystals can be described effectively in terms of 
an orientational potential only. In this section we start from 
equation (2.6) and derive an effective orientational potential and 
partition function. This is then compared with the mean field 
results of Humphries-James-Luckhurst (33) and Maier-Saupe (57), 
which were obtained in a more ad hoc way. One of the aims will be 
to define the various approximations involved to obtain these results 
and to indicate possible generalizations and improvements. The spirit 
is similar to that adopted in a paper by Schultz (71), with the 
difference however that he essentially poses the problem and then turns 
to a specific model where the molecules can only assume discrete 
positions and orientations. Instead we shall aim to be as general as 


possible. 


z.u — 

As a starting point let us assume that the total potential can 
be separated into a part Ug depending on the position of the molecules 
only together with an anisotropic part depending on the distance 
between molecules as well as their orientation: 

U = U c + U A . (2.7 

S A 

This assumption allows us to rewrite the partition function as 
Q n = (1/N! ) jexp(~8Ug) exp(-gU A ){dR}{dn} 

- (l/N!) Je*p(-6U S > (dRlJtdti) V ldR 1 

exp (“BUgHdR} 

= . j{dfl} exp(-6U A ), (2.8 

S , . . . 

where is the partition function relating to the position of the 

molecules (their centres of mass say): 

= (l/N!) jdRj dR 2 ...dR N exp{-$U s (R 1 ...R N )} f (2.9 

and the upper bar in equation (2.8) indicates a positional average. 

We see then that even after separation of the potential in two parts 
we are not yet able to factorize the totalpartitionfunction into 
positional and anisotropic contributionsas we should like. To achieve 
this Schultz (71) proves that 



£- / 



We see then from eauations (2.8) and (2.13) that the inequality (2.10) 
is indeed verified, being the quantity (correction terms) a non nega¬ 
tive function. 

We would now like to retain, for the sake of simplicity, only 
the first non-zero cumulant, i.e. U^, which amounts essentially to a 




high temperature approximation. This does not seem reasonable at first 
sight, since is an energy referred to the whole system and so is 
proportional to the number of particles N. On the other hand it is 
clear that the validity of the expansion i.e. the structure of the 
equations that we get should be independent of the number of molecules 
in the system. The solution to this apparent paradox often appearing 

in this sort of calculation is simply that the second term in the se- 

~2 —2 2 

ries.(U^ - U^),is still proportional to N and not to N because of 

the finite range of the intermolecular potential (73). It then turns 

out the condition of applicability of the expansion is that U A /N, i.e. 

the single particle energy, is small compared to kT. For the moment 

we shall then retain only the first term in equation (2.13) (but others 

terms can be treated in the same way) and try to give in more detail. 

Let us first introduce a pairwise additivity assumption: 


u A = E U(i,j), 
i<j 


(2.14 


which allows U. to be written as: 
A 


U A 


l d 3 R. U(i,j) 


1<J 


J?,i d \ exp{ ~ gy R i-v } 


d^.-.d 3 !^ exp{- BU s (R 1 ...R n )} 


(2.15 


l fi\ Aj 


U(i , j ) f (i, j ) 


1<J 


(2.16 


The pair distribution f(i,j) gives the probability of finding 


the centres of mass of molecules i and j respectively at fL and 
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Since we have assumed U to he isotropic, f(i,j) will depend only on 

O 

o 

| - Rj|= R^j and we could write f(i,j) = g(R„)/V , where V is the 

volume of the system and g(R) the pair distribution function of the 

2 

centres of mass. For large separations f(i,j) goes as 1/V . Notice 
that the isotropy in f(i,j) comes directly from equation (2.7) without 
further assumptions or, viceversa, that within that approximation it 
is not possible to allow for deviations from spherical symmetry of 
g(R) even though such anisotropy exists in liquid crystals. 

Since the explicit forms of U(i,j) and f(i,j) in equation (2.16) 
are not known, to proceed further we shall try to extract as much 
information as possible from their general properties. We can write 
the anisotropic intermolecular potential U(i,j) in a (complete) product 
basis set of Wigner rotation matrices (39) in the orientations of the 
two molecules i.e. as a series of terms of increasing tensorial rank 
(66,67). Choosing a coordinate system, which we call the intermole¬ 
cular or IM system, with the z axis along the R.j = R^ - R^ vector, 
this has the form 


U(i,j) 


Iu T _ (R. .) D Ll (fi. _ ) D L2 

u LiL 2 iunin 2 ij m,ni l ~ _ m,n2 


<Vr..>’ 

ij 


(2,i: 


where (2. defines the orientation of molecule i in this system and 
1”K 

the prime on the summation indicates that the = 0, “ 0 term is 

excluded, since it constitutes part of the scalar potential. 

To avoid the inconvenience of every pair of molecules defining 
a different coordinate system we now rewrite equation (2,17) in a 
laboratory system using the relation (see Appendix 1) 


m,n 


(n i-R. 


) - I d; 


ij 


l* 

q,m 


(n R..-L> Yn (n i-L)’ 
IJ 


(2,18 



where £1 and £1. represent a rotation from a laboratory system 

IV • a “ *JLl 1 “L 

1J 

to the IM system and to a coordinate system fixed on molecule i 


respectively. This gives: 


U(i,j) - T'u. . (R..) D L 1 (B. T ) D L 2 (£1. .) 


(2.19 


* d Li *(«„ D^ 2 * (n D _ T ). 


q,,m R.,-L “q^-m' R..-I/ 


l l lj 


1J 


Substituting in equation (2.16) we obtain 


U A = V 


J V (L , (R.. 

.“. . , L fl.Lmn.n, lj 

kj L 1 ,L 2 ,n 1 ,n 2 1 12 12 


) f(R..) R?.dR..} 
ij ij 1J 




{[d Li * 1^2 * d o 

J q^*® q 2 » -m K ij -L 


} ( 2.20 


x d*; 1 (o. ) d*; 2 (a. ). 

^l ,n l 1 " ^2 ,n 2 ^ ^ 


Since f(R^) is isotropic, the integral over the two Wigner 


rotation matrices gives immediately 


q + m Att , * x 

1 ’ (2L +1) q^,~q 2 m * m L 1» L 2 ’ 


( 2.2 


and we can rewrite the positionally averaged anisotropic potential as 


(j). (2.22 


i r 

i<j L.n^ 


: y (-) q D 1 ' (i) D L 

Ln^^ L q -q,n 2 


Here we have defined 


'Ln^n2 = 


4 ttV 
2L+1 


l <->" f 


m 


LLmn^ 


(R) f(R) R dR, 


(2.23 


and written (ft. T ) as (i), since from now on all the orientations 
l-L 

are measured from the laboratory frame and no confusion should arise. 

From equation (2.22) we see that is a purely orientational 
potential: each molecule can interact, through its orientation, with 
any other molecule, with an average "strength of interaction" coeffi¬ 
cient c. coming from the smear over all the molecules of the 

Ln l n 2 

local interactions. Thus, even very short range interactions give 
rise to an orientational potential not dependent on distance i.e. 
of infinite range, while of course the coefficients c, are infi- 
nitesimal in magnitude. For instance, if we suppose the molecule can 
interact with only z neighbours and with an average potential e, the 


coefficients 


Cu, l n 2 


will be of order ze/2N. 


Substitution of the expression for in equation (2.11) shows 
that the problem of the evaluation of the orientational partition 
function has been reduced to calculating 

= jdft 1 dft 2 ...dft N exp{-0 l U(ft.,ftj)}, (2.24 

i<j 

where 


U(ft. ,ft.) 
i J 


I ■ <-> q 

» n 2 


c D L (i) D L (j). 

Ln^ q.nj -q,n 2 J 


(2.25 
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The similarity between this potential and that used to treat orientatio¬ 
nal and magnetic transitions on lattices is quite striking (74-76). 
However, a difference with the lattice treatments that is worth pointing 
out is that they usually restrict the orientational (or spin) intera¬ 
ctions to nearest neighbours only, i.e. they have a local interaction, 
while here, since we have averaged over the distances, only infinite 
range orientational interactions appear. 

We could now go on with the potential in equation (2.22) and 
keep the treatment quite general but, since this complicates the book¬ 
keeping without changing the physics of the problem, we shall restrict 
ourselves to the case of cylindrically symmetric molecules (i.e. n^ 3 ^ 3 ^) 
with second rank interactions only (i.e. L=2). Then the anisotropic 
potential will be 

" A ' I l <-> q c 2 ■to' 1 * <,,„«> < 2 ' 26 

i< j q 

■ '2 ll D q,0 (i > D q>> < 2 ' 26 ' 

q i<j 

= c 2 l P 2 (cosS Vj ) - 
i<j 

Here C 2 = c 2oo anc * t * ie angle in equation (2.26’) is the difference 

of the two rotations ft. and ft^. In equation (2.26*) we have rewritten 
the potential by using the spherical harmonics addition theorem in a 
form that puts in evidence that the potential is invariant under rota¬ 
tion of the laboratory system and of the long axis of the molecules. 

We now write the Wigner rotation matrices in equation (2.26) as 



JJ - 


D 2 (i) = D - a) (i), 
q,0 q 


(2.27 



parameter to be determined variationally later on. We shall show 
that this provides a way of decoupling the potential into a "molecu- 

l 

lar field" part, in agreement with the Humphries-Jdmes-Luckhurst one 
(33),and into a correction term. We shall also show that the separation 
parameter D is, in the molecular field limit, just an orientational 
average of a Wigner rotation matrix i.e. an order parameter. To begin 
with we substitute equation (2.27) in equation (2.26) and find 


U. 


+ U, 


(2.28 


where 


U = c 
0 c 2 


l ("> q 0 l Di q>0 (i) + D l D q, 0 (i) ] } ’ (2 - 29 


provided the number of particles is so large that N(N-I)^N , and 


U 1 - c 2 l <-> q l « q (i) ^_ q <i)- <2.30 

q i<j 

We shall take U Q , which is linear in the orientational functions, as 
our zeroth order or unperturbed potential and U^, which is quadratic 
in the same quantities, as a correction. Accordingly the orientational 
partition function can be written as: 



- J4 - 


= jexp{-BU A > {dfl} 

Jexp{-BU 0 ) exp{-gUj} {dfl} 

= |exp{-gU 0 } {dO} - (2.31 

J exp{-gU 0 } {dft} 

= $ N exp (-gup. 

where the upper bar indicates the average over the unperturbed state 
and where we have defined a zeroth order N-particle partition function 


$ N = jexp{-gU 0 } {dfl} . (2.32 

Let us now investigate more closely the parameter D. We first notice 
that to obtain an effective potential for a mesophase having cylindri¬ 
cal symmetry we have somehow to break the full rotational symmetry of 

the potential U.. A natural way of doing this is to include an external 

N 

field term of the proper symmetry as c J P~(cosg.) where g. is the 

i=l 1 1 

angle between the field direction, taken as the laboratory z axis, and 
the i-th molecule symmetry axis (75-77). Of course, since we are 
interested in the field-free system, we then let the coupling parameter 
C go to zero. If we assume that this limiting process does not inter¬ 
fere with our successive mathematical steps, we can take the limit at 
once i.e. in practice avoid writing down the field term, remembering 
however that the potential has now to conform to the broken symmetry. 

In our case this means to be invariant under rotation about the z axis 
and under reflection in a plane perpendicular to it. After symmetriza- 



tion the potential becomes 


" c 2 J 1)2 * c 2 N D I P 2 (cosf^), 
i=l 


(2.33 


This in turn gives 


w 


<t = exp{-gu } {dn} 


exp(3c 2 j D 2 }(4> N , 


(2.34 


where {<f>} is the N-th power of the single particle pseudo-partition 
function 


exp{-c 2 P 2 (cos6 i )} dcosB^. 


(2.35 


The value of D is now determined by minimizing the appropriate unpertur¬ 
bed thermodynamic potential which, since we suppose the volume is 
essentially constant, is the Helmoholtz free energy A. 

From the definition of A we have 


A = -kT ln<& 
A 0 lnv N 


(2.36 


- ^ c 2 N 2 D 2 - NkT ln<)>. 


(2.37 


A necessary condition for a minimum is that 


3D" ' - C 2 N2D * c 2 N2f> 2 ' °- 


(2.38 



Thus, to minimize the free energy, the separation parameter D must 
satisfy the consistency equation 

D = ? 2 = jp^ (x) exp{-(Bc 2 N P 2 <x)} dx , (2.39 

or, in other words, the best separation parameter is just the order 
parameter P 2 - To confirm that D=P 2 is a true minimum the second de¬ 
rivative of the free energy_ nas to be calculated and this 

\ 3D 2 /D=P 2 

will give the region of stability of the phase. We shall see an 
example in Section 2.4. 

Let us now consider the correction term exp(-f}U^). In lattice 
problems a method for successive systematic approximations to this 
type of integral has been developed (74,75). The physical idea is to 
consider progressively larger clustersof molecules in the molecular 
field produced by all the others. The unperturbed potential U Q 
corresponds from this point of view to a cluster of just one molecule 
i.e. to the mean field approximation. The other terms allow for the 
short range effects, since the interactions within a cluster are treated 
in detail. Since our orientational potential is of infinite range 
we do not expect the correction term to be important, by analogy with, 
for example, ferromagnetic problems, where it has been shown by Kac 
(78) that the mean field solution becomes exact when the range of the 
interaction goes to infinity. 

To check this point consider the first cumulant of i.e. 
approximate exp(-pUj) by exp(-BUj). If we rewrite as 

U 1 = I C 2 {N2l)2 “ 2 ^ l P 2 (i) + i £ ( ~ )qD qO (i)D -qO (i)} ’ (2-4 ° 

i*j q 


l 



we see immediately that U^O since at this level of approximation , 
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as we have shown before. One should then retain at least terms up to 8 
and since we have rejected terms of this type during the course of 
the derivation, to be consistent we shall content ourselves in retai¬ 
ning only the unperturbed term as an approximation to the potential U^. 

In summary, we have shown that starting from the total partition 
function for a system of N particles with anisotropic interactions 
it is possible to obtain an approximate orientational partition function 
Corresponding to this partition function we have an effective 
orientational potential U Q which, in the case of cylindrically symmetric 
molecules with second rank interactions only is 

1 2-? - N 

U 0 = “ I C 2 N P 2 + C 2 NP 2 I P 2 (cosB.), (2.41 

i=l 

where c ^ (or Nc 2 ) gives the average strength of the interaction, is 
an order parameter and is the angle between the symmetry axis of 
the molecule and that of the mesophase. 

It is maybe worth pointing out here that the decomposition 
scheme we have adopted to reduce the total partition function is not 
the only possible one. As an example we could have introduced the 
separation parameters before the integration over the intermolecular 
separations. In general it is difficult to compare different decompo¬ 
sition schemes because they usually do not involve just a different 
ordering of the same operations but different sorts of approximations 
or physical assumptions as well. In any case, at the first level of 
approximation which is the one explicitly treated in this Chapter,the 
molecular field result is obtained anyway. 
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We are now in a position to compare our results with those 

obtained by Humphries, James and Luckhurst (H-J-L) and by Maier and 

Saupe. To do this, we notice first that H-J-L have not only retained 

terms with L=2, but, at least formally, all the terms with even L in 

the potential. We can do the same, by having separation parameters 

which are L dependent i.e. D L » instead of D. The symmetry argument 

and the minimization of the free energy would then give D - P , 

L Li 

where L must be even if the mesophase has uniaxial symmetry. If we 

1/N 

consider now the single particle partition function and 

define U 25 C 2 N we see that we have essentially recovered H-J-L results 
i.e. the effective potential of a molecule in the field of all the 
others is 


U(x) = u L P L P L (x). (L even) (2.42 

L 

Restriction to the first term, L**2, as we did before gives a potential 
formally analogous to the M-S one. Having established the equivalence 
of the treatments at this level of approximation, we may underline 
some other points. In this treatment we have obtained H-J-L as a first 
order approximation in a rather general expansion. In particular, 
our partition function is the proper one in the sense that we can obtain 
all the thermodynamic quantities from it by using the standard formulae 
(70,73). This sheds light on a sort of paradox of the usual mean field 
theory, where one has sometimes to use what seem ad hoc relations. 

In fact, since the mean field treatment concentrates on the single 
molecule potential U(x), what is naturally evaluated is not the partition 
function but what we have already called the pseudopartition function 



i.e. $ and not 


$ 1 * <V 1/N = exp(6u 2 p 2 /2) <ti . 

Thus, as an example, the quantity -kT ln<fi will not be just the free 
energy, but a term or in general £ u l P i / 2 wil1 have to be 

added as is usually done in practice. b 


2.3 The Maier - Saupe Model 

The simplest effective potential that gives a reasonably good 
description of the liquid crystalline phase has the angular dependence 
of a second Legendre polynomial. We have shown in the previous section 
how this potential can be obtained in a rather general way. The 
potential is formally equivalent to the one originally obtained by 
Maier and Saupe (57) and it could also be obtained on purely heuristic 
grounds (c.f. Section 2.1) or by using the Humphries, James and 
Luckhurst method (33). In this section we treat in some detail the 
problem of calculating the partition function, the moments and the 
order parameters for this Maier-Saupe like model. Since this has 
already been studied, the justification for the following Section lies 
in the fact that many of the results, particularly on the order para¬ 
meters, will be used again in this thesis and that to obtain them we 
give a number of new and hopefully useful expressions. 

The single particle partition function for the model is 

♦ x " exp(Bu 2 l>2/2H, (2.43 

<p * 2exp(gu 2 P 2 /2)Z Q , (2.44 

where we have defined a (reduced) pseudopartition function as 

Z Q = jexp(£x 2 )dx, (2.45 

J 0 

and £* j* u 2 P 2 = ^ c 2^2 * s the i nteract ^ on parameter, which is positive 
if the pairwise potential favours parallel alignment. Even for this 
simple angular dependence the pseudopartition function is not an 
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elementary function of the interaction £. We can however express 
in terms of known functions. Rewriting as 


Z Q - Jexp(£x 2 )dx 
0 

- fexp(Ct)t~ 1/2 (1-t) 3 ' 2 * 1/2 -»dt. 
0 


(2.46 


we see that this is just the integral representation of the Rummer 
confluent hypergeometric function M(a,b,z) (79,80) if we take 
a*l/2, b=3/2, z*£. In other words 


Z Q - M(l/2, 3/2, C). 


(2.47 


As well as its compactness, this expression has the advantage of 

allowing us to use the large amount of results existing for these 

special functions, as we shall see in the following. Some of the 

most useful properties of the Rummer functions are listed in Appendix 2. 

Since, as we said, Z Q is not an elementary function of £ it is 

difficult to have a feeling for how it behaves. To this purpose a 

simple approximate expression for Z^ would be useful, and one is obtained 

in the following way. Notice first that a trivial power expansion 

2 

of Zq is obtained from equation (2.45) by expanding exp(£x ) and 
integrating term by term i.e. 

Z Q * l £ m /[(2m+l) • m!] * 1 + £/3 + £ 2 /10 +... (2.48 

m=0 
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is given in figure 2.2. 

A compact expression for all the even moments of the orientatio- 

2 

nal pseudodistribution f(x)=exp(£x )/Z Q can be obtained using the pro¬ 
perties of the Kummer functions. By definition (see Chapter 1) the 
moments are 



(2.54 


It is easy to see that 





M(l/2, 3/2, O, 


or using the differentiation formula for M(a,b,z), 


(2.55 


Z 2n = M(n+1/2 > n+3/2 f 5)/(2n+l). (2.56 

Substitution in equation (2.54) gives the desired expression for the 
moments. An expression for all the Legendre polynomial averages can 
be obtained at once from equations (1.7) and (2.54). 

Asymptotic expressions for high values of the interaction 
parameter can be obtained from the corresponding formulae for the Rummer 
functions reported in Appendix 2. At high order 

z 2n = l (“>“ <» + 1/2 - ni) m f*, U - »), (2.57 

m=0 

where (C) = C(C+1)...(C+n-1) and ( C) n = 

n u 


1 . 




Figure 2.2 The orientational order parameter P^Cx), calculated for a 
pseudopotential - l s plotted as a function of 

the interaction parameter £. 






In particular one obtains 


x 2 = 1 -1/5 - l/(25 2 ) +0(5 3 ), 

x* = 1 - 2/5 + 1/5 2 + 0(5 _3 ). 

Or, equivalently, 

P 2 = 1 - 3/(25) - 3/(45 2 ) + 0(5" 3 ), 
P 4 = 1 - 5/5 + 25/(45 2 ) + 0(5~ 3 ) . 


(2.58 


(2.59 


We shall return to the very high order limit in Chapter 5, where 
we show that the same asymptotic results for the order parameters 
can be obtained also for the more general potential U(x) = £ C P (x) 
Even retaining just the term in 5 3 the asymptotic expression works 
rather well, especially for P^, when5> 6. In Table 2.1 some values 
of 1*2 and P^ obtained in this way are compared with values obtained 
by numerical integration. 

Low order expansions for the parameters P can be obtained easi 

L 

from equation (2.48). We just give the result for 

P 2 = 25/15 + 45 2 /315 + 0(5 3 ) (2.60 


When the limiting expressions for high and low order are not appli¬ 
cable one has in general to resort to numerical calculation to 
determine the order parameters. Since it is sometimes necessary 
(e.g. in linewidth studies, c.f. Chapter 5 and 6) to generate a large 
number of these order parameters P^ even of relatively high rank 
(up to L=30 say) the problem of their efficient 


TABLE 2.1 The order parameters P ^ and calculated by numerical 

integration and from the asymptotic expansion equation(2.59) 
truncated at the first term 


c 

? 2 

(asympt.) 

P 4 

P^(asympt 

5 

0.646 

0.700 

0.287 

0 

6 

0.712 

0.750 

0.361 

0.167 

7 

0.759 

0.786 

0.425 

0.286 

8 

0.793 

0.812 

0.480 

0.375 

9 

0.819 

0.833 

0.526 

0.444 

10 

0.839 

0.864 

0.566 

0.500 

11 

0.855 

0.864 

0.600 

0.545 

12 

0.868 

0.875 

0.629 

0.583 

13 

0.879 

0.885 

0.654 

0.615 

14 

0.888 

0.893 

0.676 

0.643 

15 

0.896 

0.900 

0.695 

0.667 

16 

0.903 

0.906 

0.713 

0.687 

17 

0.909 

0.912 

0.728 

0.706 

18 

0.914 

0.917 

0.742 

0.722 

19 

0.919 

0.921 

0.755 

0.737 

20 

0.923 

0.925 

0.766 

0.750 

21 

0.927 

0.929 

0.776 

0.762 

22 

0.930 

0.932 

0.786 

0.773 

23 

0.933 

0.935 

0.795 

0.783 

24 

0.936 

0.937 

0.803 

0.792 

25 

0.939 

0.940 

0.810 

0.800 

30 

0.949 

0.950 

0.840 

0.833 

35 

0.957 

0.957 

0.862 

0.857 

40 

0.962 

0.962 

0.879 

0.875 
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computation poses. One possibility is to generate the moments x 1 " in 

a recursive way, which is very easy to do,and then use equation (1.7) 

to obtain the P from these. However this is not so convenient from 

the numerical point of view since it implies the evaluation of the 

(small) coefficients a in equation (1.7) which are ratios of large 

L, r 

factorials. Moreover it is more satisfactory to have a recursion 
relation for order parameters without intermediaries. This prompted 
us to derive such a recurrence relation as we shall now describe. 

As it is well known (37) 


L Pj(x) = (2L-l)x P^U) “ (L-l) P L _ 2 (x). (2.61 


Multiplication of both sides by the distribution function and 
integration gives 


L P L - (2L-1)XP L . 1 - (L-l) P L _ 2 . 

The problem is then to evaluate xP , where L' = L-l is odd. 

Li 

Integrating by parts, using the fact that P , (0)=0 for odd L', and 

Lt 

remembering the relation (37) 


d 

dx 


P L' (X) 


(2L , -1)P l ,_ 1 + (2L'-5)P l ,_ 3 + 


(2L'-9)P L ,_ 5 +...+P 0 , 


(L' odd) (2.62 


we finally obtain 


■l'-TET f**P<»' z o 


(I-2)/2 

l 

n=0 


f(2L-3-4n)+ 


2g(L-l) 

(2L-1) 



(L even) 


(2.63 
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This allows recurrent computation of every in terms of Z^. As an 
example we can write down the explicit expression for the first three 
order parameters i.e. 

? 2 -! e ‘Hf - (1 ’ < 2 - 64 

- <5 P 2 - 1 ’- < 2 - 65 

? 6 ‘ H? < - < 9 * TT^ ? 4 - 5? 2 - 1 ’’ <2 ‘ 66 

According to the Maier-Saupe model the knowledge of one order parameter 
permits the prediction of all the others, (e.g. by using equation (2.63). 
In figures 2.3, 2.4 the Maier-Saupe orientational order parameters P^ 
and are plotted one against the other for positive and negative 
values of the interaction parameter £. It may be noticed that P^ is 
always predicted to be positive even when P^ is negative (e.g., for a 
disclike solute molecule whose symmetry axis tends to be perpendicular 


to the director). 




Figure 2.3 The relation between the order parameters and 
to the Maier-Saupe model for positive P^. 


according 



to the Maier-Saupe model, for negative P^. 
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2.4 The Nematic - Isotropic Phase Transition 

In the previous section we have studied the order parameter 
defined by equation (2.52), as a function of the interaction parameter 
5. We have seen that P 2 (5) is a continuous and differentiable function 
that vanishes only when 5=0. It is apparent that we cannot just look 
at its graph (figure 2.2) to locate where the nematic-isotropic transi¬ 
tion occurs. What we should do instead is to remember (c.f. Chapter 1) 
that a first order transition occurs when the slope of the free energy 
of the system has a discontinuity. Since we calculate only the 
orientational contribution to the thermodynamical quantities we can 
equivalently say that the nematic-isotropic phase transition occurs 
when the free energy of the mesophase becomes higher than the free 
energy of the isotropic reference state, which is just taken as zero. 

In figure 2.5 we plot the orientational free energy of the Maier-Saupe 
like model (continuous line). It is now clear from this figure that 
a phase transition does indeed occur when 5^5^ * 2.903, since there the 
curve crosses the abscissa. Remembering that 5= ~u^?^/kT and noticing 
at this point is 0.43, we obtain that 

= (-u 2 /kT R ) = 4.541, 

where T^ is the transition temperature. The transition is of first 
order, since the slope of the free energy, when it crosses the abscissa, 
is not zero, and at that point the order parameter jumps discontinuously 
to zero. The peak in the free energy at 5=2.2 correspond physically 
to the maximum temperature, T', at which the mesophase can be thermo¬ 
dynamically stable, even if it is unfavoured as compared to the isotro¬ 
pic phase since T>T R . For T>T' the mesophase is absolutely unstable (81) 
T' represents the superheating limit of the mesophase. 


100A/RT 



Figure 2.5 The orientational molar free energy for the Maier-Saupe 

model. The continuous line is the free energy of the hypo¬ 
thetical mesophase and the dashed line that of the disordered 


reference state. 
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CHAPTER THREE 


A MOLECULAR FIELD THEORY FOR UNIAXIAL NEMATIC LIQUID 
CRYSTALS FORMED BY NON-CYLINDRICALLY SYMMETRIC MOLECULES 


3.1 Introduction 

An essential characteristic of compounds forming liquid crystals 
is the rod-like shape of their constituent molecules, with an attendant 
high length to breadth ratio. This common feature has prompted both 
experimentalist and theoreticians to suppose that the molecules possess 
cylindrical symmetry; an assumption which has many attractions. For 
example, the ordering matrix, often used to describe the partial 
alignment in a mesophase, contains just one independent element and 
this can be determined by a variety of techniques (43,49). The assumptioi 
of molecular cylindrical symmetry is appealing to a statistical mecha¬ 
nician because the pairwise anisotropic intermolecular potential requi¬ 
red in any calculation is particularly simple for such particles (66). 
However, the molecules forming liquid crystals are, in fact lath-like 
and certainly do not possess the high symmetry which is frequently 
assumed for them. 

Two principal components of the ordering matrix are therefore 
required to describe the orientational order of a uniaxial mesophase 
composed of lath-like molecules although most investigators are content 
to assume that one is sufficient. However Alben et al. (82) have esti¬ 
mated the ordering matrix for 4,4'-dimethoxyazoxybenzene and find the 
deviation of this matrix from cylindrical symmetry to be significant. 
Their analysis is open to question (83) and it is regrettable that 
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there have been no further attempts to determine the entire ordering 
matrix for real nematogens. Fortunately the importance of deviations 
from cylindrical symmetry may be inferred from unambiguous determina¬ 
tions of the ordering matrix for rodlike molecules, such as 4,4'-dichlo- 
robiphenyl (84), dissolved in a nematic mesophase. In addition these 
matrices are found to be comparable to that estimated by Alben et al. 
for a pure mesophase. 

In unusual contrast to the experimentalists, theoreticians have 

been more concerned with exploring the consequences of deviations from 

molecular cylindrical symmetry. Thus Frieser (85) has shown that a 

system composed of particles with a lower symmetry than D , is capable 

°°n 

of existing either as a uniaxial or a biaxial liquid crystal. The 
possible existence of a biaxial phase has been studied in some detail 
for a system of hard rectangular plates using a lattice model (86), 
the Landau approach (87) and the molecular field approximation (88). 
These calculations confirm the original prediction of the eventual 
transition to a biaxial phase although this transition has yet to be 
discovered. Indeed, all known nematics are uniaxial so presumably 
the mesophase freezes before the predicted transition to the biaxial 
phase can occur. None-the-less although deviations from molecular 
cylindrical symmetry have yet to be evidenced by the observation of 
a biaxial phase such deviations should influence the order parameters 
for the liquid crystal in the uniaxial phase. For example, the devia¬ 
tion of the ordering matrix from cylindrical symmetry is clearly deter¬ 
mined by the molecular symmetry and the element of the intermolecular 
potential. Straley (88) has reported calculations of the ordering 
matrix for an ensemble of hard rectangular particles. These calcula- 
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tions are of some interest but they may not be particularly realistic, 
especially as recent experiments have indicated that dispersion forces 
may make a dominant contribution to the anisotropic intermolecular 
potential (89). In this Chapter therefore we shall develop, within 
the molecular field approximation, a theory for non-cylindrically 
symmetric particles interacting via a completely general intermolecular 
potential. We then obtain a series expansion for the orientational 
pseudo-potential in the uniaxial phase. The expansion coefficients 
are subsequently determined by assuming that dispersion forces make 
the only contribution to the anisotropic intermolecular potential. 

We are then able to investigate the influence of deviations from 
molecular cylindrical symmetry on the nematic-isotropic transition 
and the various order parameters. Finally we compare these predictions 
with the data available for the nematogen 4,4'-dimethoxyazoxybenzene. 
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3.2 The orientational Pseudo-Potential 

As we have seen in the previous Chapter the simple procedure 
employed by Humphries, James and Luckhurst (33) to obtain the orienta¬ 
tional pseudo-potential yields the same results,at the first level 
of approximation, as a more rigorous approch based on the factoriza¬ 
tion of the total partition function. For the sake of simplicity we 
shall then use in this Chapter the H-J-L procedure, which amounts to 
taking three averages of the intermolecular potential. The first of 
these is over all orientations of the intermolecular vector, the 
second over all orientations of one particle and finally over the 
intermolecular separation. 

We begin with the intermolecular potential which, for particles 
of general shape, we expand in a product basis of Wigner rotation 
matrices (34) as 


12 


“ I u 


L 1 L 2 m 1 in 2 n 1 n 2 


(R 


12 


) D 


-L 2 


„ D m n <S-R ) * 

I 1 K 9 n2 t K 


(3.1 


where is the separation between molecules 1 and 2 (67). The 
orientation of molecule i in a coordinate system containing the inter 
molecular vector as the z axis is denoted by The energy 

is clearly invariant under rotation of the coordinate system about 
this z axis and so we restrict the summation to those terms with m^ 
equal to -m 2 (67): 


U 12 ^ U L 1 L 2 mn 1 n 2 (R 12 ) ( °1-R* D -m,n 2 ^2-R* * 


, L 2 


where u 


LjL 2 mn^n 2 " ; \L 2 m-mn^n 2 . 


(3.2 
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It is convenient to define the molecular orientation in terms of a 
common coordinate system and so we transform the potential to a 

laboratory frame. The choice of this coordinate system is determined 
by the symmetry of the liquid crystal phase and for a uniaxial mesopha- 
se we shall take the laboratory z-axis to be parallel to the symmetry 
axis of the mesophase i.e. to the director. The transformation of 
is carried out by performing the rotation from the intermolecular 
vector to the molecule in two steps using the relationship 


m,n IK 


l D L *(ft„ ) D L (ft. ) 

L q,m R-L q,n 1-L 


(3.3 


Here the subscript R-L denotes the rotation from the laboratory to 
the intermolecular frame while 1-L indicates that from the laboratory 
to the molecule coordinate system. The intermolecular potential 
becomes 


12 


K 


T mn „ 

^L2mn^n2 12 




( <W D 


n 2 (n ? T> D n 1 T> 

q 2 ,n 2 2 “ L q i ,m RL 


if* 

q 2 


,-m^R-L^ 


9 


(3.4 


and we may now take the average over the orientation of the inter¬ 
molecular vector. If the distribution function for the intermolecular 
vector is orientation independent (see the discussion in Section 2.2) 
then we may take advantage of the orthogonality of the rotation matrices 

to evaluate the ensemble average ) d^ 2 * (ft ) 3S 

q^.m'' R-L^ q2»“m' R-L^ 
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D n 1 * m < fl p-r> D « 2 “- m (ft R-i> = 6 - - 6 ’ * /(2L, + 1). (3.5 

q^ ,m K—L. q^, m K-L. 


.Lo* 


= /_nq2 +m 


V q 2 L 1* L 2 1 


The partially averaged potential may then be written as 


r-n q+m i i 

U 12 = 2 (2L+1) u LLmn 1 n 2 (R 12 ) ^.n^l-L* D -q,n 2 (fi 2-L ) 


(3.6 


The next stage is to average over the orientations adopted by 
particle 2 and so we require . Within the molecular field 

4 >^2 * 

approximation this average is taken to be independent of the orienta¬ 
tion of molecule 1 and is to be identified with the normal ensemble 

average. Since we are only concerned with a D , uniaxial mesophase the 
_ °°h 

average ^ will vanish unless q is zero and L is even. The poten¬ 
tial is now reduced to 



(-l) m 

(2L+1) 


Uy _ (R.») D_ 

LLmn^n 2 12 0,n 2 



(fli_ L ); 


(3.7 


where the prime denotes a summation restricted to even values of L 
greater than zero. The scalar contribution to the potential is ignored 
because we are only concerned with the orientational properties of 
the mesophase. The final average is over all intermolecular separations 
and gives the orientational pseudo-potential, for molecule 1, as 


v i-l 


, m 


(2L+1) “LLmn^ °0,n 2 Vn^l-L* 


(3.8 


The expansion coefficients are defined by 


U LLmn^n 2 


= / U LLmn 1 n 2 (R 12 ) expf-UgO^.. .R N >/kT} d 3 ^ d 3 R 2 


/exp{-U g (R 1 -R^) /kT} d^-d 3 ^ 


(3.9 
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where Ug(R 

ensemble (c.f. Chapter 2). The pseudo-potential is conveniently written as 

U(BY) - V d£ >p D^ q (Br), <3.10 


. R n ) is the scalar potential of the N molecules in the 


where 


"Lqp ’ l <-*>" 

m 

the Euler angles (f},y) define the orientation of the director in the 
molecular coordinate system. 

We have three checks on this form of the pseudo-potential. The 
first occurs for an ensemble of cylindrically symmetric molecules for 
then the summation in equation (3.2) for the pairwise intermolecular 
potential is restricted to terms for which both n^ and ^ are zero as 
we assumed in the previous chapter. With these restrictions the 
pseudo-potential, given in equation (3.10), reduces to 

U(g) = p u LQ0 P L P l (cos6), (3.11 

in complete agreement with the result derived by Humphries, James 
and Luckhurst (33). The second particular case is an infinitely 
dilute solution of molecules of arbitrary shape in a liquid of cylin¬ 
drically symmetric particles. The orientational pseudo-potential for 
the non-cylindrically symmetric molecule in such a system is obtained 
formally from equation (3.10) by setting p equal to zero and yielding 
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U(ST) ' I' u Lqo P L D 0,q< 8T)> 


(3.12 


in accord with the result obtained by Humphries and Luckhurst (90). 
Finally if the summation in equation (3.10) is restricted to terms 
with L equal to 2, then the pseudo-potential is identical in form to 
that proposed by Straley on the basis of symmetry arguments (88). 

We conclude this section by considering the symmetry properties 
of the coefficients u^ . ‘ For identical particles the permutation 
symmetry of the pairwise intermolecular potential ensures that (67) 


U T T mn n < R 1 9> = (-D^ 2 U_ . n (R.,), 

L^2niu^ 1^2 1^ L 2 1-* 


(3.13 


and so 


^qp ” U ^pq’ 


(3.14 


The pair potential must also be real and this requires that 


U T T mn n (R 1? ) = <~ 1 > ni+n2 U T T -m-n -n (R 19>*5 < 3 * 15 


consequently 


Lqp 


(-i) q * p 


u * 
L-q-p 


(3.16 
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3.3 Parameterisation 

One of our major objectives is to investigate the influence 
of deviations from molecular cylindrical symmetry on the various 
orientational order parameters for a uniaxial nematic mesophase. 

We therefore seek to minimise the number of variables in the calcula¬ 
tion while maintaining the essential physics. As a first approxima¬ 
tion we shall consider only those terms with L equal to 2 in the 
expansion of the pseudo-potential. We expect the neglect of terms 
higher than quadratic to be a realistic approximation since similar 
assumptions for cylindrically synmetric particles have led to results 
in reasonable agreement with experiment (32,33). We can reduce the 
number of expansion coefficients still further by appealing to some 
specific model for the interactions or by imposing symmetry restrictions 
on the molecules. For example if we follow Straley (88) and take the 
particles to be hard rectangular parallelopipeds then the coefficients 
U£qp are independent of the sign of either p or q and zero if either 
of these subscripts is odd. In fact 

u 200 = { _2b ( w2+l2 > " 2W(L 2 +B 2 ) + L(W 2 +B 2 ) + 8WBL}/3, (3.17 

U 220 = ( l2 -BW)(B-W)//6, (3.18 

and 

U 222 = -L (W-B) 2 /2, (3.19 


where L is thelength, B the breadth and W the width of the particle. 
This particular parameterisation is only approximate (88) and in 
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addition anisotropic repulsive forces are probably not dominant in 
determining the behaviour of real liquid crystals. We have therefore 
derived expressions for the u 2 qp f° r molecules interacting via dispersion 
forces and, as we show in Section 3.3(b) these coefficients again vanish 
if either p or q is odd and are also independent of the sign of p or q. 

It is possible to place similar restrictions on the expansion 
coefficients , without appealing to specific forms of the inter- 
molecular interactions, by using more formal arguments based on the 
molecular symmetry and its influence on the pair potential. For example 
if each molecule has a centre of symmetry then (67) 


“l T nrn n (Hi?) “ (-l) Ll_ni U (R ) 

L l L 2 mn l n 2 ^ J l Ij 2 m-n l n 2 


(—l) L 2 _n 2 


U L 1 L 2 mn 1 -n 2 (R 12 ) 


(3.20 


consequently 


( - 1 )^l + ^ 2-n l~ n 2 (R ). 

V * u L 1 L 2 m-n 1 -n 2 lK 12 ; » 


u = (-) q u = (-l) p u * (-l) q+p 
Lqp L-qp Lq-p 


Vo- 


q-p 


(3.21 


but for this to be compatible with equation (3.16) the coefficients 
must be real. Further, if the molecules also possess a plane of symmetry 
orthogonal to their z-axes then 


u, t „ n ( R n o) = <"l) Ll u. n (R.,) - (-l)" 2 u 

* J l L ’2 nin l n 2 L^L^m-n^^ “2 


L * " (R„) 


= (-l) Ll+Lz u 


l 1 L 2 m-n 1 -n 2 (R 12 } 


(3.22 



- 59 - 


and so 

u T = (~1) L u = u = u_ . (3.23 

MP T^q-p L-q-p L q p 

However, equations (3.21) and (3.23) can only be consistent if both 
p and q are restricted to even values; in addition the coefficients 
are independent of the sign of p and q. 

Given these general restrictions on the coefficients we 

may write the pseudo-potential as 

l %p<p t ^p HD o, q <6 Y)*'>o,- q (6’r))/ 

L,|p|,|q| 

(even) 

(1 * V (1 * V- < 3 - 24 

If we now limit the summation to those terms with L equal to 2 we have 


U(6y) t u 200 d 0,0 + 2u 220 d 02 * cos2y ^ d o,0^ e) 

+ < 2 “200 d 0,0" 4 “222 4,2 COs 2 Y ) d 0,2 (6>COS2 Y ’ < 3 ' 25 

where d 2 (g) is a reduced Wigner rotation matrix (34). The second 
m,n 

rank order parameters d^ ^ and d 2 2 cos 2Y occurring in this expression 
are directly related to the principal elements of the ordering matrix 
(c.f. equations 1.9). In particular 


S 

zz 



9 


(3.26 
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S - S = /6 d'* _ cos2y- 
xx yy 0,2 ' 


(3.27 


We see therefore that the order parameter d^ _ cos2y reflects the 

u, ^ 

deviation of the ordering matrix from cylindrical symmetry . Other 
order parameters have been introduced to measure the deviation of the 
ordering matrix from cylindrical symmetry although they only differ 
by simple numerical factors; thus Alben et al. (82) use the symbol D 
to denote /6 d^ ^ cos2y while Straley (88) has chosen U which is equi¬ 
valent to /8/3 d^ ^ cos2y . For calculations it is convenient to write 
the pseudo-potential as 


U(8y) = -kT (a dg o ( B ) + b d2, 2 ( 6 ) C os2 Y } , (3.28 


where 


3 ^ U 200 d 0,0 + 2u 220 d 0,2 cos 2 y ) / kT > 


(3.29 


b = " {2u 220 d 0,0 + 4u 222 d 0,2 COs 2 Y > 


(3.30 


The pseudo-partition function for this single particle potential is then 


2 x 

Z(a,b) = j | exp{a d^ Q (g) + b d* 2 (0) cos2 Y } singdBdy , 


0 0 


= 2 7T 


| I Q {(b d^ 2 (g)} exp{a d£ Q (g)} sinBdg, 


(3.31 
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where I (x) is an nth order modified Bessel function (80). The orienta- 
n 

tional molar potential energy U is just one half the total average 
energy U (33), 


U - -(RT/2) (a d 2 >() + b d? cos 2 y ), (3.32 

the molar entropy S is 


S = - R ( a d o o + b d 0 2 cos2 Y> + R ln Z(a,b), (3.33 


and so the orientational molar Helmholtz function is 


A = (RT/2)(a d* + b d 2 cos 2 y ) -RT In Z(a,b) (3.34 
m u, u u, z 


This free energy will be a minimum provided the consistency equations 
for the order parameters 

2 7r 

d 0 0 = [ | d 0 0^ GX Pta d 2 q (b) + b d 2 2 (B) cos2 Y )}sinBdBd Y /Z(a,b) 

Q 0 

= 2 it j d 2 Q (g) I Q { b d 2 2 (p)} exp{a d 2 Q (&)} sinBdB/Z (a,b), 


(3.35 


and 


d£ 0 cos2 Y 


0,2 


2* TT 

’ | d 2 2 (e) cos2 Y exp(a d 2 Q (6) + b d 2 2 (B) cos2y) 


U 0 


sinBd8d Y /Z(a,b), 


= 2tt f d 2 2 (6) Ij(b d 2 2 (6)} exp{a d 2 Q (B)} singdB/Z(a,b) , 

* n 


are satisfied. 


(3.36 
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The two consistency equations may be solved in a variety of 
ways. For example we might first fix the three variables u__ n /kT, 
u 22 c/kT and u 222 /kT and then evaluate the integrals in equations (3.35) 
and (3.36) for trial values of the order parameters d^ ^ and d^ ^cosly. 
The trial values would then be varied, using a non-linear least squares 
fitting routine (33), until the consistency conditions are satisfied. 

An alternative, and often more convenient approach, is to tabulate the 
order parameters and thermodynamic functions for a range of values of 

the variables a and b. The ratios u 22c/ u 200 and u 222^ u 200 t * le 
expansion coefficients may then be extracted from the variables with 
the aid of equations (3.29) and (3.30). This exercise is particularly 
straightforward when the particles interact via dispersion forces because 

the ratio u 222^ u 200 j ust ^ U 220^ U 200^ Z aS we sha11 now show * 

Section 3.3(b) Parameterisation for Dispersion Forces 

Here we derive expressions for the expansion coefficients 
u T _ (R,_) when the molecules interact via dispersion forces (91). 

Our starting point is the dispersion interaction energy for a pair of 
molecules in their electronic ground states (92), 


U 12 d 12 A a6 B o6’ 


(3.37 


where the tensor convention of summation over repeated Greek subscripts 
is adopted. The proportionality constant d^ is given by 



-e 1 e 2 /(e 1 + e 2 ). 


(3.38 


where e^ is an energy parameter often identified with an ionization 
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potential. The tensorsA and B are defined in terms of the polariza¬ 
bility tensorsa^ and c^, for the two molecules, by 



V B (1/R 12> 


'66 


(3.39 


and 


B = a. V V.(l/R 1? ). 
a6 1 y 6 12 

ay 


(3.40 


The scalar product in equation (3.37) may also be written in irreducible 
spherical tensor form (34) as 


U 12 = d 12 l A (L ’ m) B (L ’ m) *. (3.41 

L,m 

where L may only take values 0,1 and 2 because A and B are second 
rank tensors. The irreducible components A^ ,m ^ are constructed from 
the tensors T and a according to the prescription (91b) 


A (L,m) = J {(2L 1 *1)(2L 2 +1))* W(L^L^.;1L) C(L^Ljm-q ,q) 


L 1 L 2 


xT (Ll ’ q) a 2 (L2,m ‘ q) , (3.42 


where W(abcdjef) and C(abc;de) are Racah and Clebsch-Gordan 
coefficients respectively (34). The values of L- and q may be 
restricted by working in a coordinate system containing the intermo- 
lecular vector as z-axis; in this case the only non-vanishing 
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component is 


T (2 ’°) = /3/2 T « /6/R3 . 

zz 12 


The same restrictions hold for and so 


U 


_ y _ /d \ (Lj,m)* (i> 2 >®) 

12 A ? L L m (R 12 ) °1 “2 

L 1 L 2 


(3.A3 


m 


where 


C L L m (R 12 ) = 'T d 12 ^ 2L i +1)(2L 2 +1)}i I W(L 2 121 5 1L)W(21 L i; 1L > 

1 2 R 12 L 


x C(2L 1 L;Om)C(L 2 2L;mO) (3.AA 


Transforming the polarizability tensorsto their respective molecular 
coordinate system gives the dispersion energy as 

U 12 “ IH - I )""" 1 '" 2 “ l Ll ’ ni ) *'- 2 L2, " 2) *) 


D -», n2 (!! 2-R> 


(3.A5 


This result has the same form as the general expansion of the pairwise 
intermolecular potential given in equation (3.2) and allows us to 
identify the coefficients in that expansion as 


I^I^mn^ „2 


„ «.,) - (R ) „< L 2> n 2>*; 

n 0 12 L^L^m 12 1 2 


(3.46 
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of course and may only take values 0 or 2 because the 

polarizability a is a symmetric second rank tensor. Consequently the 

surviving coefficients u T are found to be 

Lqp 


- d, 2 a <2 -"> a (2 - p) /5 R* 2 , 


(3.47 


after evaluating the vector coupling coefficients. 
Equations (3.29) and (3.30) now become 


and 


where 



(3.48 


(3.49 


(3.50 


and X is the ratio u 220^ u 200‘ 

We may therefore fix the parameter X, which is a measure of the 
deviation from cylindrical symmetry, at the start of the calculation 
and from the results, tabulated simply as a function of a , obtain the 
behaviour of the ensemble as the temperature is varied. In the next 
section we shall describe the results of such calculations and compare 
these theoretical predictions with experiment. For this final section 
we shall denote the order d^ ^ by the more usual P ^ and the secondary 
parameter d^ ^ cos2y by ^. 
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3.4 Predictions and Reality 

The first stage in the calculation is the identification of the 

transition from the isotropic phase to the uniaxial nematic phase. 

This may be accomplished for a given molecular asymmetry by determining 

the value of a at which the orientational Helmholtz function vanishes, 

provided of course there is no volume change at the transition. The 

value, dp. is then related to the transition temperature T by 
K K 

a K = -u 20(/ kT K‘ (3.51 


The results of such calculations are listed in Table 3.1 for a range of 

X values including A=0 which corresponds to the Maier-Saupe limit for 

cylindrically symmetric particles. We can see immediately that increasing 

the deviation from cylindrical symmetry decreases a and so increases 

K. 

the nematic-isotropic transition temperature for constant u 2qo’ ^t ^ rst 
sight this trend is somewhat surprising although it is in accord with 
calculations based on a quite different model (88). To appreciate the 


Table 3.1 The order parameters and entropy at the nematic-isotropic 
transition. 


X 

“k 

P (K) 

2 

d 2(k) 

0,2 

P (K) 

4 

S<«/R 

0 

4.541 

0.429 

0 

0.120 

0.417 

0.1 

4.505 

0.408 

0.017 

0.109 

0.384 

0.2 

4.386 

0.341 

0.035 

0.075 

0.275 


0.3 


4.160 0.207 


0.041 


0.029 


0.112 
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origin of this effect it is important to realise that varying X while 

holding fixed corresponds to keeping the component a’ of the 

(JU zz 

traceless polarizability tensor constant, increasing a' by a given 

XX 

amount and decreasing by the same amount. Under these conditions 
the orientational energy of a pair of particles, with parallel principal 
co-ordinate systems,becomes more negative as the deviation from 
cylindrical symmetry increases. It is this stabilization of the parallel 
configuration which is responsible for the increased range of the 
nematic mesophase. 

The values of several order parameters, at the transition 

— (K) 

temperature, are also given in the table. The parameter is seen 

to decrease with increasing X, a behaviour which might appear to be 

incompatible with the greater tendency for parallel molecular alignment. 

However, in the limit that the components a' and a 1 of the traceless 

zz xx 

polarizability tensor become identical, as for a disc-like particle, 
it is the y-axis which will tend to be parallel to the director. 
Consequently the z-axis will tend to be orthogonal to the director and 
?2 will be small or even negative; the order parameter must there¬ 
fore decrease with increasing deviation from cylindrical symmetry so 
as to approach this limit. Similar arguments hold for the higher para- 
meter P^ as X is increased. In contrast the order parameter ^ 
which reflectsthe deviation from cylindrical symmetry, shows the 
expected increase with increasing X. Finally, the entropy change at 

the order-disorder transition, also listed in the table, decreases 

—(K) 

with increasing X, in accord with the decrease in P^ 

The temperature dependence of the order parameters and 
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for the long molecular axis is shown in figure 3.1 for several values 
of X. The results are plotted as a function of the reduced variable 
(a/ct^) 1 which is identical to the reduced temperature T/T provided 
the coefficients u 2 are themselves independent of temperature. The 
lines labelled A=0 correspond to the Maier-Saupe limit and the effect 
ofintroducing deviations from cylindrical symmetry is to lower these 
curves as well as changing their slopes. The order parameter is 
found to be positive for all reasonable values of A and so deviations 
from cylindrical symmetry cannot be invoked as an explanation for the 
negative values of determined for 4-n-butyloxybenzylidene-4'-cyano 
aniline dissolved in 4-methoxybenzylidene-4'-n-butylaniline (51). 


The secondary order parameter D Q 2 is plotted as a function of P 2 in 
figure 3.2 for various A values. These results exhibit an unusual 


behaviour for the order parameter 2 is observed to increase with 

increasing P 2 pass through a maximum and then decrease, going finally 

to zero when the z-axis is completely ordered. The temperature 

~2 

dependence of the secondary order parameter D is shown in figure 3.3 

u 

for some values of A. 

We now turn to a detailed comparison of these predictions with 
experiment. The most direct evidence for the influence of deviations 
from molecular cylindrical symmetry on the behaviour of the mesophase is 


the observation of the secondary order parameter „ because this 
is identically zero for cylindrically symmetric particles. Of course 
the other order parameters will also depart from their Maier-Saupe 
values but similar deviations may also be observed for cylindrically 
symmetric particles if quadratic and quartic terras are retained in the 
pseudo-potential given by equation (3.11) (33,93). Unfortunately the 









8 



Figure 3.3 The temperature dependence of the secondary order parameter 

~ 2 ~ 

_ = (S - S )//6. The curves are calculated for 
0,2 xx yy 

A equal to (a)0.1, (b)0.2 , (c)0.3 . 
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only attempt to determine D _ for a pure nematogen was when Alben et al. 

u, / 

(82) showed how the nuclear magnetic resonance experiments performed by 
Rowell et al. (94) on partially deuterated samples of 4,4'-dimethoxy- 


azoxybenzene might be interpreted to yield D Q ^ as well as Although 

these experiments were performed throughout the nematic range the ana¬ 
lysis was restricted to a single temperature (82). We have therefore 
applied an analogous procedure to all of the data obtained by Rowell 


et al. (94) and found the values of ? 2 and D Q 2 plotted in figure 3.2. 
The results are reasonably well accounted for by the curve obtained 


with X equal to 0.2, especially as the values of D are particularly 

0,2 

sensitive to the fine details of the geometry assumed for the nematogen. 

Given a X of 0.2 we predict the entropy of transition S 7 /R to be 0.27 

in surprisingly good agreement with the observed value of 0.28, although 

be 

this does appear to subject to some experimental uncertainty (81). 

This prediction is a considerable improvement on the Maier-Saupe value 

for S /R of 0.42. The order parameter P 2 at the transition point 

may also be calculated for this particular X and we find 0.34 in good 
. the 

agreement with observed value of about 0.36. This prediction should 
be contrasted with the Maier-Saupe value of 0.43. The higher parame¬ 
ter is predicted to be 0.08 at the order-disorder transition but 

as yet there are no determinations of this quantity although it is 
available from light scattering experiments (50). 

As a final test of the theory we compare the predicted temperature 
dependence of the order parameter P 2 with the observed dependence. 

This comparison is not straightforward because the measurements are 
almost always made at constant pressure and as a consequence the volume 
will also change as the temperature is varied. The pair distribution 
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function employed in evaluating the average will then be temperature 
dependent and so the coefficients u 2 qp occurring in the theory will 
change with temperature. It is possible to allow for this variation 
by assuming a simple volume dependence of the form 


u„ = u§ /V Y , 
2qp 2qp 


where the u^ i s temperature independent (33). However, even though 
the temperature dependence of the volume is known, such an assumption 
introduces an additional variable y into the calculation. This pro¬ 
blem can be avoided if the order parameter is measured at constant 
volume for then the coefficients u^^ should, to a good approximation, 
be independent of temperature. Such measurements are available for 
4,4'-dimethoxyazoxybenzene (61) and the order parameter is plotted 
as a function of the reduced temperature in figure 3.4. The solid 
line also shown in the figure was obtained from the theory with \ set 
equal to 0.2. The agreement with experiment is good although the slope 
of the theoretical curve is slightly greater than that observed. 

The theory based on deviations from cylindrical symmetry and 
treated within the molecular field approximation would appear to provi¬ 
de an acceptable account of the orientational properties of the nema- 
togen 4,4'-dimethoxyazoxybenzene. However it is worth to recall that 
the same properties have been interpreted, with comparable success, 
by the pseudo-potential, for cylindrically symnetric particles, given 
by equation(3.11)with L restricted to 2 and 4 (33). The order parame- 
ter is predicted to be 0.3, the entropy change S '/R at the 

transition is calculated to be 0.32 and the temperature dependence of 



100 t/t k 


Figure 3.4 The temperature dependence of observed for 4,4'-dimethoxy 
azoxybenzene at constant volume; the curve is calculated wit 
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the order parameter P^ at constant volume is in complete accord with 
experiment (69). Of course the theory fails to predict the order 


parameter D 7 since for cylindrically symmetric particles this must 
u, z 

be zero. On present evidence it is hard to distinguish between the 

— (K) 

two theories although it may be significant that the parameter 

is usually less than the Maier-Saupe limit (69) which is a natural 

result of deviations from molecular cylindrical symmetry. Clearly 

future experimental studies of the orientational order in nematics 

~2 

should concentrate on the secondary order parameter D In the 

event that this is shown, unambiguously, to be non-zero, then the 
appropriate terms must be included in the pseudo-potential. However, 
such observations may not exclude the retention of the quartic terms 
in the expansion. Indeed when P^ exceeds the Maier-Saupe value of 
0.43, as it does for 4,4'-diethoxyazoxybenzene (68), then quartic terms 
will have to be retained because such an effect cannot be accounted 
for simply in terms of the pseudo-potential for non-cylindrically 
symmetric molecules when this is restricted to quadratic terms. We 
shall not attempt calculations with the extended potential here because 
the available experimental evidence does not warrant such an extension. 
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CHAPTER FOUR 


SPIN RELAXATION IN LIQUID CRYSTALS 


4.1 Introduction 

In this Chapter and in the following ones some static and dyna¬ 
mic aspects of liquid crystals are investigated using electron reso¬ 
nance spectroscopy. This technique can in fact often be used to 
investigate a diamagnetic system simply by adding a radical, or spin 
probe, whose interaction with its environment is reflected in the 
electron resonance spectrum. The nature of the system should, of 
course, play a dominant role in the choice of spin probe. For example 
if the system to be studied is a liquid crystal, or a membrane, then 
two properties of the spin probe are of prime importance. The molecules 
must deviate significantly from spherical symmetry so that they may be 
highly oriented by their anisotropic environment. Secondly there must 
be considerable anisotropy in the magnetic interactions in order for 
the extent of the partial alignment to be readily discernable in the 
electron resonance spectrum (38). These anisotropic interactions may 
also influence the widths of the lines as well as their positions because 
the anisotropy coupled to the molecular reorientation invariably 
constitutes a powerful spin relaxation process (38). Consequently the 
spectrum of the spin probe may often exhibit a characteristic linewidth 
variation and, because of the anisotropic nature of the environment, 
this variation will depend on the sample orientation with respect to 
the magnetic field. An understanding of such linewidth effects is 
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important not only as an aid to a correct spectral analysis but also 

because the linewidths contain valuable information relating to both 

the dynamic and static behaviour of the system (38). The static information 

obtainable includes the orientational order parameters, as 

and P^. Their values and their temperature dependence provide a test 

for the statistical theories of the mesophase (c.f. the previous Chapters). 

Several theories have been advanced to describe molecular reo¬ 
rientation in the presence of an ordering potential and they range 
from those in which the motion proceeds via small jumps (95-97) to 
the strong collision model where the jump angle can take any value 
(38,AO,52). 

In the next Chapters both these models will be used to interpret 
some experimental results and some pieces of theory necessary to connect 
the microscopic process to the observable spectrum will be developed. 

Our aim will be twofold. Firstly to obtain a qualitative description 
of the reorientation and try to distinguish between the models. 

Secondly to obtain values for the physically significant parameters of 
the models, such as the correlation times for reorientation of the long 
axis and about the long axis in the strong collision model, or the 

parallel and perpendicular components of the rotational diffusion 
tensor in the diffusion model. We shall treat in some detail the 

problems encountered in this type of analysis later on. In this 
Chapter we shall instead discuss some background material necessary 
for the successive developments. 
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4.2 Spin Relaxation Theory 

This section gives a brief summary of the theory of spin rela¬ 
xation originally due to Bloch (99), Redfield (100) etc..(101). To 
start with, it is useful to recall that the state of any system, 
classical or quantum, can be described in terms of an operator called 
the density operator or the density matrix of that system (102). In 
fact any observable A(t) can be calculated from the density operator 
p(t) at a given time t as (103) 


A(t) = Tr{p(t)A}, 


(4.1 


where Tr{...} indicates an ordinary trace or an integral over the phase 
space as appropriate. The time evolution of p is governed by the 
Liouville equation 

p(t) =frp(t) =~i[H,p(t)] 


(4.2 


= -iH X p (t). 


Here H is the hamiltonian of the system; the square brackets denote a 

commutator in the quantum case (ft = 1 is used here) and /-I = i times 

. . x 

the Poisson brackets in the classical case. The operator H acts on 

operators producing their commutator with the hamiltonian. Operators 

acting on operators to give other operators are called superoperators 

x 

and H is often called the Liouville superoperator. The algebra of 
superoperators is treated in detail elsewhere e.g. by Muus (104); here 
W e just mention two properties which will be used later, firstly that 

x 

the eigenvalues of H are energy differences (as directly observed in 
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spectroscopy) and secondly that 

exp(H) A exp(-H) = exp(H X ) A. (4.3 

In a magnetic resonance experiment (105) the observable can 
be related to a component of the complex magnetization at a certain 
frequency (or time). This in turn can be related to the expectation 
value of a particular spin operator S^, where a=x,y,z . As a simple 
example. 


M a (t) = g6Tr{S a p(t)} 

= ggTr„ (S o(t)}. 

b O 


(4.4 


where g is the g factor and g the electron Bohr magneton. The last 
equation illustrates the obvious fact that to calculate spin properties 
we do not need the complete density matrix but just a spin densi ty 
matrix o; this is obtained from the complete density matrix by summing 
over all degrees of freedom except the spin. Tr^ indicates the remaining 
trace over the spin variables. In the following we shall therefore 
only concern ourselves with the problem of finding an approximate 
expression for a(t). A natural starting point to do this is to divide 
the total system in two parts: a spin system and a bath. Neglecting 
intermolecular spin couplings, the spin system will be a single radical. 
The bath (or lattice) is weakly coupled to the spin system; it affects 
the spins and as we shall see later, causes their relaxation, but it is 
negligibly affected by the spins. Corresponding to this separation the 
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total hamiltonian is 


H = H 0 + «b + V 

where is the hamiltonian for the spin system in the presence of 
static magnetic field, is the hamiltonian of the bath, dependent on 
position and orientation of the particles, is the (weak) coupling 
term. In accord with these definitions 

H o "b 2 [V"b] ■ 0; 

H 0 H 1 - I H 0’ H l] * 0; 

K H i 5 [“b'M * 

In the following, it should be remembered that we treat the bath as 
a classical system and the spins as a quantum one, even though a uni¬ 
form notation is adopted. The Liouville equation for the total density 
matrix becomes 


P = -i(H* + + H*)p . 


(4.5 


It is convenient to transform equation (4.5) to a new representation 

x x X 

by writing e.g. p=exp{i(H Q + H^tjp . Equation (4.5) then becomes 


.1 

P 



I 

P . 
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This has the formal solution 

p I (t) = exp{-i fdt' (t 1 ) X }p 1 (0) (4.6 

0 

We now assume, as usual, that at time zero the spin system and the bath 
are not coupled (106,107) so that the total density is 

p(0) = p^o(0) or p (0) = p^ o^(0), (4.7 

where is the equilibrium density of the bath i.e. its Boltzmann 
distribution. Substituting in equation (4.6) and remembering that the 
trace of a product of operators is independent on their representation 
we find 


cr + (t) = exp{-i| dt' H*(t')*}o + (0), (4.8 

1 0 

where we have introduced the definitions o =exp(iHQt)a and 
Hj = exp(il^ t)lL. This representation is often called the interaction 
representation. The upper bar indicates an average over the bath 
variables i.e. Tr^ {(...)p^ } and in practice for our case this will 
mean an ensemble average over the orientations of the radical. We 
are now left with the problemof evaluating the right hand side of 
equation (4.8), To do this we shall follow Freed (108) and use Kubo's 
cumulant method (72). According to cumulant theory 

exp{>-i Jdt' Hj(t') X }= exp K(t) , 

' 0 


(4.9 
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The cumulant function K(t) is given by (108) 
t t 

KCt) - l i, | dtj... f dt n <T (tj) X H*(t 2 ) X H*(t 2 ) X ...«;(t ii ) X > c 
n =l 0 0 

* r t r tn -' + x + x 

= l dt x dt 2 ,,. I dt n <H 1 (t 1 ) ...Hj(t^) > c (4.10 

n=l J 0 0 

where is a cumulant average and T is a time ordering operator, 

which ensures that the time dependent operators on its right are 
arranged in chronological order, Substitution of equation (4.9) in 
equation (4.8) and differentiation gives the equation of motion 

o + Ct) = K(t) o + (t). (4.11 

The terms in equation (4.10) vanish if any of the H^(t) are 
uncorrelated, so that the only nonvanishing contributions come from 
times less than a characteristic correlation time t. It is then rea¬ 
sonable to look for an approximate asymptotic solution K(t-K*>) valid 
for t>>r. If this condition is satisfied the upper limits in the 
time integration can be taken as infinity. Since we can define 
HQ=HQ+H^(t) and H^(t)=Hj(t)-H^(t) we shall take Hj(t)=H*(t)=0 in 
the following. Using this fact and cumulant algebra it is not 
difficult to show that, to second order, 

K(o») « - ] dt' H*(t) X H*(t-t') X (4.13 

0 1 1 

We shall retain only this term in the expansion, an approximation 
which is expected to be valid only if J(t) | 2 r 2 «l where r is a 
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characteristic correlation time. Transforming back from the interaction 
representation we find 


£(») - expCiHp t) X R, 


(4.14 


where 


R - r dt' H 1 (t) X exp(iH X t») H 1 (t-t*) X . (4.15 

0 

The (super) matrix R defined in this way is just the celebrated 
Redfield matrix. It is easy to show that in a basis of eigenvectors 
of Hq, where (a > = E a |a >, the matrix elements of R are explicitly 


R 


aa’bb' 


L abb , a'^ab ) 


+ ^tbb'a^Va' 5 " 6 a'b* 


^ L accb^ u, cb^ 


^ab ^ ^b'cca’ ^b'c^ * 
c 

where 6 , ia a Kronecker delta and 
ab 


(4.16 


L abb r a ? ^ E <a | (t) | bxb' | Hj (t-t') | a '> exp(-ioit ’ )dt', (4.17 

J 0 

The quantity 'a' ^ w ab^ * 8 t ^ ien defined 38 the Fourier-Laplace 

transform at frequency ui a ^= (E a '-E^)/ll of the dynamic hamiltonian 
correlation function. In practice the imaginary part of the transform 
can be neglected (109) and the Redfield matrix is just 


^aa'bb’ ^aba 'b' ^°ab^ ^ab ^ J acbc^ U cb^ 

c 

^ 5 , y j, * . (id. * ), 

ab L b'ca'c b'c 


(4.18 


c 
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To derive equation (4.18) we have used the fact that is an Hermitian 
operator and introduced the usual notation 

J aba'b’ (a>) = l 0 <a|H 1 (t)|bxa , |H 1 (t-t»)|b->^osa,t'dt’. (4.19 

The following relations are easily proved 


J aba'b' (a)) = J b’a'ba (a)) = J aba'b'^ = J bab'a* (u) * 


(4.20 


If we now consider again the equation for the spin density matrix the 
following results are obtained immediately. In the interaction repre¬ 
sentation 


• + 

a aa 


,(t) = l exp[i(a )aa , -o. bb ,)t] R aa , bb , o bb ,(t) (4.21 


bb' 


and in the ordinary representation 


o aa ,(t) = -ia> aa , o aa ,(t) + l R aa , bb , a bb ,(t). 

bb' 


(4.22 


In matrix form we have 


a = (-iH* + R)o. (4.23 

If the relaxation term R is zero the solution is just an oscillation 
about , if R is non zero the oscillation becomes damped, as we might 

expect. The real part of R is the damping term which relaxes the ini¬ 

tially non-equilibrium situation. It is obvious from its physical 
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significance and also verifiable explicitly that R does not act on 
the equilibrium density matrix i.e. Ro e ^=0. We could then write also 

a = ’-iH* o+ R(a-o e q), (4.24 

as is sometimes found in the literature (110). 

Equation (4,23) describes the evolution of the spin density matrix 
■under the action of the bath. In an electron resonance experiment 
We are usually interested in measuring the complex magnetization in 
a plane perpendicular to the static magnetic field B^. A weak radio¬ 
frequency field perpendicular to B Q is applied and the rate of change 
in the magnetization of the sample is monitored. To keep into account 
the interaction with the radiofrequency field B^ we assume that, being 
weak, it produces only a linear response in the density matrix. We 
then have 


9o 

at 


• /8cr s 

_1 V + ( at ^relax. 


+ (— ) r 

V 3t r,f . 


(4.25 


= {-iH* + R -iH 


X 

r ,f 


(t)}o. 


(4.26 


where R is the relaxation matrix previously determined and H is 

r »i i 

the perturbing Hamiltonian: 


H , = uj- (S exp(-iwt) + S exp(iwt)}, (4.27 

r »x * 1+ *■ % 


tu 1 = ggBj/h. 


(4.28 
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By Fourier transformation of equation 4.26 we find 




o(cd) ={ iaj-R + iH Q } {-iai 1 S + a(2u>) -iuj 1 S_o(0)} (4.29 


where 


o(u>) - J 


a(t) exp(-iu>t)dt 


(4.30 


and o(0) is the equilibrium spin density matrix. The non-resonant 

x . 

term u)^S + o(2u)) on the right can usually be neglected and, if we assume 
B^<< B q and gBB^ « RT, the other term becomes (111) 


-io)S X o(0)= -iCS , 


(4.31 


where C is a real constaht. For an unsaturated absorption experiment 
the line shape is (111-113) 


S(u>) lm{Tr (S x o_j) } = Im £ <a | | a-lxa-11 a | a> 


Im{S + .o}, 


(4.32 


where S + is the (2S+1) dimensional vector formed with the elements 
<a|S + |b>. Substitution of the expression already found for a(m) gives 


S (w) « -Re {S_ MU)S_ } , 


(4.33 


where we have defined 


M = -iH* -irn+R 


(4.34 
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M 


aa *bb' 


i (u> . -oi) 6 


a a 


ab 


'b' 


+ R 


aa 


'bb' 


In electro nee the first derivative of S (m) is 

measured. In .. -t formalism this is just 


S'U)=-x. ' }. 


Let us now examine the spectrum limiting 
the R matrix is diagonal. With this resti. ’ is 
and we can write 


M a-l,a,b-l,b = {iU a “ “ )+r a }6 ab* 


It follows that 


S(u.) - Re { l (S_)^ l a /[i(« a -m) ♦rj} , 

a 

= l 8 a L(m,a) a ,r a ), 

a 

where g_ = (S+a) (S-a+1) and L(w,jd .T ) is a lorentzian 

ci cl 3. 


LCtDjOjQjlJ) = — Re{i(o)-o)Q) + } 


In our case 


L(o),(j) ,r ) = -Re{l/[i(u> +ImT ) -iu> + Rer ]} 

a a d a 


(4.34' 

usually 

(4.35 

situation that 
diagonal as well 

(4.36 


(4.37 

line shape i.e. 

(4.38 


(4.39 
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Thus, when R is diagonal, the spectrum is a sum of lorentzians; 

(ReT ) = Re R . is a linewidth and ImT is a (dynamic) 

3 cl L^SLycllySl 3 

frequency shift. A common more complex situation arises when a group 

of transitions is degenerate i.e. they have the same transition 

frequency but these are well separated from the others or, in other 

words, when the M matrix is block diagonal. The block corresponding 

to the degenerate set has to be inverted explicitly. Sometimes the 

inversion of the block can be done analytically (c.f. Chapter 7) but 

in general it is necessary to invert numerically the matrix for every 

frequency at which the spectrum has to be calculated. However, the 

situation, even in the general case of degenerate and non-degenerate 

transitions is not as complicated as it may seem. This is firstly 

because we are only interested in calculating certain elements of M 
. -1 . . 

i.e. the M , , , , . An inspection of the M matrix shows that it 

is possible to permute its rows and columns and rearrange all of these 
elements in a block matrix of dimension 2S*2S, which is quite a dramatic 
reduction in size from the original space of dimensions (2S+1) 2 x(2S+1) 2 . 
For a triplet, as an example, the reduction is from a 9*9 matrix to a 
2x2. The contraction of the total Liouville space to the 2S dimensional 
"transition" space can be achieved because of certain restrictions on 
the subscripts of the relaxation matrix coming, physically, from energy 
conservation(100 ). 

Summarizing*the Liouville space formalism offers a compact and 
general way of calculating the spectrum. In the fast motion limit the 
relaxation operator has been shown to coincide with the Redfield matrix. 
A contraction of the total Liouville space is possible reducing its 
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dimension to 2S i.e. to the number of allowed transitions. It is 
worth pointing out that the formalism does not depend on having 
degenerate or nondegenerate transitions as does the usual Redfield 
treatment. An expression is given for the total spectrum i.e. 

S(u>)“Re(S + M *S_) and it just reduces to the proper limiting cases 
for degenerate and non-degenerate transitions. The method is particu¬ 
larly useful in problems like the one of the angular dependence of 
the electron resonance spectrum treated in this thesis since by 
changing the angle the transitions can pass from non-degenerate to 
degenerate and viceversa. 



- 86 - 


4.3 The Static Spin Hamiltonian 

The calculation of the spectrum outlined in the previous 
Section requires the knowledge of the transition frequencies and of 
the relaxation matrix of the system under consideration. The calcu¬ 
lation of the transition frequencies is treated here. 

Consider the total time dependent spin hamiltonian H(t),that 
is the unperturbed spin hamiltonian together with the spin-bath 
coupling hamiltonian H^. Since in the derivation of the Redfield 
matrix we have supposed that the dynamic hamiltonian has zero average, 
in the following we shall define H = H(t) (static hamiltonian) and 
H^(t) = H(t)-H(t) (dynamic hamiltonian). We shall use the notation 
H instead of H^. The transition frequencies are then obtained from 
H and the elements of the relaxation matrix from the correlation 
function of H^(t). The total spin hamiltonian H(t), since it contains 
the coupling between the spin system and the thermal batl\ will depend 
both on spin and spatial variables. We shall find it convenient to 
write the hamiltonian H(t) as a contraction of irreducible tensors, 
which at a certain instant of time is just, 

H =7 (_)P F ^* p V (L ’- p) . (4.40 

L p p 

In this expression p denotes the nature of the interaction, L its 
rank and F its strength, while T is the corresponding spin tensor. 
Since we shall assume that the molecular reorientation is the only 
relaxation mechanism for the spins, equation (4.40) gives also the 
hamiltonian at a given orientation. The choice of using irreducible 
tensors appears therefore as a very natural one, given their simple 
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transformation properties under rotation. The time dependence can 
be conveniently separated out by writing the tensorsT in the molecular 
coordinate system in terms of the same tensors in a coordinate system 
fixed in the laboratory. It is convenient to choose the z axis of 
this laboratory system,which we call the d frame,along the director. 
The spin hamiltonian is then 


H(t) = Y(-) p F (L * p) D L (t) T (L,q) . (4.41 

L U q,-p p 

Provided the molecular reorientation is fast on our experimental 
time scale, the static spin hamiltonian can be obtained from equation 
(4.41) by taking an ensemble average over the orientations. 


H - Y(-) P f ( l,p) 

L V 


D L T a ’ q) . 

q.-p v 


(4.42 


The symmetry of the mesophase can now be employed to restrict the 
number of non-vanishing order parameters. If the mesophase is uniaxial, 
then 


D 


L 

q»-p 




> 


and 


H = Y(-) P f ( l,p) 
L y 



t ( l * 0) . 

y 


(4.43 


If the director is not parallel to the stationary magnetic field we 
require one more transformation to express the spin operators in terms 
of their components in a frame with the z axis along the magnetic field 
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(the B frame). The result is 


H(t) = £(-) P F 


(L>p) qL 


(t) D L T (t,r) 

q,-p r,q p 


(4.44 


and 


H 


I<-> 


-^P f < l »p) d l 


L T (L,r) 
0,-p r,0 y 


(4.45 


For all the practical cases only the angle y between the z axes of 
the B frame and the d frame needs to be specified, so only the 
rotation (OyO) needs to be considered. 

r »q 

The static spin hamiltonian in equation (4.45) determines of 
course the line positions in the spectrum. It can be rewritten 
formally as the static hamiltonian for a cylindrically symmetric species 
i.e. 


H -y D L „ T tt - r) 

L p r,0 p 


(4.46 


where 

F a ’ 0) - V(-) P F (L,p) dJ; , (4.47 

P L p 0,-p 

is called a partially averaged tensor (38). This useful result shows 
that,irrespective of the symmetry of the probe,the line positions are 
formally that of a pseudomolecule with the symmetry of the mesophase 
(c.f. Chapter 1). The spectrum of a probe dissolved in a liquid crystal 
can therefore be quite different from that in the isotropic phase. 

In fact, taking the isotropic limit of equation (4.46) i.e. letting 
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the order parameter go to zero, we see that 

H = £f (0,0) t (0,0) . (4.48 

L y w 

In other words, in an isotropic solvent the anisotropy in the magnetic 
interactions is completely averaged to zero, in the fast motion 
limit that we have assumed so far. As a trivial consequence the 
spectrum is angular independent (isotropic). By contrast,in an 
ordered mesophase non-scalar interactions can survive even in the fast 
motion limit and the spectrum of an oriented probe shows a dramatic 
change from the isotropic one. The number of lines, their position, 
their shape can change; furthermore the spectrum can exibit varia¬ 
tions when the angle between field and director is changed. As far 
as the line positions are concerned, it is not necessary to go into 
details of these changes: having reduced the static hamiltonian to 
that of a cylindrically symmetric pseudomolecule, these details can 
be found in many textbooks (114,115). We shall just report the 
results for the important case of a spin probe containing one unpaired 
electron which interacts with a single magnetic nucleus of spin I. 

The observed resonant fields are, to first order in |a/ uJq|» 

B r = B q - a(y) hm/{g(y)B} , (4.49 

B 0 = » (4.50 

where is the operating frequency of the spectrometer,and 
neglecting any nuclear quadrupole interaction. The electron 
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resonance spectrum contains (21+1) hyperfine lines. The angular 
dependent coupling constant and g factor are 


g(y) ={g 2 + (g 2 “ g 2 ) cos 2 y}* 

and 


(4.51 


a(y) ={g 2 a 2 + (g 2 A 2 - g 2 A 2 ) cos 2 y}^/g(y)- (4.52 

11 II II 11 


The partially averaged components of the hyperfine tensor A are 
defined as 


A - a + (2/3)* l (-) p A (2,p) D 2 _ , (4.53 

II P 

P 

A - a - (1/6)* l (-) p A (2,p) D 2 , (4.54 

1 u 0,-p 

P 


where a is the scalar coupling constant. 
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4.4 The Dynamic Spin Hamiltonian 

The general form for the dynamic hamiltonian H^(t), in the 
laboratory frame,is obtained by subtracting equation (4.44) and 
equation(4.45), as 


Vt) - I(-) p f ( l,p) d l 


-r, 


.{ D 


q,-p 


(t) - D 


0,-p Oq 


} T 


(L-r) 


(4.55 


= l(~) T F L ’ r (t) T L “ r . 
L V U 


(4.56 


We have defined 


.(L,r) 


[(-) P+r F (L ’ P) { X) L (t) - 6 } D L . (4.57 

" y q,-p U,~p uq —r,q 


The spectral densities to be calculated are 


J aba’b'^ = [ <a | H ! | b><a ' | H i ( t+t ') l b ' >COSajt ' dt ' * (4.58 

' n 


Substituting equation (4.56) we find 


J aba 'b' = l <->***' <a l T(L r) ! b >< a, | Ta, ^ (u) » (4 * 59 


where 


J PP ,(m) 

rr' 


F ( L » r )( t ) F^ ,,r,) (t+t')* coswt'dt'. 


(4.60 


The problem of evaluating the Redfield matrix elements separates 
therefore in a spin part consisting in the calculation of certain 
matrix elements and in a properly dynamic part i.e. the calculation 
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of correlation function transforms The spin part can be dealt 

with quite simply by using angular momentum techniques. In practice 
the calculation of the spin part is the same as for isotropic systems 
provided the eigenstates of H are still the same. Therefore quite 
often the results for the spin part of the problems can be found in 
the standard works on spin relaxation in isotropic systems (109,116). 
We just mention here some well known results for a probe with a single 
electronic spin interacting with a nucleus of spin I. This includes 
the case of the popular nitroxide spin probes whose spectra we shall 
study in the next Chapters. As mentioned before,the calculation of 
the matrix elements of the Redfield matrix ought normally to be 
performed in a basis of eigenvectors of the static spin hamiltonian H. 
It is convenient however to invoke the so called high-field approxima¬ 
tion which corresponds to retaining only that part of the static 
hamiltonian which commutes with and to work in the simple product 
basis |S,mg> |l,nLj.> (see however Chapter 6). When the correlation 
times for the reorientation of the probe is greater than the inverse 
of the resonant frequency the terms which do not commute with S 

Lt 

(i.e. the non-secular terms) can be neglected also in the dynamic 
spin hamiltonian. In this case the linewidth of the transition 
labelled by m^ * m is just the inverse of a Redfield matrix element 
and can be written simply as (116) 


T 2 ’ (m) 


2 

A + Bi+ Cm 


(A.61 


where 
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A = (2o>q /3g 2 ) J 88 + {I(I+l)/4} (4.62 

B = (4to 0 /3g) (4.63 

C = (2/3) - (1/4) (4.64 

To obtain these expressions only contributions from the hyperfine A 
tensor and from the g tensor have been considered. A simplified notation 
for the spectral densities ,e.g., instead of J>> is adopted. 


4.4 (b) The Correlation Function 

In the previous sections we have shown that in the so called 
fast motion limit, i.e. when the bath relax much more quickly than 
the spin system, the linewidths and in general the shape of the spectrum 

lili * 

depends on the spectral densities J^_, . From equations (4.60) and 
(4.57) it is easy to see that the only dynamical quantity to be conside¬ 
red is the correlation function of the Wigner rotation matrices, that 
we write in general,for a stationary random process, as 


&v<‘> 


(dJ; (0) - D L ) (D L ! ,(t) - D L ,’ ,) 

q,p q,p q ,p q f p 


D n n (0) 

q»p 


JL’ 


.(t)* - D“ n 

p q.p 


d l ; *, 

q ,p 


(4.65 


In practise we shall only need here correlation functions with L=L'=2 
so the superscripts L,L' will be left implied. It is possible to 
restrict the number of subscripts in ^qpq'p' by employing the symmetry 
of the mesophase and that of the probe (c.f. Chapter 5). As an example, 
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for a cylindrically symmetric probe and mesophase, the only case 
treated in this thesis. 


t> , , (t) = $ 6,6, 

qpq p qpqp qq PP 


* (t) 

qp 


(4.66 


The calculation of the correlation function requires some knowledge 
of the microscopic dynamic reorientation process. We assume it is 
a stochastic, Markov process (117), so that 

D n n (0) D o T,* (t) A d «n P( V P < fi nl nt) D n < 4 * 67 

q,p q,p J 0 0 q,p 0 J O' q,p 


where P(ft^) is the equilibrium orientational distribution function 
and P(ft Q |ftt) is the probability that, if the molecular orientation 
is ft Q at time zero, it will be ft at time t; it is called the conditio¬ 
nal probability. 

The equilibrium probability is given by 


P(ft) = exp{-U(ft)/kT }/ exp{-U(ft)/kT }, 


(4.68 


where U(ft) is the pseudo potential for the probe molecule. 

According to a molecular field model of binary mixtures of rod-like 
particles (118) the pseudopotential for a solvent and solute molecule is 


U (1) (cos8) = l { (1-x) u^ U) P^ L) + x u^ 12) P< 2) } P L (cos8), (4.69 


L 

(even) 
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U (2) (cos8) = l { (1-x) u t (12) + x u< 22) P^ 2) } P l (cos6), (4.70 

L 

(even) 

where the labels (1) and (2) indicate the solvent and the solute, 
x is the mole fraction of solute, the coefficients u^^^are interaction 
parameters between molecules of species i and j. In electron resonance 
the solute (probe) concentration is vanishingly small. Taking the 
limit for x~K) of equations (4.69), (4.70) one finds (119,120) 

U (1) (cosf3) = l u^ U) P^ l) P L (cosB) 

(even) 

* kT l A ^ 1} P L (cosg) 

L 

U (2) (cosS) = l u< 12 > P^ 1} P l (cos6) 

(even) 

= kT l \ < 12) P L (cosB) 

L 

(even) 

Thus, when the probe concentration is very low, the two simultaneous 

equations (4.69).(4.70) become decoupled. The solute pseudopotential 

( 2 ) . . ( 12 ) 

U depends on the solvent-solute interaction u and on the solvent 

L 

order parameters the solvent pseudopotential is the same as that 

of the pure mesophase. The quantities directly measurable in electron 

- ( 2 ) 

resonance are the solute order parameters P (L=2,4). They can be 

L 

used, if the solvent order parameters are known, to extract some informa¬ 
tion on the solute-solvent interaction (118). By choosing the probe 


(4.71 


(4.72 
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to be as similar as possible to a nematogen, in the sense that 

u (11) some insight can be gained on the mesophase itself (83). 

1 _« 

From a more formal point of view, notice that the solute pseudopotential 
in equation (4.72) has the same angular dependence as the solvent one. 

In particular by truncating the series at the first term (L=2) the 
Maier-Saupe like potential studied in Section 2.3 is recovered. 

In this and the next Chapters the superscripts labelling solute 
and solvent will be omitted unless confusion can arise. 

Let us now turn to the calculation of the conditional probability 
P(fi Q |fit). It is known that the conditional probability for a Markov 
process obeys an integral master equation and that under certain condi¬ 
tions (117) this can be reduced to a differential equation. We shall 
assume this to be the case and consider in brief two rather different 
limiting situations. 

4.4 (c) The Strong Collision Model 

The strong collision model (121,122) represents a process where 
the molecule of interest undergoes sudden changes in its orientation 
at time intervals t. The time taken for the transition from one 

orientation to the other is supposed to be negligible; the orientational 

\ 

probabilities before and after the sudden change (collision) are assumed 
to be given by the equilibrium Boltzmann distribution. Under these 
conditions the rate of change of the conditional probability P^^lftt) is 
simply 


( 12 ) 

"L 


3 “ P(« 0 |nt) = i (p(n) - p(n 0 |nt)}. 


(4.73 
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The general integral of this first order linear differential equation, 
subject to the initial condition P(ftg|fiO) = 6(ft^ -ft) is 


P(ft Q |ftt) *6(ft-ftg) exp(-t/x) + P(ft){ 1-exp(~t/x)} 

= P(ft) + (5(ft-ft 0 ) - P(ft)} exp(~t/x) . (4.74 


This conditional probability obviously reduces to the equilibrium 
probability for very long times, as it has to be (117). Thus, in 
general, the quantity 


D n n (0) D n* n ^ t ^ 

q.p q.p 


L 2 

will reduce to |d | at very long times and the correlation function 

q.p 

<t> (t) defined by equation(4.65) will then decay to zero. It is 

q.p 

easy to calculate the correlation function <f> at any instant of time 

q.p 

t employing equation (4.74). The calculation shows that for a cylindri- 
cally symmetric probe and mesophase equation (4.66) is obeyed, and that 


♦ _ _(t) =( T> 2 T> 2 * ~ P~ 6 6 } exp (-t/x). 

q.p q.pq.p z uq up 


(4.75 


Although this function has been employed to study the electron resonance 
spectrum of vanadyl chelates dissolved in a nematic mesophase (40), the 
use of a single correlation time ignores the anisotropic nature of the so¬ 
lute. This deficiency has been removed by assuming (52) that the corre¬ 
lation time depends on p but not q and so 
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♦_ (t) ={ D^* - 6 6 } exp(-t/x ) 

q,p q,p q,p 2 Oq Op r p 


(4.76 


In linewidth theory (c.f. the previous section) the second subscript 
of p turns out to indicate the component of a certain (2nd rank) 
tensorial interaction so that this generalized strong collision 
model allows each component to decay with its own correlation time t 

\ 

Another form of equation (4.76) can be obtained by coupling the two 
Wigner rotation matrices as (34) 


2 2 * 

D D 

q.p q.p 


,q -p _2 _2 


(-) 


D D~ 

q.p -q,-p 


(-) q " P l C(22L;q-q) C(22L;p-p) . (4.77 


This gives 


*,,,(« ■ ' P l < C(22L;,-q) C(22L;p-p) P^La 6 0, V“ P< ' t/T p > 

L 


(even) 


(4.78 


For linewidth studies it is often useful to have the zero frequency 

Fourier transforms of $ (t), i.e. the spectral densities i (0). 

q.p P J qp 

They are just 


i„p<°> - <-> q ' P l ( C(22L;q-q) C(22L;p-p) V^lA, %> >T p • 

L 

(even) (4.79 

After evaluation of the Clebsch-Gordan coefficients we find the explicit 
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expressions listed in Table 4.1. 


TABLE 4.1 


The spectral densities i 

qp 


for the strong collision model 


Spectral density j^(0) 


f— + — P + -I?- p -P^lr 
*5 7 P 2 35 P 4 P 2 } 0 


±1 


{ I + Ip - li p }T 

l 5 7 z 35 V T 0 


±2 


<T ' f P 2 + 3T V T 0 


±2 


{ 5 ' 7 P 2 + 35 V T 2 


±1 


±2 


{ 5 ” 7 P 2 ' 35 P 4 }t 2 


±2 


±2 


<T * 7 F 2 * 75 V T 2 


Notice that the following symmetry relations are obeyed: 


qp 


q»-p 


J_ 


q,p 


J -q.-p 


(4.80 
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4.4 (d) The Diffusion Model 

The equation of motion for the conditional probability for a 
particle undergoing a rotational Brownian motion subjected to a poten¬ 
tial U(0) is, in coordinate free form (123), 


- P(n 0 |nt) = r P(n Q |nt) 


(4.81 


= ~J{ (D J j^r) P(Q 0 |0t)}- J D J p(o 0 |nt). 


This is subjected to the initial boundary condition 


P(fi 0 |fl 0 ) = 6(n-o 0 ). 


(4.82 


In equation (4.81) T is called the stochastic operator, J is an 
angular momentum operator in a particle fixed frame and in dimensionless 
form and finally D is the rotational diffusion tensor of the particle. 

The single particle potential U is a function of the molecular orientation. 
When U is a constant,equation (4.81) reduces to the well known equation 
of rotational diffusion in a isotropic medium; the eigenvalues and 
eigenfunctions of the stochastic operator are in that case (124) 


E Jp = D a U(J+ 1) + (D |( /D x -l)p 2 }. 


(4.83 


and 



(4.84 


if the diffusion tensor is cylindrically symmetric, as we shall assume 
from now on. The presence of the term containing the potential in 
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equation (4.81) causes in general T to be non hermitian. However, 
it can be shown that, if the detailed balance principle holds, r can 
be made hermitian by the transformation 


r 




(4.85 


where P is the equilibrium orientational distribution function and 
T is the hermitian stochastic operator. The existence of this 
transformation guarantees that T can be diagonalized and that its 
eigenvalues are real. The eigenvalues of r constitute a non-negative 
sequence; the reciprocal of the eigenvalues are the relaxation times 
of the process. The first eigenvalue of r is always zero, corresponding 
to the equilibrium probability which, being unaffected by the stochastic 
operator, has an infinite relaxation time. We shall employ the symme¬ 
trized version of the diffusion equation (98) in Chapter 5; by now 
we shall instead examine briefly the unsymmetrized problem, originally 
solved by Nordio and his group (96). They express r in Euler angles, 
where it has the form 


r p(n 0 |nt) = v* P(n 0 |fit) + (£ ^ {sin e .p( fi() |nt) ^ , ( 4.86 


2 

where 7 ^ is the Laplacian in terms of the Euler angles (agy) descri¬ 
bing the orientation of the molecule with respect to the director. 

The operator r is then given a representation in a Wigner rotation 

2 

matrices basis (i.e., as we have seen, the eigenfunctions of J D J =v^). 
For a cylindrically symmetric probe reorienting in a uniaxial mesophase 
it is clear that the angular momentum projection quantum numbers q 
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and p are still good labels for the eigenvalues of r. In other words 
the total representation of T factorizes in diagonal blocks labelled 
by q and p. In Nordio's notation the (q,p) block is called -D R qp 
and one has 



<L r qp 

dt J 


I <« qP > JJ , c q ? . 

J' 


where the vector 


c qp 


is defined through the expansion 


(4.87 


Pta 0 |nt) - [ c’P(t) D^ pte ). 


(4.88 


The explicit expression for the R 


qp 


matrix is 


(96) 




- -[J(J+1) + (D |( /D a -1) P 2 J6 jj( (4.89 

-J l A L [J(J+1)-J'(J'+l)+L(L+l)]c(LJ'J;Oq) 

L 

* C(LJ’J;0p). 


The parameters A are the expansion coefficients for the orientational 

L 

potential of a molecule, 

U(B) = kT A P (cosg). (4.90 

L L 

L 

(even) 

The axis of symmetry of the diffusion and of the ordering tensors 
are assumed to coincide. It has been shownby Nordio et al (96) that 
the solution to equation (4.87) subject to the initial condition 
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equation (4.82) is 


P« 0 |at) - £(21+1) (X CJP ) JK (X qP )^ e xp(D 1 rJ“> t) D- p (B 0 ) D^ p «) 


9P n 1 * 


(4.91 


where X is the transformation that diagonalizes R: 


X’ 1 R X = r. 


(4.92 


A lower case notation is introduced for the diagonal eigenvalues 
matrix. The zero frequency spectral densities calculated with this 
conditional probability can be written as (96) 


3 qp“ 3 qp^ 0) ^ l D q»p D q,p 6 qo ^po D o,o D o,o }T jqp’ (4 ‘ 93 

j=o 


.2 J 


where x T are time constants defined by 

Jqp 


2J+1 rR qp-i-i 
T Jqp Td7 ( } 2J 


(4.94 


The spectral densities obtained from this model obey certain 
symmetry rules because of the properties of the Clebsch-Gordan 
coefficients, in fact (96) 


J = J = j = j 

q*p -q,-p q,-p J -q,p 


(4.95 
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4.4 (e) Computational Problems 

The calculation of the spectral densities j can be performed 
by inverting the matrix R qp in equation (4.89) to determine the time 
constants Tj a nd then substituting in equation (4.93). The dimension 
of R corresponds to the numbers of terms retained in the series 
expansion (4.88) and is chosen to ensure convergence. A typical value 
is 30. The inversion of R has to be repeated for every (q,p), for 
any D| ( /D A and, of course for any order parameter one is interested 
in. Computation time requirements can then be quite heavy for a 
complete set of calculations. Hemminga (128,129) quotes about 4 hours 
on a CDC CYBER 74.16 computer. Here we shall investigate how to 
reduce dramatically this time while getting a better insight into the 
theory. As a starting point we decompose R in three parts. 


R = a + d + XB , (4.96 

where the matrices a, d, B are defined as 

ajjl = -J(J+1) 6 JJf , (4.97 

d jj - ' - ( §7 - ^ { jj - - <‘- 98 

Bjj, m {J(J+l)-J'(J'+l)+6} C(2J'J;0q) C(2J'j;0p). (4.99 

We use a lower case notation for diagonal matrices. We assume, for 
the sake of simplicity, that the pseudopotential is of the Maier-Saupe 
type U(x) = kTA (x). The following arguments can be generalized to 
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more complex potentials. 

Suppose now we want to change the diffusion tensor. If p=0 

(i.e. if we calculate j _) there is no dependence on D.. /D. and then 

qU 11 i 

no problem. If p+0,notice that d^, is a multiple of the identity 
matrix and that it is the only part of R that contains D^/D^. We can 
write the matrix R as 

R = (a + XB) + d 


(A.100 


= C + d. 


Now,if Q is the matrix that diagonalizes C, 


then 


CQ = Qc , 

Q _1 RQ = c + Q -1 dQ 
■ c + d , 


(A.101 


(A.102 


since d commutes with Q. Therefore the inverse of R is 


so that 


R _1 = Q(c+d) _1 Q _1 


(A.103 


T Jqp 


2J+1 

5 


Q 2K Q KJ 


C K " ( V 1)p2 


(A.10A 


in units of D ^ . It is then not necessary to repeat the calculation 

of R 1 for every D (j /D^ : we just have to diagonalize the matrix C and 
invert its eigenvectors matrix Q once. Let us now return to the spectral 
densities j . Substituting equation (A.10A) in (A.93) and exchanging 

qp 
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the order of the summations on K and J we find 


qp 


-I 


Z K P 


K 


_ / Du 1 .2 
*TC ( D A 1 )P 


where the array Z qP is 



(-l) q+P 

5 



(4.105 


x l (2J+1){ l C(22J\q-q) C(2JJ’,p-p)P Jf 

J=0 J' 

(even) 


- ? 2 ? J 6 0q W 


(4.106 


Equation (4.105) shows very clearly how the ratio D^/D^ affects the 
spectral densities. It allows a very efficient computation of the 
spectral densities and then of the linewidth coefficients in the case 
when the order parameter is kept constant while D (| /D^ is changed, 
whichis what is often required in practice. For example, the order 
parameter can be determined from the line positions assuming a 
Maier-Saupe potential the interaction parameter X can be determined 
e.g. from figure 2.1(notice there {•■-j * is used) and the array Z qP 
is calculated once and for all. The diffusion tensor components ratio 
D 11 /D 1 is then determined by changing it until the linewidth coefficients 
agree with experiment. A slightly more complicated problem arises 
when we want to change the order parameter P^. This is because neither 
matrix B nor matrix a are a multiple of the identity so that the 
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commutativity assumed in equation (A.102) is not guaranteed. We can 
overcome the problem operating some transformations on the matrix R i.e. 


(a+d) * R(a+d) * = 1 + A(a+d) - * B(a+d)~* 


(4.107 


= 1 + AG . 


Now, if S is the transformation that diagonalizes G, i.e. GS«=Sg, then 
after some manipulations 


R _1 = (a+d)“* S(1+Ag) -1 S -1 (a+d)“* , 
and , in units of D^, 

2J+1 ,-i .-i p S 2k S kJ 

T Jqp " 5 22 d JJ 1 (1 + Ag kk ) * 


(4.108 


(4.109 


Therefore, to calculate the time constants t t at various values of 

Jqp 

A (and then of we have just to set up the matrix G, diagonalize 
it and invert its eigenvector matrix once; then perform the summation 
in equation (4.109) for every A of interest. Notice that the formalism 
described here is even more convenient when the symmetrized version 
(98) of relaxation theory is used for then the inverse of the eigenvector 
matrix is just its hermitian conjugate. 
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CHAPTER FIVE 


ANGUI.AR DEPENDENT LINEWIDTHS 


5.1 Introduction 

The magnetic field present in the electron resonance spectrometer 
is usually employed to align the director in the mesophase. It is, 
however, possible to vary the angle between the field and the director 
by proper application of an additional constraint to counterbalance 
the orienting effect of the magnetic field. A number of techniques 
have been developed for this purpose. They include spinning the sample 
(130-132) and the application of an electric field (133-135). As the 
angle between director and magnetic field changes, both the positions 
and the widths of the spectral lines can vary (52-55). From an 
analysis of the angular dependent linewidths in terms of the strong 
collision and diffusion model some static and dynamic informations can 
be obtained. The development of such an analysis should be of value 
also for the interpretation of spectra of spin probes dissolved in 
biological membranes (128, 129, 136-141) in view of the similarity 
between these anisotropic system and a smectic mesophase. 

In this Chapter we report a study of the angular dependence 
of the electron resonance spectrum of a spin probe dissolved in a 
smectic A mesophase. The spectra are analysed to obtain the angular 
variation of the linewidth coefficients and hence the angular linewidth 
coefficients. These angular coefficients may be interpreted by 
employing a theory of spin relaxation based on the strong collision 
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model for molecular reorientation. Such an analysis yields values 
for the correlation times for reorientation about the long and short 
axes of the spin probe as well as the magnitude of the solute order 
parameter P^. The possibility of interpreting the angular linewidth 
coefficients in terms of a theory based on the diffusion model for 
molecular reorientation is also explored. Relatively straighforward 
expressions for these coefficients can only be obtained from the 
diffusion model when the orientational order is high, as it is in a 
smectic phase. The analysis of the angular linewidth coefficients then 
provides values for the components of the diffusion tensor parallel 
and perpendicular to the long axis of the spin probe. The large 
anisotropy in the diffusion tensor may be indicative of the unhindered 
molecular motion about the long axis. 

Both models of molecular reorientation in a highly ordered 
system appear to provide a satisfactory account of the observed angular 


linewidth coefficients. 
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5.2 Experimental 

The mesogen employed in this investigation was 3-N-(4'~ 
ethoxybenzilideneamino)-6-n-butylpyridine. This is reported to melt 
at about 29°C to form a smectic phase which passes to a nematic meso- 
phase at 49°C; the isotropic fluid is subsequently formed at approxi¬ 
mately 54° C (143). The spin probe, (3-spiro-[2'N-oxyl-3',3'-dimethyl- 
oxazolidine ]) 5a-cholestane was synthesised from 5a -cholestan-3-one 
using the standard route (144). The molecular structure of this spin 
probe is shown in figure 5.1 together with the principal coordinate 
system for the magnetic interaction tensors. The transition tempera¬ 
tures for the mesogen were depressed several degrees by the addition 
of the spin probe but this will not affect our analysis of the line- 
width coefficients. A monodomain sample of the smectic mesophase was 
prepared in the following way. The sample was first heated, within 
the cavity of a Varian E-3 spectrometer, until the nematic phase was 
formed and at this point the director, or preferred molecular orienta¬ 
tion, is aligned parallel to the magnetic field of the spectrometer. 

The field was then set to its maximum value of 6.5 kG and the tempera¬ 
ture lowered slowly until the smectic A phase was formed. The surface 
forces, which can inhibit the production of a monodomain sample, were 
reduced by coating the tube with plastic (142). The electron resonance 
spectrum of the sample at 28°C was then measured as a function of its 
orientation relative to that of the magnetic field. The spectra obtained 
for all orientations contained the three hyperfine lines expected for a 
monodomain sample doped with a nitroxide spin probe. In addition the 
angular dependence of the nitrogen hyperfine spacing shows that the 
director remains parallel to the magnetic field on passing from the 
nematic to the smectic phase and so shows this to be smectic A (142). 





Figure 5.1 


The structure of the spin probe and the principal 
coordinate system for the magnetic tensors. 




Ill 


The line shape is essentially gaussian because of inhomogeneous 
broadening caused by unresolved proton hyperfine structure. The 
linewidth was equated with the peak-to-peak separation. The widths 
of the three hyperfine lines were determined from that of the central 
line and the heights of the outer two lines relative to the central 
peak. These linewidths were then used to determine the linewidth 
coefficients in the expression 


T 2 T (m) 


A + Bm + Cm 


2 


(5.1 


where m is the nuclear quantum number. The observed angular 
dependence of these linewidth coefficients is shown in figure 5.2 
and it is this dependence which we shall now examine in greater detail. 


5.3 Angular Dependent Linewidths 

The dominant relaxation process for a spin probe, similar to 
that employed in this study, is its molecular reorientation, with 
respect to the director, coupled to the anisotropy in the magnetic 
interactions. The observation of motional narrowed spectra from this 
system is important for it shows that the linewidths are directly 
related to the dynamic spin hamiltonian. In fact they are determined 
by the Fourier transform of the correlation function of the dynamic 
perturbation. For this problem the dynamic spin hamiltonian is conve¬ 
niently written in irreducible tensor notation as (52): 


H t (t) = 


I (-D 


-I'iP n 2 


D |D (t) -D^ 6_ | T V 

r,q 1 q,~P 0,-p Oq ’ p 


(2,r) 


VJP.q.r 


(5.2 
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The time dependence in H^(t) is contained entirely in the Wigner 
2 

rotation matrix D (t) relating the molecular coordinate system to 

q,“P 

a space fixed system containing the director. This space fixed coordi¬ 
nate system is connected to a second, which contains the magnetic 

, 2 

field, by the rotation matrix D . The orientational order of the 

r,q 

spin probe in the mesophase is given by the ensemble average of the 
. 2 

third rotation matrix, D_ . For douplet state species the only ani- 

0,-p 

sotropic magnetic interactions of any importance are the second rank 

Zeemen and hyperfine terms. The strength of these interactions is 
.<2,p) 


denoted by F 


where the subscript y indicates a particular term. 


(2 r) 

Finally the relevant spin operators are represented by T^ ’ . The 

dynamic aspects of the problem are contained therefore in the correla¬ 
tion function 


. 2 * 


. 2 * 


<j> , ,(t) = D (0) D , , (t)-D“' D" , 6 6 (5.3 

q,-p,-q ,-p’ q,-p q ,-p’ 0,-p 0,-p 0q 0q 


This may be simplified to 


t (t) * D 2 (0) D 2 * (t) - P 2 6 n 6_ 6 

q»-p q»-p q»-p 2 0q Op qq 


pp 


(5.4 

by invoking the cylindrical symmetry of both the mesophase and the 
spin probe as we shall see later on. The complete cross correlation 
function ^ for the interactions y and y' is then 


iJ> (mJ,) (t) = l (-l) r ~ q F (2>p) F ( ?’ p) *C(22L;r-r)C(22L;q-q)<|> (t)P_(cosY), 
r M P m * P L 

L 

(even) (5 . S 

p»q 
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where C(22L;r-r) is a Clebsch-Gordan coefficient (34). The restriction 

of the summation to even values of L results from the D symmetry 

“>h 

of the smectic A phase and the properties of the Clebsch-Gordan 

coefficients require that the maximum value is four. The angular 

dependence of the correlation function is contained in the Legendre 

polynomial P (cosy) where y is the angle between the director and the 
L 

magnetic field. Since the linewidths are proportional to the Fourier 
transform of t|i(t) their angular dependence and hence that of the line- 
width coefficients will be of the form (52) 

A = A q + A 2 (cosy) + A^ P^(cosy). (5.6 

It is important to realise that this relatively simple angular 
dependence only obtains if the nuclear spin is quantised in the 
direction of the magnetic field (137). For the moment assume this 
to be the case. In general however, due to the effect of the pseudo¬ 
secular hyperfine terms in the hamiltonian, the direction of quantiza¬ 
tion of the nuclear spin will change as we change the angle between 
the director and the field. The effect that this has on the linewidths 
is examined in detail in the next Chapter. 

The observed angular linewidth coefficients, listed in Table 5.1, 
were obtained by a least-squares fit of the angular dependence of the 
linewidth coefficients to the theoretical expression in equation (5.6). 
The theoretical angular dependence, calculated with these values, 
is plotted as the solid lines in Figure 5.2. The agreement between 
theory and experiment is seen to be extremely good for both the A and B 



TABLE 5.1 The angular linewidth coefficients (Mrad.s *) 


Component 

L 

a l 

b l 

C L 

0 

52.7 

-24.2 

40.1 

2 

-10.3 

24.6 

-11.2 

4 

8.4 

-0.6 

-12.8 


coefficients but not quite as impressive for C. The discrepancy 
might be attributed to imperfections in the monodomain sample as 
well as to the nuclear spin quantization effect and we shall return 
to this point later. In the following two Sections we shall see how 
these coefficients may be analysed using first the strong collision 
model and then the diffusion model. 
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5.4 The Strong Collision Model 

As we saw in Section 4.4 (c) the 2nd rank Wigner rotation matrices 
correlation function may be written, for the generalized strong colli¬ 
sion model, as 


b (t) = ( dJ n D 2 * - P 2 6 n 6 }exp(-t/x ) (5.7 

q,-p q,-p q,-p 2 Oq Op p 


The angular linewidth coefficients predicted by this relatively crude 
correlation function are (52) 


A 0 = 2u, 0 { 8 <2,0>2 t 0 (1 " P 2 ) + 2 8 (2 ’ 2) T 2 }/15g2+ I(I+1){ A (2,0)2 t 0 (1-P 2 ) 

+ 2A (2,2)2 t 2 }/20, (5.8 

A 2 = 4w 0 ? 2 fg (2,0)2 t 0 (1-P 2 ) " 2 8 (2,2)2 T 2 )/21g 2 + 

1(I+1)P 2 {A (2,0)2 t q (1-P 2 ) - 2A (2,2)2 t 2 }/28, 

(5.9 


\ = H {3 8 (2 ’° )2 t 0 ( V P 2 } + g(2 ’ 2)2 T 2 V /35g2 


-1(1+1) {3A 2,0)2 t 0 (P 4 -P 2 )+A (2,2)2 x 2 P 4 }/35, 


B 0 = Ao> 0 {g (2 ’ 0) A (2 ’ 0) T 0 (l-P 2 )+2g (2 ’ 2) A (2 » 2) T 2 )/15g, 


(5.10 

(5.11 


B 2 = % P 2 {g (2, ° )A(2,0) x 0 (1 -P 2 ) *2g (2 ’ 2) A (2 ’ 2) T 2 }/21g, (5.12 
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B 4 = 8 Wo {3g (2 »° ) A (2 » 0) T n (P / -P?) + g (2 » 2) A (2,2) T,P A )/35g, (5.13 


l 0' 4 2 


2 4 


C Q = {A (2j0)2 t 0 (1-P 2 ) + 2a (2,2)2 t 2 }/12, 


(5.14 


C 2 = 13P 2 ( A (2,0)2 t Q (l-P 2 ) - 2A (2 ’ 2)2 t 2 }/84, 


(5.15 


C 4 = (3A (2 » 0)2 t 0 (F a -P 2 ) + A (2 ’ 2) t 2 P 4 } 17. 


(5.16 


In these equations u>q is the microwave frequency of the spectrometer, 
g is the g tensor, A is the hyperfine tensor and g is the partially 
averaged g factor (c.f. (4.51)). The correlation time is associated 
with fluctuations in the orientation of the long axis of the spin 
probe and x 2 with rotations about this axis. The expressions for the 
angular A coefficients differ slightly from those of Reference (52) 
which are in error. 

McFarland and McConnell (137)were the first to attempt a calcula 

tion of the angular dependence of the linewidth resulting from rotatio 

nal diffusion subject to an anisotropic ordering potential. However. 

it is difficult to compare their results with ours because they adopt 

a more general form for the correlation function. In fact they assume 

an exponential decay for $ (t) as in equation (5.7) but where the 

Q»P 

correlation time is taken to depend on q as well as on p. We shall 
not therefore attempt a comparison here. 

Hemminga and Berendsen 0.40 ) have also implicitly employed the 
strong collision model in their calculation of the angular dependence 
of the linewidths for a spin probe. Their results are not directly 
comparable with equations (5.8) to (5.16) since they assume complete 
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ordering and neglect the pseudo-secular hyperfine terms in the dynamic 
spin hamiltonian. We can achieve some comparisons by setting both P 2 
and equal to one in our equations. In this limit our results for 
the angular B coefficients are in complete accord with those obtained 
by Hemminga and Berendsen (140); there are however differences in the 
A and C coefficients. These result from their neglect of the pseudo¬ 
secular terms and so our results can be employed to gauge the magni¬ 
tude of the errors resulting from this neglect. Thus their expressions 
for the angular A coefficients do not contain the contributions involving 
the hyperfine tensor and their C coefficients are, in our notation. 


C Q = 4A (2 ’ 2) ? 2 /15, 
C 2 = -8A (2,2) ? 2 /21, 

C 4 * 4A (2 ’ 2) % 2 /35. 


(5.17 

(5.18 

(5.19 


We see then that ignoring the pseudo-secular terms creates considerable 

error in C„ but is more successful for C_ and C,. 

0 2 4 

We must now consider the analysis of the observed coefficients 
using the equations obtained from the strong collision model. Of the 
various parameters in equations (5.8) - (5.16) only three, Tq,t 2 anc * 

P^ are unknown. Thus the g and hyperfine tensors have been determi¬ 
ned for (3-spiro- [2'N-oxyl-3',3'-dimethyloxazolidine]) 5a-androstane- 
178 A ol, whose structure is virtually identical to that of the spin 
probe (136) and it is reasonable to employ these tensors. The compo¬ 
nents of the g tensor, in the coordinate system shown in figure 5.1, are 
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g = 2.0089, 
x 

g = 2.0058, 

y 

g z = 2.0021, 

and those for the nitrogen hyperfine tensor are 

A = A = 101.78 Mrad.s * 
x y 

and A = 540.35 Mrad.s 1 . 

z 

When calculating the irreducible components F^ 2,p ^ from these 
results it is important to relise that they are required in the 
coordinate system which contains the supposed molecular symmetry axis. 
Since the y axis makes an angle of about 20° with this symmetry axis 
the required components of the irreducible tensors are (120): 

g (2 ’ 0) = 0.0007, 

(2,2) = -0.0032, 

8 

The remaining parameter is the solute order Y^ which is found to be 
0.89 from the angular dependence of the hyperfine spacing. These 
parameters were then used to extract the two correlation times and 
t 2 as well as the solute order from the three angular C coefficients 
with the results shown in Table 5.2. A second set of values was 


a ( 2 ,°) „ _ 171>5 Mrad.s" 1 , 
A (2,2) = 216.1 Mrad.s" 1 . 
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determined from the B linewidth coefficients and these results are 
also included in Table 5.2. 

TABLE 5.2 The strong collision parameters 


Parameter 

From B 

From C 

Tq/ s*rad 1 

1.4 xl(f 8 

3.9 x10~ 8 

x 2 / s*rad 1 

2.0 xio" 9 

2.4 x10~ 9 

*4 

0.63 

0.73 


There is reasonable agreement between these two sets although it is 
not possible, at this stage, to attribute the differences either to 
a failure of the strong collision model or to experimental error. 

The angular A linewidth coefficients have not been treated in this 
manner because they do contain a contribution resulting from unresolved 
proton hyperfine structure. We have however calculated the A coeffi¬ 
cients using equations (5.8-10)and the two sets of parameters given 
in Table 5.2; these calculated coefficients are compared with the 
observed values in Table 5.3. 

TABLE 5.3 The theoretical and experimental A coefficients (Mrad.s 


Coefficients 

Experimental 

From B 

From C 

A o 

52.7 

32.8 

54.0 

A 2 

-10.3 

-14.0 

-11.6 

A 4 

8.4 

8.3 

6.5 








120 - 


The calculated values of the A_ and A.linewidth coefficients are in 

2 4 

reasonable agreement with their experimental counterparts. In 
contrast there is a large discrepancy between the observed and mean 
predicted values of A^. However this difference, which is about 10 
Mrad.s \ is comparable to the inhomogeneous line broadening caused 
by the unresolved proton hyperfine structure, and so does not indicate 
the presence of other relaxation processes. 

One of the important consequences of such linewidth studies 
is the ability to determine the solute order which is extremely 
difficult to determine in any other way. A knowledge of this additional 
parameter is important because it provides another test of the stati¬ 
stical theories of liquid crystal mixtures and, by implication of 
liquid crystals themselves. The cylindrical symmetry of both the 
spin probe and the smectic A phase dictate that the orientational 
energy of the probe is given by 


U(p) =kT ^ * L P l (cosB), (5.20 

L 

(even) 

where /3 is the angle between the molecular symmetry axis and the 
director. According to theory (c.f. Section 4.4(b)) when the solute 
concentration is vanishingly small the expansion coefficients are 


X L ■ f L U ' kT - 


(5.21 


where P^^ is the solvent order and u^^ f s a measure of the solute- 
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solvent interaction. The solute order P is obtained by taking the 

Li 

appropriate Boltzmann average: 

P L (cos8)exp{-U(B)/kT}sinBd8 

_ . (5.22 

f exp{-U(8)/kT}sin8d8 

0 

Consequently if both P and P^ are known it should be possible to 
determine two of the expansion coefficients X and so test the 

L 

expected (118-120) convergence of equation (5.20). Unfortunately 
this determination is not possible when the solute order is high as 
we shall now discover. 

If the orientational order is high, as it often is in a smectic 
phase, then the Boltzmann factor will be a rapidly decreasing function 
of 8- Accordingly, in this high order limit, only small values of 8 
make an important contribution to the integrand and we therefore 
obtain a simple expression for P by expanding the Legendre polynomials 

L. 

2 

retaining terms up to 8 . The potential energy becomes 



U(8) 

kT 


c 0 * ee 2 


(5.23 


where 


and 


L/2 

£ = ~ l l (L/2 - r)X L a L>r /kT, (5.24 

L r=0 
(even) 


a 


L,r 


(-l) r (2L - 2r)! 


2 L rI (L-r)'(L-2r)I 


(5.25 
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The ensemble averages P are independent of £ and so we shall not be 

L U 

concerned with its magnitude. Since the integrand is insensitive to 
high values of 6 the upper limit of the integrals may be extended to 
infinity and this then gives 


P 2 = 1 - 3/2$ + 0(£ 2 ), (5.26 

and 

P 4 = 1 - 5/5 + 0(C _2 ). (5.27 

In this high order limit the magnitudes of both P 2 and P^ are 
determined by the single parameter £ and so there is apparently little 
to be gained from knowing P^ as well as P^. Notice as an example 
that equations (5.26), (5.27) are in complete accord with the asymptotic 
formulae (2.59) obtained for a simple Maier-Saupe potential such as 

U(B) = X 2 kT P 2 (cos6), (5.28 

with X=-2E,/3. In Table 2.1 we compared, not unfavourably, values 
of P 2 and P^ obtained from the asymptotic formulae with values obtained 
from numerical integration. For example when £ is 14 the approximate 
P 2 is less than 1% too high while P^ is about 4% too small. The 
observed solute order of 0.89 corresponds to £ equal to 13.0 and so 
the limiting relations for P 2 and P^ should be relatively good. 

Indeed the value of P^ predicted by equation (5.27) is 0.62 which 
agrees well with the experimental values of 0.63 and 0.73. We may 
therefore conclude that U(8) is indeed a rapidly decreasing function 
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of 6 as we might have anticipated. 

The two correlation times, given in Table 5.2 show that the 

neglect of non-secular terms in deriving equations (5.8)-(5.16) is 

2 2 

a valid approximation for this system for t^> 1 (52), although the 
extension to include these terms is quite straight forward (14). 

However it is difficult to comment further on these times since the 
strong collision model does not indicate how they should depend on the 
molecular structure or environment. In contrast the diffusion model 
does predict the dependence of the correlation functions on both the 
diffusion tensor and the orientational order and we shall now analyse 
our results in terms of these predictions. 

5.5 The Diffusion Model 

The solution of the rotational diffusion, in presence of an 
ordering potential, is usually expressed (96) in a basis of eigen¬ 
functions of the diffusion operator for an isotropic potential. 

As mentioned in Section 4.4 (d), these eigenfunctions are just 
Wigner rotation matrices. This method of solution is clearly perfectly 
general but in practise the number of terms to be retained in the rele¬ 
vant series expansions (4.91), (4.93) increases as the system becomes 
more ordered and the potential more anisotropic. However in this case, 
Polnaszek et al. (98) have obtained a simple limiting solution. They 
start by a symmetrized version of the equation of motion (4.81) for 
the conditional probability, i.e. 


I 7 P(n 0 |nO = f P(n 0 |nt), 


(5.29 
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where P(n Q |nt) = P * P (« Q |nt) , T = P W and P = P(n) is the 

A 

equilibrium angular distribution. The zeroth eigenfunction of T : Gq(^) 
has a zero eigenvalue and is just 


G 0 (n) « P*. 


(5.30 


A formal solution to equation (5.29) is 

p(n 0 |nt) = e _ ^ t p _i (Q 0 ) 6(n- n 0 ), (5.31 

when the orientation at zero time is Since the eigenfunctions 

A 

of r form a complete set they may be used as a representation of the 
6 function and so 


where 


P(n Q |nt) = P 5 (fl 0 ) l e V| X;fi><X;n 0 | , (5.32 

x 


G (fi) H IX;Q> . 
x 


(5.33 


The formalism used up to here can be applied to a quite general 
Markov process. As an example, it can include the strong collision 
model involving a single correlation time t. In this case the eigenva¬ 
lues of the stochastic operator are 


= ( I - 


4 0x> ' T 


(5.34 



125 - 


As another aside we can prove, before passing to describe the 
asymptotic solution, the simplification of the correlation function 
^qpq'p' in e q uat * on (5.3) to the $ in equation (5.4). From the 
formal expression for the conditional probability we find 


<f> , i (t) = 

qpq p 


I 

x*0 


-E 
e x 


C <X;n ()l D q,p ( Vl °»V <0; l D qt jp ,(fi)|X;fi> . 

(5.35 


where the matrix notation, introduced by Polnaszek et al. (98) 
is defined by 


<X;fl|f (fl)|x';n> 


= G (Q) f(fl) G (fi)dfl 

X X 


(5.36 


Since the rotation matrices form a complete set of functions which span 
the space of the Euler angles we can expand |X;S1> as 


where 


x;n> = l 

L;m,n 


XL 

a 

mn 


m, n 


(5.37 


XL 

i 

ran 


= <D (ft) X;ft>. 
m,n 1 


(5.38 


Of course the eigenfunctions of T must transform according to an 
irreducible representation of the system which in this case is 
comprised of the liquid crystal and the spin probe. When both of 
these possess cylindrical symmetry the eigenfunctions are of the form 


lm a 


G X (B) e in 


X;ft> = e 


(5.39 
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Consequently the expression for the eigenfunction in equation (5.37) 
reduces to 


|X;fl> = J a XL D L (n), 
1 L mn m,n 


(5.40 


and the zeroth eigenfunction is simply 


' l D o,o (!!) 


(5.41 


The correlation function is then 


<i> , , (t) = y e E x t y 

qpq p L L 


OL* XL' XL"* OL'" 
3 00 3 mn 3 mn a 00 


x*0 


L,L',L",L"' 


<D 0,0 < Vl D q,p te 0 ) l D ^n ( V > 


('.42 

The matrix elements are just integrals of three rotation matrices and 
are given (34) 


<D 0,0 ( V l D q,p < Vl<’n to o) > * 


8tt 2 


« n C(L'2L;mq)C(L'2L;np) 


2L+1 q+m,0 p+n,0 


(5.43 


and 


T 11 T 1,1 

<D »i," (n)|D q',p' fe)|D 0,0 (1! > > - 


TiTir' 1 " c <'-"'2I-"iO-,')ca'"2L";0-p') 


(5.44 
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Combining the 6 functions from these two matrix elements gives 
the desidered result. 

When the spin probe is not cylindrically symmetric the eigen¬ 
functions of T take the form 

| X;fl> = e lm ' a G x (B,y). (5.45 

Consequently the correlation function can only be reduced to 


<t> , , (t) = <J> , , (t) 6 , 

qpq p qpq p qq 


(5.46 


for such a spin probe. 

Let us return now to the asymptotic solution of equation (5.29) 
In view of equation (5.39) we may restrict our attention to the 
operator 


86 


r . - D { _ii + cot g_9. - P +q : 2qp cosg _ ( Jl. } 2 

ft 1 1 3g2 P 3 ft sin^ft V F 


D 


* Er slnB fe <si " 6 ^ » > (5 - 47 


which may be written as 


V" v 


a 2 

W 


+ 11- 
6 36 


_ p-ql 2 

e 2 


®|| 2 3 

-pq - % - l)p^2? (2+6 1^)} 


(5.48 


in the limit of high order. This result is obtained using the complete 
expansion of the ordering potential; none the less it is identical in 



form to the stochastic operator derived by Polnaszek et al. (98) who 
restricted the potential to just the first term. The eigenfunctions 


of this limiting form of T are found (98) to be 

p 

|M,p,q;g> = /(2£ )N^ exp(-£B 2 )(£*6 ) V l”(£B 2 ), 


where v is |p—q|; the normalization factor is 


(5.49 


Njj =* {2£ M! / (M +v )!} (5.50 

and the function L^(£ g 2 ) is a generalized Laguerre polynomial 
(80). The eigenvalues of the stochastic operator in the high-order 
limit are (98) 


E = - {2C D, (2M+v) + pqD +p 2 (D, -D )}. (5.51 

“ > V 1 I II i 


The normalized distribution function obtained from equation (5.49) for 
the zeroth eigenfunction is 


P 0 (8) = 2£ exp(-CB 2 ), 


(5.52 


which is, of course, identical to that resulting from the limiting 
form of the ordering potential given in equation (5.23). The matrix 
elements of the rotation matrices required to determine the correlation 
function 4>^_p(t) have been calculated by Polnaszek et al.(98); 
they find 
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<M f p,q;fi|D^ m |M’,p ’,q ’; 


^M,M' ^pjp'+m* ^q,q'+m ^v,v ,+ ^p,p'+m' ^q,q'+in ^v, v '+1 ^M,M' *^ M+V ) 

/(M+1 ^ +6 v,v'-1< 6 M,M’ /(M+ v+ 1) “ 6 M,M i + 1 ^ K'.m+lV”’” 0 


m ,m-I — 


(5.53 


where 


f (M,m) ={ M(M+1) - m(m±l)} . 


(5.54 


This expression for the matrix element differs from that in (5.23) 
where, presumably as the result of a printer's error, a factor 1/2? 
occurs instead of 1/2 

We have evaluated the correlation functions <P (t), required 

q ,~P 

in the calculation of the angular linewidth coefficients, from equations 
(5.35^(5.49) and (5.51). The zero-frequency Fourier transforms j _p(0) 
of these functions are given in Table 5.4; the terms with p=±l are not 
shown because the principal axes of the magnetic tensors are taken 
to be parallel to those of the diffusion tensor for this spin probe. 

We can now calculate the angular linewidth coefficients in terms of the 
three non-zero spectral densities j 10 » ji2 and J22J the results are 

A o = 4 “o 2 <g (2,0)2 iio ♦ g (2,2)2 (ii2 + j 22 )>/15i 2 

+ 1(1 + l){A (2 ’ O)2 j 10 + a ( 2 * 2)2 (j 12 + j 22 )}/10 , (5.55 

^2 * 4u> o 2 {g^ 2,0 ^ j 10 + g^ 2 ’ 2 ^ (ji 2 - 2j 22 ) }/21g 2 

+ 1(1 + l){A (2 ’° )2 j 10 + a ( 2 » 2)2 (j 12 - 2j 22 )}/28 , (5.56 



TABLE 5.4 The spectral densities j (0) for the diffusion model 

q.p 

in the limit of high order. 


q 

p 

(0) 

q»p 

0 

0 

0 

±i 

0 

3/4? 2 D A 

±2 

0 

0 

0 

±2 

0 

±i 

±2 

1/2£{2D |( + U - 1)D X } 

±2 

±2 

1/4 D„ 


> 

II 

- 4m o 2 {4g (2 ’ O)2 j 10 + 

g <2,2)2 (Aj 12 - j 22 )}/35g 2 



+ 1(1 + l){4A (2,0)2 j 10 

+ A (2 ’ 2)2 (4j 12 - j 22 )}/35 , 

(5.57 

B o = 

8u> o (g (2 * 0) A (2,0) j 10 + 

g (2,2) A (2,2) (ji2 + j22) } /15 g , 

(5.58 

B 2 = 

n , (2,0) (2,0). 

8(D o (g A y J 10 + 

g(2.2) A ( 2,2) (jl2 - 2j 22 )>/2lg , 

(5.59 

B 4 = 

- 8a )o {4g (2 »° ) A (2 ’ O) j 10 

+ g (2,2) A (2,2) (4j 12 - j 22 )}/35g , 

(5.60 

c o " 

u (2 ' 0,! i„ 4 a (2 - 2)2 

Cj 1 2 + j 22 )}/6 , 

(5.61 

C 2 " 

13{A (2,0)2 j 1 o + A (2 ’ 2)2 (j 12 - 2j 22 )>/84 , 

(5.62 


= - (4A (2,0)/ j 10 + a ( 2,2)2 (4j 12 - j 22 )}/7 • 


(5.63 


Although the angular linewidth coefficients are a function of three 


spectral densities there are only two unknowns namely the principal 
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components of the diffusion tensor. Consequently the problem is over 
determined when all three components of an angular linewidth coefficient 
are known. We have therefore sought to obtain the best values for D () 
and D I in the following way. The three spectral densities cannot be 
evalulated from either the B or the C coefficients because the determinant 
of the appropriate coefficients vanishes. We therefore define new variables 
which for the C linewidth coefficients are 

(2 2) 2 

a = 9 * 322 , (5.64 

, , (2,0) 2 . A (2,2) 2 . /c . 

and b = A J10 + ^ Jl2» (5.65 

similar variables are employed in the analysis of the B coefficients. The 
best values for a and b are now obtained from the linewidth coefficients 

t 

using a linear least squares fit. As a check we use these least-squares 
fitted parameters to recalculate the coefficients; the results of such 
calculations are given in Table 5.5. The agreement between theory and 

TABLE 5.5 The observed and calculated values for the angular B and C 
linewidth coefficients 


Component 

L 





Observed 

Calculated 

Observed 

Calculated 

0 

-24.2 

-24.1 

40.1 

25.3 

2 

24.6 

24.6 

-11.2 

-21.0 

4 

-0.6 

-0.6 

-12.8 

-18.2 


experiment is extremely good for the B coefficients but clearly less 
impressive for the C linewidth coefficients; we shall return to this 
discrepancy shortly. Given the value for the variable a the spectral 
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density 322 can be calculated and from this we can evaluate the parallel 
component of the diffusion tensor. Armed with this component it is 
possible to extract the other component from a knowledge of b. The 
values obtained for D ( and from both the B and C linewidth coefficients 
are given in Table 5.6. The agreement between the values extracted from 
the two sets of linewidth coefficients is encouraging; we shall return to 
the large anisotropy in the diffusion tensor shortly. 

TABLE 5.6 The diffusion tensor 


Diffusion tensor 

From B 

From C 

(Mrad s *) 

137 

122 

D (Mrad s 

2.6 

2.6 

D ,i 

53 

46 


As a check on the validity of the diffusion model, in the limit of 
high order, the components of the diffusion tensor, given in Table 5.6, 
have been used to calculate two sets of the angular A linewidth coefficients. 
The results of these calculations are listed in Table 5.7 together with 
the experimental values of A^. The agreement with experiment is seen to 

TABLE 5.7 The experimental and theoretical A linewidth coefficients 


(Mrad s 

-b 



Coefficient 

Experimental 

From B 

From C 

A o 

52.7 

33.0 

35.6 

A 2 

-10.3 

-13.8 

-16.0 

\ 

8.4 

8.8 

8.2 
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be good especially as the inhomogeneous broadening caused by unresolved 
proton hyperfine structure might be expected to make a contribution of 
about 15 Mrad.s * to A Q . On balance this limiting form of the diffusion 
model appears to provide the most satisfactory quantitative interpretation 
of all three experimental linewidth coefficients. However, it may be 
significant that the asymptotic spectral densities calculated with this 
Polnaszek et al (98) method appear not to satisfy all the symmetry 
requirements (96,129). This important problem is now under investigation. 

It would be unwise to attempt to differentiate between the two models on 
the basis of our experimental results. However one clear advantage of 
the diffusion model is that the dynamic parameters occurring in the theory 
may then be the subject of further rationalisation. We shall now attempt 
such an analysis. 

A combination of experimental (145) and theoretical (143) investigations 
has shown that the principal components of the diffusion tensor are often 
given by (147): 


D. = kT/87rr)o.ic. • • (5.66 

l l li 

In this expression o^ is a geometric factor which takes the values (124) 
a = 4r 2 (r 2 - r 2 )/3(2r - r 2 S) , (5.67 

X Z X z x 

and a = 4(1^ - r 1+ )/3{(2r 2 - r 2 )S - 2r } , (5.68 

Z X Z X z 

for a symmetric top molecule with semi-axes r > r = r and where 

z x y 

S = 2(r 2 - r 2 ) * log e {[r g + (r 2 - r 2 )*]/r x > . (5.69 


The parameter n has the dimensions of viscosity and is frequently 
identified with the coefficient of shear viscosity. When k is the unit 
tensor then equation (5.66) is the expression obtained by Perrin (124) for 
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the diffusion tensor of a symmetric top in the hydrodynamic limit. 
Experimentally k is observed to deviate from unity; it is however found 
to be independent of both temperature and pressure although markedly 
solvent dependent (147). These results have been rationalised by a 
theory (146) which relates k to the intermolecular torque T and the 
intermolecular force F, in fact according to this theory 

k.. = (3/4r 2 )(T?/f2) , (5.70 

11 o 1 


where is the effective hydrodynamic radius of the spin probe. The 
tensor k may be referred to as the anisotropic interaction tensor since 
its magnitude is determined by the anisotropy in the intermolecular 
potential which governs the torques whereas the force depends on the 
total intermolecular potential. The ratio of the components of the 
diffusion tensor is therefore 


d m /D i ' °i K i /o n K n 


(5.71 


which should be independent of temperature and pressure but not solvent. 

The application of this result to liquid crystals implies that the 

viscosity n occurring in equation (5.66) for the diffusion tensor is a 

scalar. At first sight such an approximation might be thought to be 

rather poor for the coefficient of shear viscosity for a nematic mesophase 

is found to be highly anisotropic (6). However this need should not be 

so for it should be recalled that the shear viscosity is a translational 

property whereas the diffusion tensor relates to the rotational motion. 

In addition any attempt to allow for an anisotropy in n would introduce a 

largely arbitrary parameter into the problem. The ratio a /a |( can be 

calculated provided the molecular dimensions are known and we have 

estimated them with the aid of molecular models; r is about 11.5 X and 

z 
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r x is approximately 3.5 X. These dimensions yield a value for /o { f 
of about 3 and so the ratio k of the components of the anisotropic 

interaction tensor is approximately 17. This value is large and will 
naturally depend on the choice of the molecular dimensions, however these 
and hence the ratio ic /k (| are unlikely to be in serious error. 
Theoretically the maximum value of ic is unity and this occurs when the 
intermolecular potential is dominated by the anisotropic component (146) . 
At the other extreme k is zero which corresponds to a purely isotropic 
intermolecular potential. Accordingly the value of 17 implies that k 
is no larger than 0.06; in other words rotation about the long molecular 
axis is essentially unhindered. This result is in accord with the 
cylindrical symmetry of the ordering matrix observed for this spin probe 
dissolved in several nematic mesophases (120). 

Also, a recent independent work by Hemminga (128) on the same spin 
probe in a highly ordered (P^ - 0.85) lyotropic system seems to confirm 
that /D^ is between 40 and 50, and then that k is very small. 

It is possible to check the ratio /D^ by measuring the relative 
linewidth coefficients when the spin probe is dissolved in the isotropic 
melt of the mesogen. We have made these measurements and find that C/B 
is approximately -2.1 and independent of temperature. This ratio is 
related to D j( /D^ by 


C 

B 


15g{A 


( 2 , 0 ) : 


48 % {g (2 ’ 0) A (2 ’ 0) 


(2D n /D A + l)/3 + 2A (2,2)2 } 
(2D || /D 1 + l)/3 + 2g (2,2) A (2 ’ 2) } 


(5.72 


provided the molecular motion is sufficiently slow to make the non-secular 
contributions to the linewidth negligible (126). By using the same 
components of the anisotropic g and nitrogen hyperfine tensor as those 
encountered previously the ratio D )( /D x is calculated to be 44. This 
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result is in excellent agreement with the value determined for the spin 
probe dissolved in the smectic phase especially as D (( /D^ is rather 
sensitive to the magnitude of C/B. The close similarity in the ratios for 
the two phases is unexpected in view of the absence of any long range 
orientational order in the isotropic phase which might be expected to 
reduce k in this phase. Our results for D /D^ presumably imply that 
k is determined by the local molecular environment which is essentially 
unchanged on passing from the smectic to the isotropic phase. Accordingly 
if this interpretation is correct then D (| /D^ should be reduced when the 
solvent is composed of molecules more symmetric than those of the mesogen. 

We have therefore measured the ratio of C/B for the spin probe dissolved 
in o-terphenyl and find it to be -1.7. This corresponds to the smaller 
value for D /B l of 19 and such a reduction is in accord with our 
expectations. 

It is appropriate, in concluding, to comment on those areas in which 
the theory has been less successful. For example the predicted angular 
dependence of the B linewidth coefficient is in complete accord wich 
experiment whereas the agreement for the C coefficient was not so impressive. 
Similarly the diffusion model was able to provide an almost quantitative 
account of the angular B coefficients but only a qualitative description 
of the C linewidth coefficients. These discrepancies can be ascribed to 
two assumptions in the analysis; the first is implicit and the second 
explicit. We have assumed implicitly that the director is uniformly 
aligned throughout the sample of the smectic mesophase. However minor 
deviations from complete alignment could result when the smectic phase is 
formed from the aligned nematic mesophase. Indeed the slight asymmetry 
in the spectral line shape hints at such deviations in uniformity (52,148). 

To first order the imperfect alignment produces inhomogeneous broadening 
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of the spectral lines. For a nitroxide spin probe the contribution of 
such broadening to the B coefficient depends on the anisotropy in the 
partially averaged g tensor g whereas the contribution to C is determined 
by the anisotropy in the partially averaged hyperfine tensor A (52). 
Experimentally g is found to be virtually isotropic for this system 
whereas A is highly anisotropic. Consequently the inhomogeneous line 
broadening will alter the angular C linewidth coefficients but not the B 
coefficients. In addition the angular dependence of C resulting from 
inhomogeneous line broadening is predicted to be different from that caused 
by spin relaxation (52). The behaviour of the B linewidth coefficient is 
therefore expected to be in better accord with theory than that of the C 
coefficient. 

The explicit assumption concerns the quantization of the nuclear 
spin in the direction of the magnetic field. This approximation is exact 
when the magnetic field is either parallel or perpendicular to the director 
but is least reliable when the field is midway between these two extremes 
as we shall discuss in the next Chapter. In addition the assumption is 
less accurate when there is a large anisotropy in the partially averaged 
hyperfine tensor. Again the consequences produced by the inaccuracy of 
the assumption will affect the C rather than the B coefficients and will 
have a different angular dependence to that given e.g. in equation (5.6). 
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CHAPTER SIX 


A MORE COMPLETE THEORY OF ANGULAR DEPENDENT LINEWIDTHS 


6 .1 Introduction 

The theory of angular dependent linewidths employed in the last 
Chapter was based on an explicit assumption, namely that when the 
angle between director and magnetic field is varied, the nuclear spin 
is still quantized along the direction of the magnetic field. This 
is equivalent to neglecting the effect on the linewidths of the 
pseudosecular hyperfine terras in the static hamiltonian. The existence 
of the problem has been pointed out in the literature (137) but the 
effect has simply been disregarded as we did in the previous Chapter. 
Part of the reason for doing so is probably the absence of a simple 
theory for the effect. In this Chapter we shall develop such a 
theory, showing that the linewidths can still be expressed as 

T2^(m) = A + Bm + Cm^. (6.1 

However, the angular dependent coefficients A,B,C can no longer be 
written in the simple form 

X(y) = X q + X 2 P 2 (cosy) + X^ P^cosy) (6.2 


where X = A,B,C. 
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6.2 Theory 

We shall start by reconsidering briefly the static hamiltonian 
for a spin probe reorienting rapidly in a liquid crystal. Fast reo¬ 
rientation is intended here in the sense that Redfield (100) type 
theories (c.f. Chapter 4) are applicable. In a molecule fixed co¬ 
ordinate system the spin hamiltonian is, in irreducible tensor notation, 


h = y <-> p f ( l - p) i ( l ‘ p) , 

L P p 

u;l,p 


(6.3 


where, as in the previous Chapters, the and are the relevant 
magnetic and spin tensors. A transformation of to the director 
frame confines the time dependence to a Wigner rotation matrix. 

A further transformation to a coordinate system with the z axis 
parallel to the static magnetic field gives the spin hamiltonian 
at time t as: 


h( t) = £(-) p f ( l,p) d l (t) dJ; (fl) T (L,r) . 
P q,-p r,q p 


(6.4 


An ensemble average gives, for a uniaxial mesophase, the static 
hamiltonian as 


H = £(-) p f ( l,p) d l 6 d l (n) T (L,r ^ 

L jj q f —p qO r,q y 


(6.5 


=I ~(L.O) d L (L,r) 

L p r,0 p 

where the partially averaged tensors ^^*0) are, as usual, 
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.(L,0) . 


= l (-> 


p „(L,p) 


D 


0,-p 


(6.7 


The eigenvalues of the spin hamiltonian (6.6) determine the line 

positions and its eigenfunctions constitute the basis to be used 

in the calculation of the Redfield matrix elements. For a nitroxide 

spin probe the so called high-field approximation, which amounts to 

neglecting that part of the static spin hamiltonian which does not 

commute with S^, can usually be invoked (109, 116), We shall then 

neglect terms such as S + I + , S+I^, BS + ,i.e. the non-secular terms, 

when considering the irreducible spin operators. In Table 6.1 we 

give the explicit expressions for these spin operators and their 

corresponding spatial tensors. An inspection of the spin hamiltonian 

matrix in the product basis |Smg>|lm I > shows that it is justified 

to neglect the non-secular terms as long as the partially averaged 

hyperfine terms are much smaller that the transition frequency. 

We shall call the static hamiltonian, after the elimination of the 

non-secular terms, the truncated hamiltonian. The time-independent 

transformation D L (n) connecting director and field involves in 
r >. 

principle two Euler angles. However, in view of the cylindrical symmetry 

about the director and field, it is sufficient to consider only 

angle between the z axis of the two frames, i.e. D L (OyO) where y 

r ,u 

is the angle between the director and the static magnetic field. 

From equation (6.5) it is easy to see that if y=0 or y=tt/ 2 the 
truncated hamiltonian is already diagonal, since 


D r,0 (000) = S rO 

D r,0<° 1 0) ’ - \ 5 r0 U I S r2 


( 6.8 



TABLE 6.1 The irreducible components of the Zeeman and hyperfine 


tensors g and A together with that of their corresponding 
spin operators. It has been assumed that the magnetic 
field B = B(001). 


g <0 - 0) 

= /3 g 

T^ 0,0) = - B S // 3 
g z 

g ( 2 .°) 

= /3/2 g' zz 

T^ 2 * 0 ) = ^2/3 B S 
g z 

g< 2,±1 ) 

- 1 KJ 

T (2*±1) = -( B s )/2 

g + ± 

g(2>i 2) 

- K* ~ & 2i 

T ( 2 * ±2 ) = Q 

g 

A (°,0) 

= /3 a 

T <°-°) = S • I//3 

a (2 ’°> 

= /3/2 A' 
zz 

T^ 2 ’ 0) »/2/3[S z I z -(S + I_ + S_I + )/4j 

a (2 ,±1) 

= ;(A' ± iA ) 

+ xz yz 

T «,±l), - (I± s 2 . Sj)/2 

a <2 ,±2) 

= (A 1 - A' ± 2iA )/2 

xx yy xy 

t a 2,±2) = : ± s ± 12 


(a) 

The prime indicates the traceless tensor components. 





141 - 


Therefore, when the director is either parallel or perpendicular to 
the magnetic field B, the set |S=l/2,nig> |l*l, m^ is a good basis 
set and can be employed for the calculation of the Redfield matrix. 
However, when 0 <y<tt/ 2, the situation is far more complicated. In 
fact an inspection of the truncated hamiltonian shows that it is 
not diagonal i.e. it does not commute with S^I^ any more. Its matrix 
consists of two 3x3 blocks corresponding to m =±1/2. The off-diagonal 

D 

elements in the two blocks are the pseudosecular hyperfine terms. 

As far as the static hamiltonian is concerned the problem is formally 
identical to that occurring in single crystal studies of the hyperfine 
structure of a radical spectrum (114, 150). There one has essentially 
two possibilities. One is to express the eigenfunctions of the truncated 
H(y) in terms of the eigenfunctions of H(0) as 

IIn>'= I IIm> D 1 
1 L ' m,n 

m 

using the fact that the eigenfunctions of an angular momentum (spin) 

I transform according to the irreducible representation (151). 

The other, equivalent, possibility is to notice that terms such as 

S z I + can be eliminated by changing to a new set of axes for the nuclear 

spin I (114). Since the partially averaged hyperfine tensor has 

cylindrical symmetry, Q, the direction of quantization of I lies in the 

plane defined by the director d and the field B (114). In figure 6.1 

the various relevant directions are shown. The angle x between Q and 

B is so defined that the coefficients in front of terms like S I, 

z ± 

in the static hamiltonian vanish after the transformation. The explicit 
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Figure 6.1 The axes systems involved in the treatment of the 
nuclear spin quantization effect. In the figure, 
d represents the director, Q the axis of quantization 
of the nuclear spin, B the direction of the static magnetic 


field. 
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expression for x is given by (114,137) 


- 2 - - 2 _ 
cosy = (g^ sin y + g^A^cos y} /g a 


siny = {g^ 


A - g^ A^lsiny cosy /ga 


t 


where 


_2 _2 

g a 




~2 ~ 2 . 

K A ) 
B 1 1 


2 '2 
cos y + g A 



(6.9 


( 6.10 


This results can be further simplified by assuming that the partially 
averaged g tensor id essentially isotropic, an approximation expected 
to be very good for all the practical cases we have to deal with (52) 
This assumption has also the important consequence that as we change 
the angle y the electron spin is constantly quantized along the field 
The nuclear spin quantization axis, instead, changes as the angle y 
changes and it deviates most from the field direction when 


2 tl a 72 . , ~2 ~2 . 

cos y = (A () k L - A ± ) / (A„ - A x ) . 


( 6.11 


At this point 


cosx = cos X max = 2A i / { (A,, + A ± )(A„ A a ) ^ }. (6.12 

The magnitude of the effect clearly depends on the anisotropy in 
the partially averaged hyperfine tensor and is therefore related to 
the order parameter. For example, for a spin probe like the one 
used in the previous Chapter (see figure 5.1), with the symmetry 
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axis of A perpendicular to the long axis of the molecule, we have 



(6.13 


and therefore (A t| - A^ = -3Aj ( P ^1 4, where P gives the orientational 
order of the long molecular axis and A ( ' is the component of the traceless 
hyperfine tensor parallel to its symmetry axis. The problem so far 
has been mainly an adaptation of known results, but of course in liquid 
crystals as compared to single crystals the interesting thing is 
that we have also anisotropic dynamic effects. To investigate the 
effect that the proper quantization of the nuclear spin has on the 
electron resonance spectrum we return to the time dependent spin 
hamiltonian in equation (6.4) ard rewrite it more explicitly as 


H(t) = H (0) + £(-) P g (2 ’ p) D 2 (t) D 2 ( Y ) T (2 ’ r) 

q,-p r,q g 

* K->’’ a<2 ' P> D q 2 ,- P <'> » 2 , q M T f ,r> . 


(6.14 


where is the scalar part of the hamiltonian. Now we decouple 

T a in its spin operator components (34) 


T A 2,r) = l C(112;n,r-n) S (1,n) I ( 1 » r_n ) f ( 6>15 

n 

and reject the non-secular terms. Then 

T A 2,r) = C(112;Or) S (1,0) I (1,r) . 


(6.16 
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We then transform I to its quantization frame which we indicate with 
an upper cap. 


’P ,r) = C(112, Or) S (1 ’ 0) l D 1 ( X ) 
A L m,r A 


(6.17 


m 


Substitution in equation (6.14) gives the truncated spin hamiltonian as 


H(t) = H (0) + [(-) P g (2,p) D 2 > _ p (t) D^ q ( Y ) T^ 2,r) <5 r() (6.18 


(1,0)1(1,m) 


+ £(-) P A (2,p) D 2 (t) D 2 ( Y ) D* ( X )C(112;0r)S (1,0) I 

u q,-p r,q m,r 


When we take an ensemble average the rotation of X causes only the 
m=0 term to survive in the 2nd line of equation (6.18) 


H - H W) * !(-)" S <2 ’ P) D 2 _ p D 2 0 <y ) T< 2 - 0) 


(6.19 


+ £(-) P A (2 ’ p) D 2 5 _ p D^ >0 (y) oJ >r ( X ) C(112;0r)S (1,0) I 


( 1 , 0 ) ( 1 , 0 ). 


The dynamic spin hamiltonian is 


H x (t) = H(t) - H 


It l M P 8 (2 ’ P) ( ^,. p (t)-D 2 ,Vo q J 


(2 ,m) 

g 


m p,q,r 


+ l [ l (-) P A (2,P) * (D 2 _„(t)-D 2 _6 n jD 2 „ (y) 


q»~P 0,-p Oq r,q 


m p,q,r 


C (112; Or) 1 -> (2,m) (6.20 

C (112 ;0m) m,r u;j A 
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where we have written 


T^ 2,m) = C(112,0m) S (1,0) 


as a tensor formally analogous to the truncated hyperfine tensor (6.16) 
With this definition we may now write the dynamic hamiltonian as 


H.(t) = yc-) m F (2 ’ m) T (2,m) . 

l L y w 


( 6.21 


m 


Since the |S m^|I m"> basiu is a good basis and since the hamilto¬ 
nian (6.21) is formally identical to that of the ordinary nitroxide 
problem (52) it follows that the linewidths are still given by the 
usual expression (6.1), with the coefficients A,B,C given by equations 
(4.62)-(4.64) in terms of the spectral densities . The tensors 

C 2 jjjN 

’ in equation (6.21) are given by 


if •- m) - J ( _)P g (2,p) (D 2^ (t) - D f p6oq} (r) Sn0 , <6.22 


p,q 


. (2,-m)_ (-) 


m 


C(112,0m) 


l (-) P A (2 ’ P) {D 2 (t)-D 2 6 } 

u q,~p 0,-p Oq 


P.q.r 


(6.23 


D 2 ( Y ) C(112,0r)D^ ( X ) 

r,q m,r 


This leads to the following spectral densities 


t 88 = V " (2,P> g (2,p) * (d 2 (y) } 2 j • 

0 ,q q-p 


4 s -1 • 


(6.24 



- 146 - 


8 (2>p) A < 2 - p> * d 2 >q (y) d 2 iq (y)C(n2,Or)dJ ir (x)i q . p <6.25 


jAA . 1 f A (2,p) A (2,p)*^2 , 


m 


C (112,mO) 2 


l A' z d^ q ( Y )d^, >q ( Y )C(112,Or)C(112,Or’) 


d^ (y) d^ i (y)j 

m,r * q-p 


(6.26 


or 


= _I_ 

“ C(112,mO) 2 


l A <2 - p) A ( 2 - p) * l j q _ p ( l d 2 >q < x )d“ >r ( x )C(112,Or)j' 


r—1 


(6.27 


It is easy to check that when there is no nuclear spin quantization 
effect, i.e. X=0 these spectral densities reduce to the ordinary ones 
(c.f. Chapter 4). In particular when the director is parallel or 
perpendicular to the field the spectrum is unchanged. The spectral 
densities (6.24-6.27) show also that another important feature of the 
angular dependent linewidth coefficients is unchanged i.e. they still 
possess the symmetry 


X( Y ) = X (tt —y), (6.28 

where X=A,B,C. This can be easily verified from the fact that 
X (y)“ - X( 7t- y) an d from the symmetry properties of the Wigner rotation 
matrices and the Clebsch-Gordan coefficients (34). In general it is 
no longer possible to write the linewidth coefficients in the form 
(6.2) as we did in the previous Chapter but, in view of the symmetry 
in equation (6.28), it is possible to expand the coefficients X(y) 
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in a series of even Legendre polynomials as 

X(y) = I \, P l/ COS ^ * ^ even (6.29 

Therefore the angular dependence in equation (6.2) can be regarded 
as a special case of equation (6.29) where the series terminates on the 
L = 4 term. 

6.3 Results 

To investigate the effect that the non-quantization of the nuclear 
spin along the magnetic field has on the linewidths we have performed a 
set of model calculations for a nitroxide spin probe with various order 
parameters. We have employed the same magnetic parameters as those used 
before (cf. page 118) for the cholestane spin probe whose structure is 
shown in figure 5.1. As far as the dynamic part of the problem is 
concerned, we have used the strong collision spectral densities given in 
Table 4.1, with the following correlation times: 


t o = 5. x 10 9 

s/rad 

T 2 = 1. X 10 -9 

s/rad 


In figures (6.2) - (6.4) we compare the angular coefficients A(y), 
B(y), C(y) obtained by allowing for the non-quantization effect with those 
obtained by neglecting this effect (x = 0). The order parameters employed 
were = 0.9, 0.8, 0.6, respectively, while the parameter was 

obtained assuming a Maier-Saupe pseudopotential (cf. Table 2.1). 

It is apparent that the effect is most important for the A and C 
coefficients as we had anticipated in Chapter 5. Therefore, in general, 
the order parameter P^ extracted from the B coefficients is likely to be 
more accurate than that obtained from the C, if the non-quantization effect 




Figure 6.2 The angular dependent linewidth coefficients A(y), E(y), C(y) 
calculated considering the non-quantization effect (continuous 

line) and neglecting it (dashed line). = 0.90, = 0.71, 

~9 -9 

xg = 5. x 10 s/rad, t 2 = !• x 10 s/rad. 
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90 


-1-5 


Figure 6.3 The angular dependent linewidth coefficients A(y), B(y), C(y) 
calculated considering the non-quantization effect (continuous 
line) and neglecting it (dashed line). = 0.79, 

_Q -9 

tq = 5. x 10 s/rad, t 2 = 1. x 10 s/rad. 


P, = 0.48, 
4 
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is neglected. 

The effect is clearly more important at high orders, since then the 
anisotropy in A is the largest. The problem of the influence that the 
non-quantization effect has on the correlation times, and on extracted 
by analysing experimental angular dependent linewidth coefficients is 
now under investigation. 
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CHAPTER SEVEN 

SPIN RELAXATION FOR 

BIRADICAL SPIN PROBES IN ANISOTROPIC ENVIRONMENTS 

7.1 Introduction 

The large zero-field splitting of certain biradicals makes them 
important candidates for spin probes of anisotropic systems such as 
liquid crystals and membranes. The electron resonance spectrum of a 
biradical dissolved in a liquid crystal may be influenced by the zero- 
field splitting in two quite distinct ways. Firstly the line positions 
are affected because the spin probe is partially oriented by the liquid 
crystal solvent. In addition the zero-field splitting coupled with the 
molecular reorientation constitutes a powerful spin relaxation process 
and so may determine the width of the spectral lines. Here we develop 
a theory of such line broadening, within the limit of fast motion, 
and show how it results in an angular dependence of the linewidths. 

The application of the theory is illustrated by studying the line 
broadening for a nitroxide biradical dissolved in a nematic mesophase. 

An analysis of the angular dependence of the linewidths suggests that 
the diffusion model provides a better account of molecular reorienta¬ 
tion in a liquid crystal than the strong-collision approach. Finally 
we draw attention to the potential value of the theory for understanding 
the linewidth variations in the deuteron magnetic resonance spectra 
of liquid crystals. 
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7.2 (a) Theory 

The majority of stable biradicals are formed by linking two 
nitroxide groups together; consequently the dominant magnetic interactions 
are the zero-field splitting together with the Zeeman and nitrogen 
hyperfine couplings. However, in developing the theory, we shall find 
it convenient to ignore all hyperfine interactions and such an assumption 
will be permissible provided the zero-field splitting exceeds the 
hyperfine interaction. In situations where the hyperfine terms are 
important it would be a straight forward matter to extend the present 
theory although the book-keeping could become somewhat involved (152). 

We begin by obtaining the static spin hamiltonian for a biradical in 
an anisotropic environment and then proceed to the spin relaxation 
problem. 

7.2 (b) Statics 

The spin hamiltonian for a given orientation of a biradical 
may be written, in irreducible spherical tensor notation (34), as 

H = H° + £(-l) P D (2,P) T (2,_p) , (7.1 

P 

where H° contains the orientation independent Zeeman and exchange 

/ n \ 

interactions. In the second term, ,P is the pth component of the 
zero-field splitting tensor and T^ 2 ’ p ^ is the appropriate combination 
of total electron spin operators. We list the components of the spatial 
and spin operators, which we use, in Table 7.1 because various definitions 
are encountered in the literature. 
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TAELE 7.1 The irreducible spherical tensor components 


P 

0 

±1 

±2 


D 


<2,p) 


(3/2) 


zz 


-(D ±iD ) 
+ xz yz 


(1/2)(D -D ±2iD ) 
xx yy xy 


t ( 2 , P ) 

(2/3)^{ S| - (1/4)(S + S_+S_S + )} 
;(1/2)(S ± S Z + S z S ± ) 

(1/2) S2 


This spin hamiltonian is written in a molecular co-ordinate system 
and so the spin operators fluctuate in time as the molecule reorients; 
it is this motion which combines with the zero-field splitting to 
yield the spin relaxation process. The time dependence may be removed 
from the spin operators by transforming to a space fixed co-ordinate 
system which we choose to contain the symmetry axes of the ordering 
potential responsible for the partial alignment of the spin probe. 

This gives 


H (t) ■ H° + l (-1) P D^ 2,P ^ T (2,q) d 2 _£t) (7.2 

p.q 

2 ..... 
where D (t) is a second rank Wigner rotation matrix which contains 
q, P 

the complete time dependence of the spin hamiltonian. The molecular 

motion is assumed to be sufficiently fast that the static spin 

hamiltonian may be obtained by taking a time or ensemble average of 

equation (7.2). In other words the relevant rotational correlation 

times are smaller than the inverse of the zero-field splitting. We 

~2 

therefore require the average and for an environment with 
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symmetry, such as a nematic or smectic A liquid crystal, this quantity 
vanishes unless q is zero. The static spin hamiltonian is then 


H = H° + D^ 2,0 ^ T^ 2, °\ 


(7.3 


where the partially averaged zero-field splitting is defined by 

D (2,0) - l <-l) P D (2 ’ P) D^_ p , (7.4 

P 

2 

and is a measure of the alignment of the spin probe by its 

0,-p 

environment. This spin hamiltonian corresponds to a triplet state 
with a cylindrically symmetric zero-field splitting tensor and the 
magnetic behaviour of such a system is completely understood (115). 

When the magnetic field is parallel to the symmetry axis of D the 
three triplet spin functions |l>, |0> and |-1> are eigenfunctions of H. 
For all other orientations of the sample these spin states are mixed 
by the zero-field splitting; however, provided D is small in comparison 
with the Zeeman splitting,then the triplet functions are good 
approximations to the eigenfunctions. Within this high-field 
approximation the two allowed transitions are 


|o> -*-*■ 11> = |+>, 


and 


(7.5 


) -1 > ■*-*■ 


1 ° > = !->• 
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The electron resonance spectrum therefore contains two lines with 
resonant fields 


B 4 = - — iL P 2 (cosY), (7.6 

gB 2g0 

where D 1( is the cartesian component of the partially averaged zero- 
field splitting tensor parallel to its symmetry axis and so to the 
director; the angle between the magnetic field and the director is y. 
The spacing between the two spectral lines changes dramatically with 
this angle and we shall now see that the linewidths should also 
exhibit a pronounced angular dependence. 

7.2 (c) Dynamics 

It is necessary to calculate the linewidths using the Liouville 
space or Binsch (112) formalism discussed in Chapter 4 because for 
certain orientations of the sample the separation between the two 
spectral components is comparable to the linewidth. Under such con¬ 
ditions Redfield's relaxation theory alone is unable to predict the 
line shape in contrast to the situation when the transitions are 
either degenerate or well resolved (100). Since the motion is assumed 
to be fast the line shape in a slow-passage, low-power electron 
resonance experiment is given by (c.f. Section 4.2) 

S(a>)« Re{S_* M _1 • S_} . (7.7 

The vector S_ is composed of the matrix elements <a|s_|a'> where 
|a>denotes an eigenfunction of H. The matrix M is related to 
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Redfield's relaxation matrix R (100) and the transition frequencies by 


M M. i i (u) | 

aa bb a a 


-w)6 


ab 


'b' 


+ R 


aa'bb' 


(7.8 


For the triplet-state problem the vector S_ is proportional to 
(1,1). The resonant frequencies may be replaced by resonant fields 
because, according to equation (7.6) there is a 1:1 correspondence 
between magnetic field and frequency. The relaxation matrix is 
calculated, in the usual way, from the dynamic spin hamiltonian H^(t) 
which is obtained by subtracting the static spin hamiltonian from its 
instaneous value H(t). Thus 


H x (t) 


l (~1) P D (2,P) (D 2 (t) - D 2 6 }T (2 ’ q) , 

L q»-p 0,-p Oq 

p,q 


(7.9 

where T' ,q depends on the orientation of the director with 
respect to the magnetic field. This angular dependence is removed 
by transforming the spin operator to a laboratory frame containing 
the field 


Hj (t) 


l (-D p d (2 ’ p) (d 2 (t) 

H 9 P 


P>q>r 


D 2 6 } D 2 T (2,r) 

q,-p Oq r,q 


(7.10 


which we write formally as 


r _(2,r) _(2,-r) 


H L (t) = l (-1) F U,r; T 


(7.11 


r 



LDD 


where F^ 2,r ^ contains both the time and angular dependence of the 
dynamic perturbation. The relaxation matrix is then found to be 


R 


< + | A 
<-| B 


-> 



(7.12 


where the basis functions correspond to the two allowed transitions: 
the matrix elements are 


A - -( 2 J 0 + 2 J 1 + J 2 ^ 


and 



(7.13 


(7.14 


These spectral densities are related to the Fourier transform of the 

(2,r) 

correlation function for the spatial operator F ’ by 


J = i 

r 


F( 2 » r ) (o) F (2,r) *(t) e ira)ot dt. (7.15 


The matrix M is then 


M H 


i (iog + <5-a)) + A 


i (a)o _ 6 _ m)+A 


(7.16 


within the same basis. The angular dependence of the spectrum is 
contained entirely within M for, as we have seen, the relaxation 
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matrix elements depend on the sample orientation and in addition 

6 = | D |( P 2 (cos y). (7.17 

The general expression for the absorption spectrum is found to be 

A-B+i (o)o - u)) 

S(o)) a -Re | (7.18 

[A+i (a>o~w) ] 2 +6 2- B 2 

This adopts a particularly simple form for three limiting situations. 
One occurs when cosy equals 1//3 for then the two transitions are 
degenerate and setting 6 equal to zero in equation (7.18) shows that 
this single line has a lorentzian shape width a width 

T~* = ~(A + B). (7.19 

At the other extreme the two lines are well resolved and so the 
line separation 6 is large compared with the linewidth whose 
magnitude is of order A. In addition inspection of equations (7.13) 
and (7.14) for A and B reveals that B is always less than A, conse¬ 
quently 6>B. The off-diagonal elements in M may therefore be 
ignored and the spectrum calculated from equation (7.16) is found to be 
a sum of two lorentzians separated by 6 and with equal widths: 

T" 1 = -A. (7.20 

For intermediate situations the entire line shape should be computed 
from equation (7.18); fortunately this is not always necessary. 
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When the molecular motion is slow in the sense that the rotational 
correlation time is large compared with the inverse of the microwave 
frequency idq then the non-secular spectral densities and are 
small in comparison with the secular contribution J^. Under these 
conditions B may be ignored in comparison with A and the line shape 
is lorentzian for all sample orientations with a width 



1 j 

2 O' 


(7.21 


In the remainder of the linewidth calculation for the electron 
resonance spectrum of the biradical we shall take this to be the case. 
The calculation demands a knowledge of the spatial operator 

(2 r) 

F ’ which is obtained by a direct comparison of equations (7.10) 
and (7.11). This operator may then be used to evaluate the spectral 
density J Q , for a particular orientation as 


J o ' 


I D< 2 '»> D 2 D 2 * , 

L 0,q 0,q 


P,P 

q,q’ 


(7.22 


i [< D > - D (0) D a* -D ,(t) ‘ D ' 2 

J q» p q > p 


. D , 6 ,}dt. 

0,-p 0,-p Oq Oq 


The simplification of this multiple summation is entirely analogous 
to the procedure adopted in calculating the angular dependence of the 
linewidth for doublet state species (c.f. Chapter 5). There it was 
shown that terms in the summation vanish unless q equals q’, because 
the mesophase is cylindrically symmetric, and so 



- 158 - 


J n = l (-l) P+P +q D (2,P) D (2,p,) * C(22L;00)C(22L;q-q)P. (cosy) 

U Li 

l;p.p' 

q (7.23 

X i I (D 2 (0) D 2 * ,(t) - D 2 D 2 , 6_ } dt, 

J q,-p q,-p o,-p o,-p Oq 


where the product of rotation matrices has been replaced with a sum 
of Legendre polynomials by resorting to the Clebsch-Gordan series (34). 


The summation is restricted to even values of L because of the D 
symmetry of the environment and so the angular dependence of the 
linewidth is found to be 


°°h 


T^ 1 = X Q + X 2 P 2 (cos y) + P^(cos y). 


(7.24 


To proceed further and obtain explicit expressions for the 
angular linewidth coefficients we shall require some model for the 
molecular reorientation process. Two quite different approaches 
(strong collision and diffusion) are available, although both find 
it particularly convenient to restrict attention to those spin probes 


for which the ordering tensor D Q is cylindrically symmetric, that 
is the only non-vanishing component has p=0. Of course the magnitude 
of the ordering tensor is determined by the solute-solvent interactions 


and we can only be certain that is cylindrically symmetric when 

the spin probe has the same symmetry. In practice there are few spin 
probes with cylindrical symmetry although this does not guarantee that 
the ordering tensor for the remainder are not cylindrically symmetric, 
but now the symmetry must be established experimentally. Providing the 


ordering tensor is cylindrically symmetric we may replace D_ , in 

u ,p 
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equation (7.23) for the spectral density J Q by P 2 6 Qp and, more impor¬ 
tantly, we can eliminate those terms with p not equal to p'; the 
expression for reduces to: 


l (~l) q |D (2,P) | 2 C(22L;00)C(22L;q-q) j q _ p , (7.25 
p»q 

where the spectral density j q _ p is 

j = i |{D 2 (0) D 2 * (t) - P 2 6. 6 } dt. (7.26 

J q-p J q,~p q,-p 2 Op Op 

— CD 

The components of the zero-field splitting tensor must now be expressed 
in the co-ordinate system which contains the symmetry axis for the 
ordering tensor as one of the axes. There are several variants of the 
first model which is based implicitly on the notion of a strong 
collision which results in instantaneous molecular reorientation through 
any angle (c.f. the previous Chapters). If, on average, the time 
between reorientational collisions is r then the correlation function, 
which is the integrand in equation (7.26), decays exponentially with 
a single time constant t; accordingly 


j q _ p "{ l (-D P+q C(22L’;q-q)C(22I/;-pp) P^ - P 2 6 Qp 6 Qq }T. 
L' 


(7.27 

As usual we introduce correlation times dependent on the 
modulus of p, t = t i i instead of x as a way to reflect the aniso- 

p Ip I 

tropy of the molecular reorientation (52). One of the major advantages 
of this approach is its mathematical simplicity, for example the angular 
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linewidth coefficients are readily calculated to be 

\ = f l (~1) P |D (2,P) | 2 C(22L;00)C(22L;-pp){P L -P 2 6 Qp }T p . 

P 

(7.28 

Evaluation of the Clebsch-Gordan coefficients (34) then gives the 
explicit expressions for X^ as 

X Q = (3/10){|d (2 ’ 0) | 2 (1-P 2 )t 0 + 2|d (2 ’ 1) | 2 t 1 + 2|d (2 ' 2) | 2 t 2 >, (7.29 

X 2 = (3/7){|d (2 ’ O) | 2 (P 2 -P 2 )t 0 + |d (2,1) | 2 P 2 t 1 - 2|D (2 ’ 2) | 2 P 2 t 2 }, 

(7.30 

and 

X 4 = (9/35){3|D (2 * 0) | 2 (P 4 -p2) To - 4 |d (2 * 1) | 2 ? 4 t 1 + |d (2 ’ 2) | 2 P 4 t 2 > . 

(7.31 

The linewidth coefficients thus depend on the three correlation 
times reduce the number of parameters we make the further 

assumption that 


T p 1= T 0 1 {4 + (T 0 T 2 X ' 1)p2} (7 * 32 

This relation provides a sort of interpolation for in terms of 
and x 9 and it is suggested by the comparison e.g. with the well know 
expressions for the correlation times for rotational diffusion in an 
isotropic system (126). Of course the validity of the assumption can 
only be justified by comparison with experiment and at present the few 
results available seem to support these intuitive expression for the 
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spectral densities in anisotropic systems. 

As we mentioned in the Introduction, electron spin relaxation 
for a triplet state is identical to nuclear spin relaxation, for a 
nucleus of spin 1, by the quadrupole interaction (153). Consequently 
a partial check of equations (7.29-7.31) is possible because Egozy 
et al (154) have calculated the deuterium relaxation times, and 
for perdeuterated benzene dissolved in a nematic mesophase. 
However, they restrict their attention to the situation when the 
director is parallel to the magnetic field, the field gradient at the 
deuterium nucleus is cylindrical symmetric and the molecular reorien¬ 
tation can be described by a single correlation time, as in the simple 
version of the strong collisional model. If we impose all of these 
restrictions on equations (7.24) and (7.29-7.31) then we obtain the 
secular contribution to the linewidth, T^which has the same form 
as that found by Egozy et al but which is twice as large. Since our 
results yield the correct limiting value for the secular part of 
when both the system and the motion is isotropic (153) we suppose their 
equation is error. At the other extreme the molecular reorientation 
is assumed to proceed via small angle jumps, that is according to the 
diffusion model discussed in Section 4.4 ( b ). According to this 
model the spectral densities at zero frequency are given by (96) 



l {(-l) P+q C(2L , L ,, ;q-q)C(2L'L ,, ;p-p)P L „ - Vl'V 0q >T L 
L'*0 
L"(even) 


»q.p‘ 


(7.33 
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The time constants x^, depend on the components of the probe diffu¬ 
sion tensor and on the orientational order; they have been are defined 
in terms of the Nordio et al (96) R qp matrix in Section 4.4 (d) . 

The diffusion tensor D is assumed to be cylindrically symmetric and its 
symmetry axis to be coincident with that of the ordering tensor. 

An explicit expression for the linewidth coefficients can be obtained, 
remembering the symmetry rules 


as 


and 



" J -q-p J q-p J -qp ’ 


<3/10){|D <2,0) | 2 (i 00 *2j 10+ 2j 20 ).2 

|d <2 * 1) | 2 (i 02 + 2j u + 2j 

* 2|D (2 ’ 2) | 2 (j 02 *2j 12 .2j 22 )), 

(7.34 

(3/7){|D (2 ' 0) | 2 <j 00 +j 10 -2j 20 )+2|D 

<2> ‘ > l 2 <W«2i> 

* 2|D (2 ’ 2) | 2 (j 02+ j 12 -2j 22 )>, 

(7.35 

(9/35)1|D (2 - 0, | 2 (3j 00 -4j 10+ j 20 > +2 

l° (2 - 1> l 2 (3j 0 r 4 in*j 

♦ 2|D <2 ’ 2) | 2 (3j 02 -4j 12 *j 22 )). 

(7.36 


In the following section we shall compare the predictions of the 
two dynamical models with the angular linewidth coefficients measured 
for a nitroxide biradical dissolved in a nematic mesophase. 



7.3 (a) The Experiment 


The biradical spin probe employed to test the theoretical 
predictions was 1,4-bis(3-spiro-[2'N-oxyl-3',3'dimethyloxazolidineJ ) 
cyclohexane which has the structure shown in figure 7.1. The spin 
probe was synthesised from cyclohexan-1,4-dione using the route 
described by Michon and Rassat (155). The nematogen, 4,4'-di-n-buty- 
loxyazoxybenzene, was purchased from Eastman Kodak Ltd. and purified 
by recrystallisation from ethanol (Tg_ N 102°C; T n _j 132°C). In the 
experiment the orientation of the director is controlled by an electric 
field, consequently the sample, made by doping the nematogen with the 
spin probe, was contained in a cell made from two glass plates approxi¬ 
mately lcm*10cm and held 125 p apart with a Teflon spacer. The 
inner surfaces of the plates were coated with tin oxide and this allows 
the application of an electric field perpendicular to the surface. 

In the absence of an electric field the director is aligned parallel 
to the spectrometer's magnetic field, however the magnetic forces 
responsible for this alignment are overcome by a 14 kHz electric field 
from an applied voltage of 90 volts. Under these conditions the 
director is aligned perpendicular to the electric field and so parallel 
to the plates of the cell because the anisotropic permittivity is 
negative (39). The only effect of the magnetic field now is to constrain 
the director to lie in the plane containing both the magnetic and 
electric fields. The angle y between the director and the magnetic field 
is then 90° minus the angle made by the electric field with the magnetic 
field. The cell was mounted in a goniometer and the angular dependence 
of the spectrum for the spin probe in the nematic mesophase was 
recorded; the measurements were made at two temperature, 127°C and 104°C. 



o 


o 



o 





Figure 7.1 The molecular structure of the biradical spin probe 
1,4-bis(3-spiro-f2' N-oxyl-3',3'-dimethoxazolidine] ) 
cyclohexane together with a co-ordinate system in which 
the ordering tensor is supposed to be cylindrically 
symmetric. 
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Typical spectra are shown in figure 7.2 for the sample at 104°C with 
the director parallel and perpendicular to the magnetic field, these 
two spectra clearly demonstrate the angular dependence of the 
linewidth and we shall now describe the determination of these widths. 

7.3 (b) Spectral Analysis 

The spectra in figure 7.2 are dominated by the dipolar splitting 
whose angular dependence is predicted by equation (7.6); for those 
orientations where the dipolar splitting is well resolved there is 
complete agreement between the observed and theoretical line separa¬ 
tions. The values of determined from the angular dependence for 
the two temperatures are given in Table 7.2. The spectrum, for the 

TABLE 7.2 The partially averaged magnetic parameters 







Temperature 

D l< 

A .l 

A 1 

D ,|/(A„-A_) 

(°C) 

(MHz) 

(MHz) 

(MHz) 



104° 

101 

29.2 

47.2 

-5.6 

127° 

64 

33.7 

44.9 

-5.7 


director orthogonal to the magnetic field, also exhibits hyperfine 
structure from the two equivalent nitrogen nuclei, the spacing 
between the lines is just one half the nitrogen coupling constant 
because the electron-electron interaction greatly exceeds the hyperfine 
interaction (156) . For many orientations the hyperfine splitting is 
not resolved but does result in inhomogeneous broadening of the spectral 





I I 100 gauss 

gure 7.2 The electron resonance spectrum of the biradical in the 
nematic mesophase of 4,4'-di-n-butyloxyazoxybenzene at 
104°Cwith the director (a) parallel and (b) perpendicular 
to the magnetic field. 
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lines, consequently it is essential to allow for this broadening when 
obtaining the true linewidth. We have therefore simulated spectra 
which incorporate both the dipolar and hyperfine splittings, the 
linewidth was then varied until agreement with the observed spectrum 
was obtained. The angular dependence of the coupling constant, required 
in the simulations, is given by (39) 


a(y) = {A 2 


cos 2 y 


A 2 sin 2 y) ^, 


(7.37 


since the anisotropy in the partially averaged g tensor is negligibly 
small. The values of the two components of the partially averaged 
nitrogen hyperfine tensor employed in the spectral simulations are 
given in Table 7.2. Since a(y) could not always be obtained from 
the spectrum of the biradical we have also measured the nitrogen 
coupling constant of the nitroxide monoradical, which is responsible 
for the central three lines in figure 7.2a. The radical is almost 
certainly the hydroxylam?ne derivative of the biradical formed by 
hydrogen abstraction from the solvent. In view of the great similarity 
of the structures of mono and biradical we expect their partially 
averaged nitrogen hyperfine tensors to be the same and this is found 
to be so when both splittings can be measured simultaneously. The 
linewidths obtained by matching the theoretical and experimental spectra 
are plotted against the director orientation in figure 7.3 for the 
sample at 104°C; a similar angular dependence was observed at 127°C. 



60 


40 

t;7mhz 

20 


Figure 7.3 




166 - 


7.3 (c) The Molecular Symmetry 

The expressions obtained for the angular linewidth coefficients 


are only valid when the ordering tensor D for the spin probe is 

u, p 

cylindrically symmetric. The spin probe employed in the experiments 
does not possess cylindrical symmetry and so we must first show that 


Dq o is the only non-zero component of the ordering tensor. There 
are five unknown order parameters and since there are only two 
independent pieces of experimental information (A^ and D^) all we 
can hope to achieve is to show that they are consistent with an assumed 
cylindrical symmetry. The spin probe may exist in one of several 
conformations but solid-state electron resonance studies indicate that 
the equatorial-equatorial chair form of the trans isomer, shown in 
figure 7.1, is preferred (157). Examination of molecular models then 
suggests that the ordering tensor for the spin probe in the nematic 
mesophase should be cylindrically symmetric about the z axis shown in 
figure 7.1, together with a projection formula for the conformer. 
Provided the assumption of cylindrical symmetry is correct then the 
components of the partially averaged tensors are given by 


and 


D (| = (2/3)* D (2 ’ 0) P 2 , 


A (| - A ± = (3/2)* A (2 ’ 0) P 2 , 


(7.38 


(7.39 


where the zeroth component of the appropriate anisotropic tensor is 
evaluated in the xyz co-ordinate system, shown in figure 7.1. The 
nitrogen hyperfine tensor, for such oxazolidines, is cylindrically 
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symmetric about the x axis (157) and so 
A (2 ’ 0) = -(3/8)* A] ( 


(7.AO 


where A^ is the component of the anisotropic tensor along this symmetry 
axis. The zero-field splitting tensor is, to a good approximation, 
also cylindrically symmetric (155,157) but about the z' axis which 
makes an angle 0 with the z axis and so 


D 


( 2 , 0 ) 


(3/2) 


i (3 cos 2 9-1) 


7 I > • 

z'z' 


(7.41 


Because the symmetry axes of the two magnetic tensors are not parallel 
the ratio D || /(A || - A x ) may only be independent of temperature if the 
ordering tensor is cylindrically symmetric. The values calculated 
for this ratio are essentially constant as we can see from the results 
given in Table 7.2. The angle 0 can be extracted from these values 
for D (| /(A^-A ; ) since, from equations (7.38-7.41) 


(3cos 2 6-1) 

2 



4(A„ -A x ) 



(7.42 


Calculation of 0 provides us with an additional check on the 
assumption of cylindrical symmetry because this angle may be estimated 
from molecular models. Thus the z' axis has been found to intersect 
the N-0 bonds at about one quarter of the bond length from nitrogen 
(157) and so 0 should be about 20°, given our assignment of the z axis. 
Since 0 is relatively small rather precise values of the parameters 
occurring in equation (7.42) are required to obtain a reliable estimate 
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of this angle. The zero-field splitting tensor is known, with some 
accuracy, from single crystal studies and D^,^, is found to be 230 MHz 
(157). Unfortunately the anisotropic hyperfine tensor is not known 
with such certainty and three reasonable values are available. The 
single crystal study of the biradical yields a value for of 52.5MHz 
although the accuracy is said not to be high (157). A more reliable 
value of 51.1 MHz has been determined for the structurally related 
monoradical 2-N-oxy1,1,1,3,3,-tetramethyloxazolidine (157). Finally 
A' is 48.7 MHz for another monoradical, (3~spiro~[2'N-oxyl-3',3'- 
dimethyloxazolidine])-5a~androstane-17g-ol (158), which, in view of 
the cyclic groups, should be most closely related to the biradical. 

The magnitudes of 0 and obtained from the single value of D^,^, 
together with these three values of A' are listed in Table 7.3 for 
both temperatures. The angle 0 is seen to range from 8° to 15° even 
though the component of the anisotropic hyperfine tensor only changes 
by an amount, 3.8 MHz, not greatly in excess of the experimental error 

in A' . 

It 

TABLE 7.3 The derived parameters 0 and 


A* /MHz 

Temperature (°C) 

►—* 
o 

O 

127* 

3 

e 

P 2 

0 

m 

52.5 

9° 

0.46 

8° 

0.29 

51.1 

12° 

0.47 

11° 

0.29 

48.7 

15° 

0.49 

15° 

0.31 
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In view of this extreme sensitivity of the angle to A’ ( we believe 
the results for 6 are also in accord with the assumption of a cylin- 
drically symmetric ordering tensor. In addition we shall take 6 to 
be 15° in subsequent calculations because this value is closest to 
our original estimate of 20° based on an examination of molecular 
models. 


7.3 (d) The Angular Linewidth Coefficients 

As we have just seen our assumption of a cylindrically symmetric 
ordering tensor is not unreasonable and so we may now turn to a detaile 
analysis of the linewidths. The first test of the theory is to see if 
the angular dependence of the linewidths is in accord with that predict! 
by equation (7.24). The theoretical angular dependence of T 2 * is 
shown as the solid line in figure 7.3 and is seen to be in good agreemei 
with experiment; comparable agreement is also found for the linewidths 
measured at 127°C. Of course such comparisons do not provide a parti¬ 
cularly severe test for the theory because the predicted angular 

dependence of the linewidth stems directly from the D , symmetry of the 

°°h 

anisotropic environment and the second rank nature of the zero-field 
splitting tensor. In fact the good agreement between theory and 
experiment may be taken as confirmation of our correct assignment of 
the dominant spin relaxation process. A more searching test of the 
theory is provided by an attempt to predict the observed angular line- 
width coefficients X ; these are listed in Table 7.4 and were obtained 

Li 

by fitting equation (7.24) for T 2 ^(y) to the experimental linewidths. 
There are, of course, two theoretical expressions for the coefficients 


and one of the major objectives of any experimental linewidth 
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TABLE 7.4 

The angular 

linewidth 

coefficients 



Temperature X Q 

X 2 

X 4 

X 0 /X 2 

X 0 /X 4 

(°C) 

(MHz) 

(MHz) 

(MHz) 



104° 

31.4±0.5 

16±1 

-26±1 

2.0±0.2 

-1.21±0.07 

127° 

29.5±0.9 

9±2 

-14±2 

3.3±0.8 

-2.1 ±0.4 


investigation is to distinguish between the two extreme models for 
molecular reorientation. However in this study we shall be primarily 
concerned with using the experimental results, which are not of high 
accuracy, simply to illustrate the theory. 

We now proceed to calculate the angular linewidth coefficients 
using both models. To make the comparison with experiment as equitable 
as possible the ordering potential, given in general by a sum of even 
Legendre polynomials (c.f. Chapter 2), will be restricted to just the 
first term, i.e. 


IcT = *2 P 2^ COS6 ^ * (7.43 

in both calculations. The magnitude of the single parameter X^ may 
then calculated from the observed value of the order parameter ; 
the results of this calculation are in Table 7.5. Because the potential 
is restricted to a single term both dynamic models lead to sets of 
equations for X^ which contain just two arbitary parameters. In addi¬ 
tion the ratios X^/X^ and X^/X^ of the angular linewidth coefficients 







171 - 


TABLE 7.5 The order parameters 


Temperature 

(°C) 

P ~2 


A 2 

104° 

0.49 

0.16 

-2.3 

127° 

0.31 

0.06 

-1.4 


are seen to depend on the ratio of the appropriate parameters; in 
the case of the strong collision model this is while for the 

diffusion model it is D (| /D^. 

We begin with the strong collision model; the evaluation of the 
linewidth coefficients requires the order parameter and this was 
calculated from the values for with the results given in Table 7.5. 
This procedure contrasts with previous linewidth studies employing 
doublet states where the angular linewidth coefficients were employed 
to determine P^ (52-54). The components of the zero-field splitting 
tensor in the principal co-ordinate system, xyz, for the ordering 
tensor were found to be 


d ( 2 »°) = 252 MHz, 

D (2,iD = g5 MHz ^ 

(2 + 2 ) 

and D v * ' = 12 MHz. 


(7.44 


The two ratios of the linewidth coefficients, X^/X^ and X^/X^, were 
calculated from equations (7.29-7.31) for a wide range of values of 
the relative correlation time t q / T 2 an attempt to find agreement 
with experiment. However both Xq/^ and X Q /X^ proved to be virtually 
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independent of t^/ and it was not possible to fit the experimental 
results. For example at 104°C, where the experimental data is most 
reliable, X^,^ was predicted to change only from 2.3 to 2.2 when 
Tq / t 2 was varied from 1 to 8. This prediction is in reasonable 
agreement with experiment, but, in marked contrast, X Q /X^ is calcu¬ 
lated to be -3.7 which is over three times the observed value. It is, 
of course, possible to force agreement with the experimental value 
of Xq/X^ by changing the magnitude of the order parameter P^, since 
this affects just the X^ angular linewidth coefficient. Agreement 
can only be achieved when P^ is virtually zero; a result which is in 
accord with small values of P^ determined by Raman studies (51). 

However we suspect that the value is unreasonably small in view of 

the success of the single parameter potential in linewidth investigations 

employing doublet state spin probes (52-54). 

The analysis of the angular linewidth coefficients obtained for the 
biradical at 127°C yielded similar discrepancies between theory and 
experiment. 

The angular linewidth coefficients were also evaluated using 
the diffusion model; again the ordering potential was restricted to 
a single term A ^ taking the values given in Table 7.5. For these 
relatively low values of Ait was found to be adequate to truncate 
the series expansion of the spectral density in equation (7.33) 
after L"= 20. The relative linewidth coefficients were obtained for 
various values of D |( /D, and although the results were not particularly 
sensitive to D/Dj it did prove possible to obtain the ratios X^/X^ 
and X Q /X^ in fair, but not complete, agreement with experiment. 

The best fit was found for a range of values for D (| /D and these are 

II 



listed in Table 7.6, together with the predicted linewidth ratios. 


TABLE 7.6 The diffusion model predictions 


Temperature 

(°C) 

V D i 

V X 2 

V X 4 

D, 

•L 

(Mrads 

• 

1.2 

1.8 

- 1.4 

600 

104 ° 

1.4 

1.8 

- 1.5 

520 


1.6 

1.8 

- 1.6 

540 


1.2 

2.3 

- 3.4 

720 

127 ° 

1.4 

2.3 

- 3.6 

700 


1.6 

2.3 

- 3.8 

680 


As we expect the agreement between theory and experiment is best for 
the more accurate angular linewidth coefficients measured at 104°C. 

The agreement for the less accurate results at 127° is not so impressive 
although still an improvement on the predictions of the strong colli- 
sional model. The ratio D^/D for the components of the diffusion 
tensor is independent of temperature as we might anticipate (53). 

We can estimate this ratio in the hydrodynamic limit if the molecular 
dimensions are known (see Section5.5 ); using molecular models we 
gauge the major semi axis to be 12 & while the minor semi axis is 
about 4.5 X. The hydrodynamic value of D^/D, is then found to be 
2.7, which is about twice as large as the experimental value. It is 
surprising to find a value below the hydrodynamic limit bacause similar 
studies of doublet-state probes have yielded value in excess of this 
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limit (53,128,141). Indeed, for these systems it proved necessary 
to introduce an anisotropic interaction tensor k (53) whose magnitude 
is related to the anisotropy in the intermolecular potential (146) 
to account for the deviations from the hydrodynamic limit. Of course 
analogous explanations cannot be invoked, for rod-like molecules, 
when the ratio D^/Dj is less than the hydrodynamic limit (53). 

The diffusion model of molecular reorientation in an anisotropic 
medium provides an adequate description of the angular linewidth 
coefficients, even though the ordering potential is restricted to a 
single term. This limited success contrasts with the apparent failure 
of the strong collision model which is only able to predict the ratio 
Xq/Xz when the ordering potential contains a single term. In the 
anisotropic system which we have studied the molecular reorientation 
would appear to occur via rotational diffusion although more results 
will be required before we can differentiate with complete certainty 
between the two models. 

7.4 Additional Spectral Features 

When the director is orthogonal to the magnetic field the five 
nitrogen hyperfine lines associated with either component of the 
dipolar doublet are well resolved. Close inspection of the appropriate 
spectrum, shown in figure 7.2, reveals a minor variation in the widths 
of the hyperfine lines. Such asymmetric line broadening cannot be 
accounted for by the theory which we have just developed because the 
anisotropic hyperfine interaction was specifically excluded from the 
dynamic perturbation. However quite straightforward arguments (159) 
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show that this line broadening can be explained by inclusion of the 
hyperfine term in H^(t), the asymmetric linewidth effect then originates 
from the cross term between the anisotropic hyperfine interaction and 
the zero-field splitting tensor. We do not feel impelled to make 
detailed analysis because this broadening is quite secondary to the 
other linewidth variation which we have observed. 

The spectral line shape is expected to be symmetric about the 
base line but examination of the spectrum for the director parallel 
to the magnetic field shows that this is not the case. Similar 
asymmetry in the line shape has been observed for doublet state spin 
probes dissolved in liquid crystals and has been attributed to non- 
uniform alignment of the director, probably caused by thermal fluctua¬ 
tions (149). A slight departure from uniform alignment of the director 
would also be expected to produce asymmetric line shapes for biradical 
spin probes and we have used spectral simulations to confirm this 
expectation. Of course we also expect the line shape to become asymmetric 
when the rate of molecular rotation is comparable to the inverse of the 
zero-field splitting (108). Unfortunately it is difficult to devise 
experiments to distinguish between these two possible explanations of 
asymmetric line shapes. However the rate of motion may be obtained 
from the inverse of the diffusion tensor; from the values listed in 
Table 7.6 we can see that the ratio of the zero-field splitting to 
either component of the diffusion tensor is just less than unity. 
Accordingly the assumption of rapid motion should be valid for this 
system and the asymmetric line shape presumably originates for non- 
uniform alignment of the director. 
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7.5 Quadrupole Relaxation 

The problem which we have just solved is formally identical to 
the calculation of the spin relaxation times, for a nucleus of spin 1, 
caused by a quadrupole interaction coupled to the molecular reorienta¬ 
tion. Consequently our results could be of value in understanding 
the spin relaxation behaviour of nuclei, such as deuterium, when they 
are incorporated in molecules forming part of a liquid crystalline 
system. However some caution must be exercised in the application of 
these results. For example although the neglect of the non-secular 
terms in the dynamic perturbation is a good approximation for electron 
resonance spectroscopy such an assumption may not be valid for nuclear 
magnetic resonance where the operating frequency is much smaller. 

If the non-secular terms do have to be retained then the entire line- 
shape will have to be calculated from equation (7.18) and a simple 
analysis is only possible when the lines are well resolved or the 
transitions are completely degenerate. In addition our assignment 
of the dominant process responsible for changes in the molecular 
orientation may be quite irrelevant for nuclear magnetic resonance. 

This situation obtains because the quadrupole interaction is consi¬ 
derably smaller than the zero-field splitting and so the motion with 
respect to the director, which we have considered, is too fast to make 
a significant contribution to the nuclear spin relaxation times (160 )• 
However, the molecular orientation is also changed by thermal fluctua¬ 
tions in the director and although this motion is slow on the electron 
resonance time scale it is sufficiently fast to influence the relaxation 
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times in nuclear magnetic resonance (161). Consequently the spectral 
densities employed in this Chapter would be inappropriate for calcula¬ 
ting the deuteron relaxation times; it is however straightforward to 
obtain the relevant spectral densities from continuum theory (161). 

The calculation of the spin lattice relaxation time also presents 
certain difficulties especially when the two transitions are nearly 
degenerate. Indeed the general angular dependence of is complicated 
even when the spectral lines are well resolved and we shall not discuss 
this aspect of the problem here. 
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APPENDIX 1 


IRREDUCIBLE TENSORS AND WIGNER ROTATION MATRICES 


A tensor of rank n is a quantity that transforms under rotation 

t h n 

as the n direct power of a vector. The 3 dimensional representation 

of the rotation group realized in this way can be decomposed in a set 

of irreducible representations D^ each of dimension (2L+1). If we 

describe the rotations in terms of Euler angles (aBy) then the matrix 

2 

elements of the irreducible representation in a basis where J , J 
(i.e. the angular momentum and its projection) are diagonal are (34) 

-ictJ -ifSJ -iyJ 

D (aSy) = <Lm |e Z e y e Z |Ln>. (Al.l 

m,n 

The matrix elements D^ (aBy) are called Wigner rotation matrices 

m,n 

or Wigner functions or generalized spherical harmonics. Combinations 
of ordinary tensor components transforming according to the representa¬ 
tion D^ are called irreducible tensor components of rank L: T^ ,m ^ 

e.g. T (L,m)' = £ T (b,n) (A1.2 

u n,m 
n 

Equation (A1.2) illustrates the main reason of the usefulness of 
irreducible tensors in problems involving rotations i.e. that their 
transformation properties are very simple. The set of (2L+1) compo¬ 
nents, T^\ is called an irreducible tensor of rank L. The Euler 
angles (aBy) in equation (A1.2) determine the rotations which carry 

the original (unprimed) coordinate system into the rotated (primed) 
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one. We follow the convention of Rose (34) for the Euler angles, 
which is illustrated in figure Al.l. From equation (Al.l) it is 
apparent that we can express ^(ciBy) as (34) 


L , „ s -imoi _L -my 

D (aBy) = e d (6) e 
m,n m,n 


(A1.3 


where the real quantities 


-iBJ 

d (B) = <Lm|e ^ |Ln> 
m, n ' ' 


(A1.4 


are called reduced Wigner matrices. 

The functions D L (aBy) constitute a complete orthogonal set 
m, n 

spanning the space of the angles a,S»Y* Some of their properties, 
which we use most often in this thesis, are: 

1. Orthogonality 


[ D L (aBy) D L , , (aBy) dasinBdBdy = 

I m,n m ,n 2L+1 m,m n,n L, 

2. Integral of three Wigner rotation matrices 


L’ 

(A1.5 


D 


L"* 


•• „n (aBy)D' 


m ,n 


T » T 

, 1 (aBy)D (aBy)dasinBdBdy 

m ,n m,n 


8tt 2 , 

2L 1, +1 m+m' ,m" 


6 . „C(LL'L";mm')C(LL'L";nn') 

n+n ,n 


(A1.6 


where C(abc;de) is a Clebsch-Gordan coefficient (34). 
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3. Closure 


1 D »,n (a l 6 l Y l ) D n,o' VaV ' D m,n>' ( “ 6r) > 


(A1.7 


n 


where («8 y) is the resultant of (cjjBjYj) then ( a 2 B 2 Y 2^‘ 
4. Symme try 


D L * (oBy) = (-) m_n D L (oBy) = D L (-y-B-o) . 

m,n -m,-n n,m 


(A1.8 


5. Products 


,, T „ L'+L" 

d , , d „ „ = y 

m ,n m ,n L 




C(L , L"L; m 'm ,, )C(L , L M L;n , n")D L , „ , 

m +m ,n +n 

(A1.9 


6. Special cases 


D m,0 (a60) ( 2L+1 ^ Y L,m (B » a) ’ 


(A1.10 


D o,o (OBO) = d o,o (e) = p l ( cosP) ’ 


(Al.ll 


where Y is a spherical harmonic and P a Legendre polynomial. 

L|I11 jL 
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APPENDIX 2 


SOME PROPERTIES OF THE RUMMER CONFLUENT HYPERGEOMETRIC FUNCTIONS 


We shall give here only the properties used in the text. 

A thorough discussion of the confluent hypergeometric function is given 
in References (79,80)• The Rummer function M(a,b,z) can be defined as 


(a) m 

M(a,b,z) = l —- J , (|z|<-, b*0, -1,-2,...), 


m=0 ^ 


(A2.1 


where z,a,b, can take real and complex values and 


(a) * a(a+l)....(a+m-1); (a) = 1. 

m U 


A useful integral representation is 


M(a,b,z) = 


T(b) 


r(b-a)r(a) 


exp(zt) t 3 1 (l-t) 15 3 ' dt. 


(A2.2 


with Re b >Re a >0. 

Differentiating equation (A2.1) n times we find 


,n (a) 

— M(a,b,z) =-— M(a+n,b+n,z) 

dz 11 (b) 

n 


(A2.3 


A general asymptotic expression for large arguments is 


(b-a) (1-a) 

M(a,b,z) = e Z z (b a) [ - — - z m +0(|z| n ’) , 

m=0 


(A2.4 
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where 

|arg z|<ir/2 -6; a,b#0, -1, -2,... 


The Kummer functions can be expanded in a series of modified Bessel 
functions 


M(a,b,z) 


z/2 

e 


r(b-a-l/2) (z/A) a b+1/2 


oo 


X l 

n=0 


(2b-2a-l) 

n 


(b-2a) (b-a-l/2+n) 

n 

n! (b) 

n 


(-11 

^ J b-a-l/2+n 


(A2.5 


with 


b 4 0, 


(z/2). 


- 1 , - 2 ,.... 
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